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development  of  a  procedure  capable  of  making  full  use  of  such  censored  informalion  would  seem  very  timely  and  neces¬ 
sary 

In  section  I  of  this  report,  we  present  a  maximum-likelihood  regression  scheme,  “MLE-CY",  which  takes  all  the 
censored  yields  into  account  to  refine  the  estimated  m„:yield  relationship.  This  regression  routine  is  very  similar  to  the 
maximum-likelihood  estimator  used  in  computing  the  optimal  network  mb  values  based  on  the  censored  station  ampli¬ 
tude  measurements  due  to  clipping  and  to  non-detection.  In  the  non-censored  case,  it  gives  results  identical  to  those 
derived  by  the  standard  least-squares  method.  Applications  of  this  scheme  to  the  explosions  from  several  test  sites  of 
different  geology  show  that  it  is  a  superior  procedure,  as  compared  to  the  conventional  least-squares  approach.  The 
same  algorithm  can  be  applied  to  other  magnitude  measurements  such  as  Ms,  Pcoda,  m^LJ,  M0,  RMS  Lg  and  DOB 
etc. 

We  have  aiso  conducted  a  systematic  analysis  of  the  magnitude.yietd  relationship  at  five  major  test  sites  using 
miscellaneous  unclassified  magnitudes.  (A  classified  annex  using  the  official  mb  values  will  be  furnished  separately.) 

Several  noteworthy  results  are  summarized  here: 

[1]  Including  the  censored  yields  in  the  regression  does  improve  the  accuracy  of  the  estimates.  In  reality,  both  the 
magnitude  and  the  yield  measurements  are  subject  to  error.  Pending  the  determination  as  to  which  of  the  two 
extreme  hypotheses,  namely  o(mb)la( Y)  «  0  and  o(mb)/cr(Y)  =  is  closer  to  the  real  situation,  we  also  included 
the  results  based  on  Ericsson's  method  with  various  a(mb  )/cr(Y).  As  expected,  we  can  see  the  smooth  transition  of 
estimated  parameters  (i.e. ,  the  slope  and  the  intercept)  as  a(mb)/o(Y)  varies.  Thus  the  censored  cases  with  non¬ 
trivial  o(mb)/o(Y)  values  oould  be  “interpolated".  Our  maximum-likelihood  regression  scheme  and  Ericsson’s 
method  represent  two  different  directions  in  extending  the  standard  least  squares. 

(2)  For  Shagan  events,  Ringdal's  RMS  Lg  provides  the  smallest  scatter  around  the  calibration  curve,  provided  that 
low-yield  events  with  mb(RMS  Lg)  <  5.5  or  yield  <  40KT  are  excluded.  Geotech's  GLM  method  gives  network  mb 
values  better  than  almost  all  other  magnitudes  based  on  the  teleseismic  P  waves  and  log^J ,  in  terms  of  both 
the  yield  estimation  and  the  mb  scaling  against  Ringdal’s  RMS  Lg.  For  all  five  test  sites  we  have  compared,  mb 
measurements  reported  by  ISC  and  NEIS  are  Diased  high  systematically  at  low  yields. 

[3J  A  direct  estimation  of  the  test  site  bias  suggests  that  Nuttli's  (1987,  1988)  Degelen  puzzle  could  be  invalid  simply 
because  ot  the  relatively  poorer  quality  mb(ISC)  used.  Our  data  indicate  that  Shagan  River  Test  Site  is  more 
efficient  in  exciting  teleseismic  P  waves  than  Degelen  Mountain,  consistent  with  our  previous  modeling  study.  Also, 
the  test  site  bias  is  yield  dependent,  in  agreement  with  other  observational  study. 

[4]  We  present  an  alternative  approach  to  derive  the  mb  adjustment  converting  cratering  shots  to  contained  explosions 
of  the  same  yield.  The  correction  derived  by  this  approach  seems  to  match  that  by  the  multichannel  deconvolution 
method  rather  well. 

(5j  Degelen  Mountain  is  the  only  test  site  that  has  a  decreasing  fogiP^JP,}  and  log (Pb/P,)  with  increasing  yields.  It 
is  also  the  only  test  site  for  which  the  phase  “a"  (i.e.,  zero-crossing  to  first  peak)  shows  the  smallest  scatter 
around  the  calibration  curve,  as  compared  to  the  phases  “b”  (i.e. ,  first  peak  to  first  trough)  and  “max”  (i.e. ,  max 
peak-to-trough  or  trough-to-peak  in  the  first  5  seconds).  Both  the  mountainous  topography  (which  causes  complex 
pP  interference)  as  well  as  the  testing  practice  (e.g. ,  the  relatively  shallow  and  abnormal  shot  depths)  oould  be 
responsible.  At  Shagan  River,  the  phase  “b”  has  the  smallest  scatter  around  the  calibration  curve.  These  observa¬ 
tions  confirm  the  conjecture  (DARPA,  1981)  that  in  a  proper  environment  the  first  cycle  could  give  better  results 
than  does  “max"  phase. 

[6]  The  scale  depth  for  Konystan  explosions  is  146±1  meters,  and  the  depth  of  burial  [DOB]  is  roughly  proportional  to 
the  quartic  root  of  the  yield,  rather  than  the  cubic  root  as  frequently  cited  at  NTS.  This  empirical  scaling  rule  is 
applicable  to  Shagan  River  region,  but  not  Degelen  Mountain.  For  Konystan  and  Shagan  regions,  the  yields 
estimated  using  depth  scaling  have  accuracy  comparable  to  those  using  mb. 
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EXECUTIVE  SUMMARY 


Conventional  methods  for  estimating  underground  explosion  yields  from  seismic  record¬ 
ings  are  based  on  the  use  of  some  appropriate  “magnitude:yield”  relationship.  One  of  the 
most  important  parameters  used  to  characterize  the  seismic  signature  of  an  underground 
explosion  is  the  body-wave  magnitude.  mb .  Thus  obtaining  an  unbiased  measurement  of  mb 
(or  auxiliary  Ms,  P coda,  mb(Lg),  M0,  and  RMS  Lg  values)  is  obviously  a  key  step  in  estimat¬ 
ing  the  yield.  During  the  past  decade,  the  mb  which  is  averaged  over  a  well-distributed  global 
network  and  which  incorporates  the  maximum-likelihood  technique  into  the  inversion  scheme 
has  become  widely  accepted  as  a  means  to  obtain  mb  estimates  that  avoid  bias  due  to  the 
detection  threshold  characteristics  of  individual  network  stations. 

Recently  Soviet  seismologists  have  published  descriptions  of  96  nuclear  explosions  con¬ 
ducted  from  1961  through  1972  at  the  Semipalatinsk  Test  Site,  in  Eastern  Kazakhstan.  With 
the  exception  of  releasing  news  about  their  "peaceful  nuclear  explosions”  [PNE],  the  Soviets 
have  never  before  published  such  a  body  of  information.  However,  out  of  the  72  Degelen 
events  with  announced  yields,  only  9  events  or  12.5%  were  of  “known”  yields.  The  remaining 
were  either  left  censored  (66.7%)  or  bounded  (20.8%).  Similar  heavy-censoring  pattern  can 
be  found  for  other  test  sites.  Thus  the  development  of  a  procedure  capable  of  making  full  use 
of  such  censored  information  would  seem  very  timely  and  necessary. 

In  section  I  of  this  report,  we  present  a  maximum-likelihood  regression  scheme,  “MLE- 
CY”,  which  takes  all  the  censored  yields  into  account  to  refine  the  estimated  m6:yield  relation¬ 
ship.  This  regression  routine  is  very  similar  to  the  maximum-likelihood  estimator  used  in  com¬ 
puting  the  optimal  network  mb  values  based  on  the  censored  station  amplitude  measure¬ 
ments  due  to  clipping  and  to  non-detection.  In  the  non-censored  case,  it  gives  results  identical 
to  those  derived  by  the  standard  least-squares  method.  Applications  of  this  scheme  to  the 
explosions  from  several  test  sites  of  different  geology  show  that  it  is  a  superior  procedure,  as 
compared  to  the  conventional  least-squares  approach.  The  same  algorithm  can  be  applied  to 
other  magnitude  measurements  such  as  Ms,  P coda,  mb(Lg),  M0,  RMS  Lg  and  DOB  etc. 

We  have  also  conducted  a  systematic  analysis  of  the  magnitude:yield  relationship  at  five 
major  test  sites  using  miscellaneous  unclassified  magnitudes.  (A  classified  annex  using  the 
official  mb  values  will  be  furnished  separately.) 

Several  noteworthy  results  are  summarized  here: 

[1]  Including  the  censored  yields  in  the  regression  does  improve  the  accuracy  of  the  esti¬ 
mates.  In  reality,  both  the  magnitude  and  the  yield  measurements  are  subject  to  error 
Pending  the  determination  as  to  which  of  the  two  extreme  hypotheses,  namely 
o(mb)/o( Y)  =  0  and  a(mb)/o{ Y)  =  is  closer  to  the  real  situation,  we  also  included  the 
results  based  on  Ericsson's  method  with  various  o(mb )/o(Y).  As  expected,  we  can  see 
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the  smooth  transition  ot  estimated  parameters  (i.e. ,  the  slope  and  the  intercept)  as 
o{mb )/o(Y)  varies.  Thus  the  censored  cases  with  nontrivial  o(mb )/o(Y)  values  could  be 
“interpolated".  Our  maximum-likelihood  regression  scheme  and  Ericsson's  method 
represent  two  different  directions  in  extending  the  standard  least  squares. 

[2]  For  Shagan  events,  Ringdal’s  RMS  Lg  provides  the  smallest  scatter  around  the  calibra¬ 
tion  curve,  provided  that  low-yield  events  with  mb(RMS  Lg)  <  5.5  or  yield  <  40KT  are 
excluded.  Geotech’s  GLM  method  gives  network  mb  values  better  than  almost  all  other 
magnitudes  based  on  the  teleseismic  P  waves  and  log('FJ ,  in  terms  of  both  the  yield 
estimation  and  the  mb  scaling  against  Ringdal's  RMS  Lg.  For  all  five  test  sites  we  have 
compared,  mb  measurements  reported  by  ISC  and  NEIS  are  biased  high  systematically 
at  low  yields. 

[3]  A  direct  estimation  of  the  test  site  bias  suggests  that  Nuttli’s  (1987,  1988)  Degelen  puz¬ 
zle  could  be  invalid  simply  because  of  the  relatively  poorer  quality  mb( ISC)  used.  Our 
data  indicate  that  Shagan  River  Test  Site  is  more  efficient  in  exciting  teleseismic  P 
waves  than  Degelen  Mountain,  consistent  with  our  previous  modeling  study.  Also,  the 
test  site  bias  is  yield  dependent,  in  agreement  with  other  observational  study. 

[4]  We  present  an  alternative  approach  to  derive  the  mb  adjustment  converting  cratering 
shots  to  contained  explosions  of  the  same  yield.  The  correction  derived  by  this 
approach  seems  to  match  that  by  the  multichannel  deconvolution  method  rather  well. 

[5]  Degelen  Mountain  is  the  only  test  site  that  has  a  decreasing  log(Pmax/Pa)  and 
log (Pb/Pa)  with  increasing  yields.  It  is  also  the  only  test  site  for  which  the  phase  "a” 
(i.e.,  zero-crossing  to  first  peak)  shows  the  smallest  scatter  around  the  calibration  curve, 
as  compared  to  the  phases  “b"  (i.e.,  first  peak  to  first  trough)  and  “max"  (i.e.,  max 
peak-to-trough  or  trough-to-peak  in  the  first  5  seconds).  Both  the  mountainous  topogra¬ 
phy  (which  causes  complex  pP  interference)  as  well  as  the  testing  practice  (e.g. ,  the 
relatively  shallow  and  abnormal  shot  depths)  could  be  responsible.  At  Shagan  River,  the 
phase  "b"  has  the  smallest  scatter  around  the  calibration  curve.  These  observations 
confirm  the  conjecture  (DARPA,  1981)  that  in  a  proper  environment  the  first  cycle  could 
give  better  results  than  does  "max"  phase. 

[6]  The  scale  depth  for  Konystan  explosions  is  146±1  meters,  and  the  depth  of  burial  [DOB] 
is  roughly  proportional  to  the  quartic  root  of  the  yield,  rather  than  the  cubic  root  as  fre¬ 
quently  cited  at  NTS.  This  empirical  scaling  rule  is  applicable  to  Shagan  River  region,  but 
not  Degelen  Mountain.  For  Konystan  and  Shagan  regions,  the  yields  estimated  using 
depth  scaling  have  accuracy  comparable  to  those  using  mb. 
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1.0  ABSTRACT 


Officially  announced  yields  of  underground  nuclear  explosions  are  often  truncated  or 
incomplete.  So  far  such  censored  information  has  not  been  fully  utilized  in  the  determination 
of  mb:yield  calibration  curves.  In  this  study,  we  present  a  maximum-likelihood  regression 
scheme  which  takes  all  the  censored  yields  into  account  to  refine  the  empirical  mb  .yield  rela¬ 
tionship.  Preliminary  applications  of  this  scheme  to  the  explosions  from  several  test  sites  of 
different  geology  show  that  it  is  a  superior  procedure,  as  compared  to  the  conventional  least- 
squares  approach.  A  joint  and  direct  inversion  reveals  that  the  mb  bias  between  Eastern 
Kazakhstan,  U.S.S.R.,  and  Nevada  Test  Site,  U.S.,  is  about  0.40  and  0.44  at  10KT  and 
100KT,  respectively.  The  same  algorithm  can  be  applied  to  other  magnitude  measurements 
such  as  Ms,  P coda,  mb{Lg),  M0  and  RMS  Lg  values  etc. 


1.1  INTRODUCTION 


Conventional  methods  for  estimating  underground  explosion  yields  from  seismic  record¬ 
ings  are  based  on  the  use  of  some  appropriate  “magnitudeiyield"  relationship.  One  of  the 
most  important  parameters  used  to  characterize  the  seismic  signature  of  an  underground 
explosion  is  the  body-wave  magnitude,  mb.  Thus  obtaining  an  unbiased  measurement  of  mb 
(or  similarly  Ms,  P coda,  mb(Lg),  M0,  or  RMS  Lg  values  etc.)  is  obviously  a  key  step  in 
estimating  the  yield.  There  are  already  many  publications  which  describe  different  procedures 
to  infer  better  estimates  of  mb:  e.g.,  Douglas  (1966),  von  Seggern  (1973),  Ringdal  (1976), 
von  Seggern  and  Rivers  (1978),  Christoffersson  and  Ringdal  (1981),  Blandford  and  Shumway 
(1982),  Blandford  et  al.  (1983),  Lilwall  (1986),  McLaughlin  et  al.  (1988b),  Lilwall  et  al.  (1988), 
and  most  recently,  Jih  and  Shumway  (1989).  During  the  past  decade,  the  mb  which  is 
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averaged  over  a  well-distributed  global  network  and  which  incorporates  the  maximum- 
likelihood  technique  into  the  inversion  scheme  has  become  widely  accepted  as  a  means  to 
obtain  mb  estimates  that  avoid  bias  due  to  the  detection  threshold  characteristics  of  individual 
network  stations. 

Officially  announced  yields  of  underground  nuclear  explosions  are  often  truncated  or 
incomplete.  In  general  there  are  four  types  of  announced  yields  available: 

[0]  W  is  known  as  y0  KT  (e.g. ,  the  Pahute  Mesa  event  KNICKERBOCKER  [5/26/67]  had  a 
yield  of  71  KT). 

[1]  W  is  left  censored,  i.e. ,  the  exact  value  of  W  is  known  only  to  be  less  than  a  certain 
level  t,  (e.g. ,  the  Konystan,  U.S.S.R.,  event  on  [8/26/72]  had  a  yield  less  than  20KT). 

[2]  W  is  right  censored,  i.e. ,  the  exact  value  of  W  is  known  only  to  be  larger  than  a  certain 
level  t2  (e.g.,  the  Pahute  Mesa  event  HANDLEY  [3/26/70]  had  a  yield  slightly  larger  than 
1000KT),  and 

[3]  W  is  known  only  to  lie  between  two  bounds,  ta  and  tb  (e.g. ,  the  Yucca  Flat  event  FLASK 
[5/26/70]  had  a  yield  between  20  and  200KT). 

Observations  of  types  1  through  3  are  censored.  Regression  with  right-censored  data  is 
an  important  topic  in  survival  analysis  as  well  as  in  quality  control  (Schmee  and  Hahn,  1979; 
Aitkin,  1981;  and  many  others),  while  some  biochemical  and  environmental  studies  involving 
the  monitoring  of  toxic  material  or  water  quality  have  inevitably  led  to  the  analysis  of  left- 
censored  data  (e.g.,  Gleit,  1985;  Shumway  eta!.,  1989;  and  many  others).  Both  left- 
censored  and  right-censored  station  recordings  due  to  the  ambient  noise  and  signal  clipping 
are  crucial  in  the  estimation  of  network  mb  (Ringdal,  1976;  von  Seggern  and  Rivers,  1978; 
Blandford  and  Shumway,  1982;  Jih  and  Shumway,  1989).  For  yield  determination,  likewise, 
neglecting  any  of  the  three  aforementioned  censoring  patterns  could  cause  serious  bias,  not 
to  mention  the  waste  of  useful  information.  For  instance,  recently  Soviet  seismologists 
(Bocharov  et  al. ,  1989)  have  published  descriptions  of  96  nuclear  explosions  conducted  from 
1961  through  1972  at  the  Semipalatinsk  Test  Site,  in  Eastern  Kazakhstan  (Vergino,  1989). 
With  the  exception  of  releasing  news  about  their  "peaceful  nuclear  explosions”  [PNE]  (Nor- 
dyke,  1974),  the  Soviets  have  never  before  published  such  a  body  of  information.  However, 
out  of  the  72  Degelen  events  with  announced  yields,  only  9  events  or  12.5%  were  of  type  0. 
The  remaining  were  either  left  censored  (66.7%)  or  bounded  (20.8%).  The  U.S.  announced 
yields  (Springer  and  Kinaman,  1971  and  1975)  reflect  a  very  similar  heavy-censoring  pattern. 
Although  many  authors  have  approached  the  subject  of  determining  yield  from  mb  or  other 
magnitude  measures  in  a  systematic  way  (e.g.,  Evernden,  1967;  Ericsson,  1971a,  1971b; 
Springer  and  Hannon,  1973;  von  Seggern,  1977;  Dahlman  and  Israelson,  1977;  Marshall 
et  al. ,  1979;  Nuttli,  1986a,  1986b,  1988;  Heasler  etal.,  1988;  etc.),  the  huge  amount  of  cen¬ 
sored  information  has  never  been  fully  utilized  in  the  determination  of  mb: yield  calibration 
curves. 
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In  this  study,  we  present  a  maximum-likelihood  regression  scheme  which  takes  all  the 
censored  yields  into  account  to  refine  the  estimated  nvyield  relationship.  This  regression 
routine  is  very  similar  to  the  maximum-likelihood  estimator  used  in  computing  the  optimal  net¬ 
work  mb  values  based  on  the  censored  station  amplitude  measurements  due  to  clipping  and 
to  non-detection  In  the  non-censored  case,  it  gives  identical  results  as  that  derived  by  the 
standard  least-squares  method.  The  same  algorithm  can  be  applied  to  other  magnitude 
measurements  such  as  Ms,  P coda,  mb(Lg),  M0  and  RMS  Lg  values  etc. 
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1.2  MAXIMUM-LIKELIHOOD  YIELD  ESTIMATOR 


The  problem  of  estimating  the  yield  of  an  explosion  from  the  seismic  magnitude  has 
been  handled  traditionally  using  the  linear  model 

X  =  a  +  p  log(W)  +  v  =  a  +  pY  +  v  [1] 

where  X  is  the  measured  magnitude,  mb,  a  and  p  are  intercept  and  slope  estimators,  W  is  the 
yield  in  kiloton  [KT],  and  v  is  an  error  term,  v  is  assumed  to  be  a  Gaussian  random  variable 
with  mean  zero  and  standard  deviation  a.  The  linear  or  piecewise-linear  relationship  between 
the  log(yield)  and  the  log(amplitude)  is  based  on  both  observational  study  and  theoretical 
prediction  (e.g. ,  Mueller  and  Murphy,  1971;  von  Seggern  and  Blandford,  1972;  Murphy, 
1977). 

One  may  then  collect  a  number  of  “calibration  events”,  estimating  a  and  p  by  least 
squares  using  a  number  of  known  yields  and  measured  magnitudes.  This  classical  calibration 
approach  leads  to  predicting  a  future  log-yield  Y  at  mb  =  ft  by  inverting  equation  [1],  i.e. , 

[2] 

The  geometrical  interpretation  of  “ regressing  X  on  Y"  is  that  the  (&,  0)  thus  estimated  would 
be  the  optimal  solution  that  minimizes  the  sum  of  the  squared  magnitude  residuals. 
£  (  X  -  ft  -  0Y  )2  (and  hence  the  name  of  "m -regression”).  Implicitly,  an  assumption  is  been 
made  that  the  independent  variable  Y  has  nearly  perfect  accuracy  and  precision  as  compared 
to  X. 


Alternately,  one  can  estimate  k  and  X  in  the  inverse  regression  model 

Yslog(W)  =  k  +  XX  +  v'  [3] 

and  then  predict  a  future  log  yield  directly  as 

Y  =  ft  +  X  *  [4] 

Likewise,  this  so-called  "  Y-regression"  approach  implicitly  assumes  that  X  has  perfect  accu¬ 
racy  and  precision.  The  optimal  estimates  (if,  X)  are  the  ones  that  would  minimize  the  sum  of 
the  squared  log  yield  residuals,  £  (  Y  -  *  -  Xx  )2.  Thus  either  the  yield  or  the  magnitude 
must  be  regarded  as  error-free  independent  variable  in  these  two  models.  In  reality,  both  the 
mb  and  the  yield  measurements  are  subject  to  error.  At  NTS,  where  the  yields  can  be  meas¬ 
ured  using  the  radiochemical  method  with  a  precision  better  than  that  of  the  seismic  method, 
a(mb)  »  o(  log  yield  )  could  be  a  reasonable  assumption.  This  may  not  be  the  case  in  gen¬ 
eral.  however.  Note  that  [3]  can  be  rewritten  in  a  form  similar  to  [1]: 

X  =  a  +  0Y  +  v"  [3] 

with  the  transformations  a  =  -  tc /X,  p  =  1A.. 
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Now  suppose  there  are  n0,  n,,  n2,  and  n3  events  for  each  type,  respectively.  We  will 
derive  the  maximum-likelihood  formulation  for  /-regression  model  first  (Equations  (3]  and 
(3’]).  The  conditio., al  likelihood  function  of  the  censored  observations  (  y0,  t,.  t2,  ta>  tb )  given 
the  intercept  a  ,  slope  p  ,  and  the  standard  deviation  c  of  error  in  log  yield  is 

L  (  y0.  ti,  t2,  ta,  tb  I  a,  p,o)  =np(  Yj  =  yoj  I  Ot.  p.  o  )  ‘  np(Vtij 

i=t  i=1 


n2 


n3 


n  p(  Y,  >  t2j  I  a,  P,  a  )  *  n  p(  taj  <  Yj  <  tbj  I  a,  p,  o  ) 

1=1  i=  i 


[5] 


and  the  log-likelihood  function  is 


In  L  (  y0,  t„  t2,  ta>  tb  | 


a,  p,  a  )  =  -  ln(2rra2)  - 


2a 


X0j  a \2  + 

P 


n, 


n2 


I  In  0(2,,)  +  £  In  0(-z2i)  +  £  In  [<D(zbj)  -  O(z^)] 
i  =  i  i  =  i  i=  i 


[6] 


a  +  Bt  -  Xj 

where  x  the  seismic  magnitudes;  Zj  s  - - -  for  i  =  1  j,  2j,  aj,  and  bj  ; 

u  P° 

4>(u)  =  -^=ex p(^-)  and  O(u)  e  £\j>(x)dx  are  the  probability  density  function  [p.d.f.]  and  the 

cumulative  distribution  function  [c.d.f.]  of  the  standard  normal  N(0,1),  respectively;  and  y0,  t, , 
t2,  ta,  and  tb  are  the  collection  of  announced  yields.  The  specific  form  of  the  rightmost  term  in 
equation  [6]  reveals  the  necessity  of  treating  the  type  3  censored  data  as  a  separate  class 
rather  than  considering  each  of  type  3  event  as  two  separate  events  of  type  1  and  2. 


Solving 


ainL 

3o 


=  0  implies  immediately  that  the  maximum-likelihood  solution  of  o  must 


satisfy  the  following  necessary  condition: 


o(log  yield)2  = 


no 


I(  yoj  - 

i=i _ 


xpj  -  a 


P 


n  V  W  „ 
rin  +  / - z«; 

°  H^ij) 


_  »(Z2j)  2?  4>(Zb))Zb|  -  4>(za|)zaj 

2'  H  °(Zbj)  "  ^(Zaj) 


[7] 


(a,  p,  a)  can  be  solved  iteratively  with  the  Expectation  Maximization  (EM)  algorithm 
(Dempster  et  al. ,  1977)  as  follows: 

•  Initialization  Step: 

Infer  the  initial  guess  of  the  unknown  parameters,  (a,  p,  a),  from  the  standard  regression 
with  the  type  0  data  alone. 
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•  E  Step: 

Replace  the  censored  yields  by  their  conditional  expectations  based  on  the  current  esti¬ 
mate  of  the  parameters. 

•  M  Step: 

Recompute  o  with  [7]  and  update  a  and  p  by  regressing  with  the  refined  pseudo  obser¬ 
vations  computed  in  the  E  step.  Then  repeat  steps  E  and  M  until  a,  p,  and  o  converge. 

Dempster  et  al.  (1977)  proved  that  such  iterative  procedure  guarantees  the  monotonic 
increase  of  the  likelihood  function  of  the  new  estimate,  which  in  turn  guarantees  the  conver¬ 
gence  of  the  whole  procedure  since  the  log-likelihood  function  defined  in  [6]  is  bounded 
above,  say,  by  0. 


The  following  prerequisite  mathematics  are  used  in  the  E  step.  Let  X  be  an  arbitrary 
Gaussian  random  variable  with  mean  p.  and  variance  a2,  p.d.f.  g,  c.d.f.  G,  then 

•  E(  X  |  X  <  a  )  =  \l  -  o2g(a)/G(a) , 

•  E(X|X>a)  =  n  +  o2g(a)/G(-a) . 

.  E(X|a<X<b)  =  p-  o2[g(b)-g(a)]/[G(b)-G(a)] . 


The  calculation  of  g(x)  and  G(x)  can  be  accomplished  easily  by  the  following  transforma¬ 
tions:  g(x)  =  <t>(-^-)/a  ,  G(x)  =  ,  as  was  done  in  Equations  [6]-[7]. 

o  a 

If  we  regress  the  magnitudes  on  the  log  yields,  Equation  [7]  becomes 


where  zt  = 


o(mb)2  = 
a  +  ptj  -  Xj 


"° 

I(  a  +  Pyoj  -  Xoj ) 
H _ 


4>(Zi|)  _  g  <t>(z2|)  £  <t>(zbj)zbj  ~  <t>(Zaj)zaj 


n0  +  -  I 


'Z2j  +  E~ 


H^21j)  H^-Zss)  1  H  ^(Zbj)  -  «KZaj) 

for  i  =  1  j,  2j,  aj,  and  bj . 


[8] 


Essentially  the  same  procedure  can  be  used  for  both  the  m-  and  /-regression  models. 
The  major  difference  in  the  M  step  is  whether  we  regress  Y  on  X  (and  then  transform  k  and  X 
to  a  and  P)  or  regress  X  on  Y  to  estimate  a  and  p  directly.  The  other  minor  difference  is  in 
the  calculation  of  o  and  z.  The  o  in  [7]  represents  the  standard  deviation  of  the  residual  log 
yield,  while  the  o  in  [8]  is  actually  that  for  the  residual  magnitude.  For  the  m  -regression 
model,  a(mb)  in  [8]  is  frequently  used  as  a  measure  of  goodness  of  fit.  If  the  /-regression 
model  is  used,  o(mb)  can  be  computed  as  o(log  yield)‘p.  The  2o  uncertainty  factor  in  yield  is 
defined  as  10**(2o/p)  and  10“(2o)  for  the  m-  and  /-regression  models,  respectively. 
Recently  the  m -regression  model  has  been  given  greater  attention  in  the  nuclear  monitoring 
study,  and  hence  examples  and  discussions  in  the  subsequent  sections  will  be  limited  to  the 
m-regression  model  for  brevity.  If  n,  =  n2  =  n3  =  0,  then  both  algorithms  presented  here 
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reduce  to  the  standard  least-squares  method,  and  o  in  [7]  and  [8]  becomes  the  simple  RMS 
residuals  in  the  usual  sense. 


1.3  ILLUSTRATIVE  EXAMPLES 


During  the  past  several  years,  WWSSN  (World-Wide  Standard  Seismograph  Network) 
mb  database  measured  at  Teledyne  Geotech  (TG)  has  been  gradually  expanded  to  124 
events,  totaling  366  usable  “a”,  “b",  and  “max”  event  phases  (Blandford  and  Shumway, 
1982;  Blandford  eta!.,  1983;  McLaughlin  etal.,  1988b;  Jih  and  Shumway,  1989;  Jih  etal., 
1990a;  1990b).  We  have  applied  the  maximum-likelihood  network  mb  estimator,  GLM  [Gen¬ 
eral  Linear  Model]  (Blandford  and  Shumway,  1982),  to  the  complete  data  set  consisting  of 
15,288  teleseismic  magnitude  measurements  in  the  distance  range  from  20  degrees  to  95 
degrees  at  127  stations  to  determine  our  best  mb  values  to  date,  which  we  denote  as 
mb( TG).  The  mb(Pmax  ,TG)  and  mb(Pb, TG)  of  the  events  from  the  same  test  site  are  then  fed 
to  the  maximum-likelihood  mb:yield  regression  scheme  we  just  proposed  to  derive  the  optimal 
calibration  curve.  The  resulting  calibration  curves  are  summarized  as  follow: 

(#1)  mb(Pm aX,TG)  =  3.747[±0.075]  +  0.857[±0.034)  log(W)  for  NTS  shots  in  high-coupling 
media;  a  =  0.091;  95%  confidence  factor  =  1.630;  i.e.,  we  are  95%  confident  that  the 
actual  yield  lies  in  the  range  from  Yest/1 .630  to  Yest*1.630.  In  this  regression 
(n0,  n,,  n2,  n3)  =  (9,2,1, 2).1  • 

(#2)  mb(Pb, TG)  =  3.484[±0.089]  +  0.866[±0.040]  log(W)  for  NTS  shots  in  high-coupling  media; 
c  =  0.108;  95%  confidence  factor  =  1.775.  (n0,  n,.  n2,  n3)  =  (9, 2, 1,2). 

(#3)  mb(Pmax, TG)  =  3.659[±0.022]  +  1 .008[±0.018]  log(W)  for  Sahara  and  NTS  shots  in  gran¬ 
ite;  a  =  0.032;  95%  confidence  factor  =  1.157;  (n0,  nt,  n2,  n3)  =  (4, 6, 1,0)  (cf.  Table  1). 

(#4)  mb(Pb, TG)  =  3.348[±0.028]  +  1.040[±0.022]  log(W)  for  Sahara  and  NTS  shots  in  granite; 
a  =  0.037;  95%  confidence  factor  =  1.178;  (n0,  n1f  n2,  n3)  =  (4, 6, 1,0)  (cf.  Table  1). 

(#5)  rrib(Pma),, TG)  =  4.1 10[±0.062]  +  0.892[±0.039]  log(W)  at  Eastern  Kazakhstan;  o  =  0.093; 
95%  confidence  factor  =  1.617.  (n0,  n1f  n2,  n3)  =  (13,3,0,5). 

(#6)  mb(Pb, TG)  =  3.837[±0.059]  +  0.924[±0.037]  log(W)  at  Eastern  Kazakhstan;  a  =  0.091; 
95%  confidence  factor  =  1.571.  (n0,  n,,  n2,  n3)  =  (13,3,0,5). 

Although  the  formulae  (#1)  through  (#6)  are  preliminary,  there  are  a  few  observations 
worth  noting.  First,  the  slope  in  (#1)  matches  Murphy's  (1977)  theoretical  prediction,  0.85, 

’9  events  of  type  0:  BILBY,  SHOAL,  HANDCAR,  REX,  CHARTREUSE,  PILEDRIVER,  SCOTCH,  BOXCAR,  and  BENHAM. 
2  events  of  type  1:  ALMENORO  and  MAST;  1  event  of  type  2:  HANDLEY;  2  events  of  type  3:  CORDUROV  and  NASH.  See  Ap¬ 
pendix  (page  64)  for  the  07b  values. 
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rather  well.  Secondly,  putting  the  representative  explosions  from  various  test  sites  recorded  at 
a  global  network  (such  as  WWSSN)  into  a  single  GLM  inversion  not  only  yields  a  consistent 
set  of  station  corrections  for  global  use,  but  it  also  provides  a  direct  estimate  of  the  mb  bias 
between  any  two  test  sites  of  interest.  For  instance,  the  mb  bias  between  Eastern  Kazakhstan 
and  NTS  can  be  estimated  easily  from  (#1)  and  (#5)  as  0.40  and  0.436  magnitude  unit  at 
1 0KT  and  100KT,  respectively.  This  value  is  very  close  to  that  based  on  some  indirect 
methods  using  Pn  velocity  or  surface  waves  (Evernden  and  Marsh,  1987),  and  slightly  larger 
than  that  in  Der  et  al.  (1985)  and  Stewart  (1988).  It  includes  the  combined  effects  of  the  net 
bias  due  to  the  clustering  of  stations  on  the  focal  sphere  (McLaughlin,  1988)  as  well  as  the 
difference  of  Q,  coupling,  and  pP  interferences  between  two  test  sites. 

In  deriving  (#3)  and  (#4),  we  have  supplemented  the  French  explosions  in  Hoggar  Mas¬ 
sif,  south  Algeria,  with  U.S.  shots  PILEDRIVER  and  SHOAL  detonated  at  Climax  Stock, 
Nevada.  The  French  Test  Site  is  in  the  volcanic  terrain,  apparently  in  an  incipient  rift  zone 
(Duclaux  and  Michaud,  1970;  Schock  et  al.,  1972;  Faure,  1972).  The  t*  of  Hoggar  Massif  as 
estimated  by  Der  et  al.  (1985)  is  0.35  sec,  which  shows  no  significant  difference  in  the 
attenuation  from  that  of  NTS  (McLaughlin  et  al. ,  1 988a).  There  exists  fair  agreement  between 
U.S.  and  French  granite  shots  in  the  yield-scaled  peak  values  of  acceleration,  velocity,  and 
displacement  (Heuze,  1983).  On  the  other  hand,  although  the  Semipalatinsk  Test  Site  of 
U.S.S.R.  has  hard-rock  geology  as  well,  it  is  inappropriate  to  include  the  Soviet  events  in  the 
same  regression  with  French  explosions  unless  care  is  taken  in  advance  to  correct  for  the  test 
site  bias.  Table  1  lists  the  regression  results  using  the  least-squares  (LS)  and  the  maximum- 
likelihood  estimator  (MLE)  along  with  the  announced  yields  of  U.S.  and  French  tests  in  granite 
taken  from  Bolt  (1976)  and  Stimpson  (1988).  The  yield  estimates  in  column  “MLE”  of  Table  1 
are  predicted  by  formulae  (#3)  and  (#4),  respectively.  Although  the  network  mb  values  we 
use  are  not  corrected  for  the  pP  interference  as  suggested  in  Marshall  et  al.  (1979),  they  fit 
the  theoretical  scaling  rather  well.  The  slope  of  the  mb  yield  curve  for  this  region  is  nearly  1 
for  these  low-yield  tests,  consistent  with  an  earlier  study  by  Blandford  and  Shumway  (1982) 
using  fewer  events.  Because  of  the  nearly  ideal  fit,  the  MLE  changes  the  yields  only  slightly 
as  estimated  by  the  standard  least  squares  method  in  this  particular  case. 
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Table  1.  Estimated  Yield  of  French  and  U.S.  Explosions  in  Granite 

Event  Announced  mb(Pmax  ,TG)  LS  MLE  mb{Pb, TG)  LS  MLE 

—  —  —  —  — 


BERYL 

>20.0 

4.986 

20.6 

20.8 

4.779 

23.8 

23.8 

CORUNDON 

<20.0 

4.214 

3.5 

3.6 

3.900 

mm 

mm 

EMERAUDE 

<20.0 

4.569 

7.9 

8.0 

4.263 

mm 

GRENAT 

<20.0 

4.766 

12.4 

12.6 

4.497 

12.7 

12.7 

OPALE 

<20.0 

3.894 

1.7 

1.7 

3.853 

3.1 

3.1 

RUBIS 

52.0 

5.432 

57.2 

57.5 

5.170 

56.4 

56.5 

SAPHIR 

120.0 

5.720 

110.7 

111.1 

5.468 

109.2 

109.2 

TOURMALINE 

<20.0 

4.646 

9.5 

4.429 

10.9 

11.0 

TURQUOISE 

<20.0 

4.223 

3.6 

3.6 

3.942 

mm 

mm 

SHOAL  12.2  4.739  11.7  11.8  4,455  11.6  11.6 

PILEDRIVER  56.0  5.436  57.7  58.0  5.195  59.7  59.7 


We  have  also  derived  the  maximum-likelihood  calibration  curves  using  Nuttli’s  (1986a) 
mb(Lg)  as  well  as  Marshall’s  (1988)  mb  values  for  NTS: 

(#7)  Marshall’s  mb  =  3.892[±0.105]  +  0.833[±0.049]  log(W)  for  high-coupling  material  at  NTS; 
o  =  0.186;  95%  confidence  factor  =  2.799.  Note  that  the  mean  slope,  0.833,  is  very 
close  to  that  in  (#1).  (n0,  n1t  n2,  n3)  =  (19,13,1,27). 

(#8)  Nuttli's  mb(Lg)  =  4.402[±0.038j  +  0.730[±0.018]  log(W)  for  high-coupling  material  at  NTS; 
a  =  0.086;  95%  confidence  factor  =  1.717.  (n0,  nlP  n2,  n3)  =  (22,14,1,30). 

(#9)  Nuttli’s  mb(Lg)  =  4.020(±0.038]  +  0.841  [±0.032]  log(W)  for  low-coupling  material  at  NTS;  o 
=  0.170;  95%  confidence  factor  =  2.536.  (n0,  nlt  n2,  n3)  =  (14,41,0,24). 

To  illustrate  the  robustness  of  the  present  approach,  we  have  tabulated  below  (Table  2) 
the  best  yield  estimate  of  several  often  analyzed  nuclear  tests  in  hard  rock  computed  using 
our  formulae  based  on  Marshall’s,  Nuttli’s,  and  our  magnitudes. 
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Table  2.  Comparison  of  Yield  Estimate  of  3  Granite  Shots 

Events 

RUBIS 

SAPHIR 

Kazakhstan  01/15/65 

Announced  Yield 

52KT 

120KT 

100-1 50 KT 

Nordyke"1 

125KT 

Dahlman  and  Israelson*2 

110KT 

Marshall  et  al.  *3 

45.4KT 

[m0=5.97] 

91.6KT 

[m0=6.29] 

68.9KT 

[mQ=6.16] 

Nuttli*4  mb{Lg) ,  Quadratic  Fit 
Nuttli*5  mb(Lg) ,  Linear  Fit 

68KT 

70KT 

[mb(Lg)=  5.72] 

110KT 

117KT 

[mb(Lff)=  5.89] 

103KT 

109KT 

[mb(Lg)=  5.87] 

Stimpson’6 

68KT 

[mfc=5.49] 

127KT 
[mb=  5.70] 

This  Study,  mfa(Pmax  ,TG) 

57.5KT 

lmb(P max)=5-432] 

111.1  KT 

[^(^max)=5.720] 

96.7KT 

[^(Pmax)=5.882] 

This  Study,  mb{Pb, TG) 

56.5KT 

[mb(Pb)=5A7Q] 

109.2KT 

[mb(Pb)=5A68] 

112.9KT 

[m<,(Pt,)=5.735] 

*1)  Nordyke  (1974):  based  on  the  crater  size. 

*2)  Dahlman  and  Israelson  (1977):  slope  =  0.74. 

*3)  Marshall  et  al.  (1979):  mQ  =  4.23[±0.15]  +  1.05[±0.06]  log(W)  for  salt  and  granite 
*4)  Nuttli  (1986a,  b):  mb{Lg)  =  3.943  +  1.124  log(W)  -  0.0829  (log(W))2  for  5.2  <  mb(Lg)  <  6.7 
*5)  Nuttli  (1986a,  b):  mb(Lg)  =  4.307  +  0.765  log(W)  for  5.2  <  mb(Lg )  <  6.7 
*6)  Stimpson  (1988):  mb  =  4.08  +  0.77  log(W)  for  hard  rock 


Patton  (1988),  Ringdal  and  Marshall  (1989),  Hansen  et  al.  (1989),  and  Ringdal  and  Han¬ 
sen  (1989)  confirmed  that  the  Lg  phase  is  very  promising  for  use  in  yield  estimation,  as  origi¬ 
nally  proposed  by  Nuttli  (1986a,  1986b).  Table  2  indicates  that  the  yields  estimated  by  our 
MLE  regression  scheme  using  our  mb  measurements  have  equally  good  or  better  accuracy 
as  does  mb(Lg) .  The  improvements  over  other  conventional  regression  schemes  can  be  attri¬ 
buted  to  two  factors: 

[1]  the  maximum-likelihood  magnitude:yield  regression  method  presented  in  this  study  is 
superior  to  the  conventional  least-squares  magnitude:yield  regression,  regardless  of  what 
magnitude  is  used;  and 

[2]  Geotech’s  GLM  method  results  in  a  smaller  bias  in  the  network  mb  estimates. 


To  explore  the  validity  of  the  first  claim,  we  analyzed  7  events  with  announced  yields 
from  the  Shagan  River  Test  Site  ( i.e .,  Balapan  region)  as  listed  in  Table  3  using  Marshall’s 
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(1987)  mb  measurements  and  those  of  Sykes  and  Ruggi  (1989). 

When  the  standard  least-squares  (LS)  is  applied  to  Marshall’s  mb  values  of  the  four 
events  with  known  yields,  the  predicted  yield  of  event  [01/15/65]  is  87.9KT.  Once  the  remain¬ 
ing  events  of  censored  yields  are  added  into  our  maximum-likelihood  regression,  the  estimate 
becomes  92.0KT.  If  Sykes’  mb  values  were  used  instead,  the  yield  estimate  would  change 
from  91.2KT  (LS)  to  94.9KT  (MLE).  Both  cases  show  an  obvious  improvement  relative  to  the 
announced  yields  by  incorporating  the  censored  yields  into  the  regression.  Furthermore,  such 
improvement  is  not  an  isolated  case.  Out  of  4  events  with  known  yield,  3  events  had 
significantly  improved  estimates. 


Table  3.  Explosions  at  Shagan  River  Area  with  Announced  Yield 

Date 

Lat 

Long 

Depth 

Yield 

NEIS 

Sykes 

Marshall 

[N] 

[E] 

[m] 

[KT] 

mb 

mb 

mb 

650115 

49.9350 

79.0094 

178 

100-150 

6.3 

5.905 

5.931 

680619 

49.9803 

78.9855 

316 

<20 

5.5 

5.350 

5.354 

691130 

49.9243 

78.9558 

472 

125 

6.0 

5.954 

6.048 

710630 

49.9460 

78.9805 

E9 

<20 

5.4 

5.290 

720210 

50.0243 

78.8781 

16 

5.5 

5.370 

721102 

49.9270 

78.8173 

521 

165 

6.2 

6.181 

6.224 

721210 

50.0270 

78.9956 

478 

140 

6.0 

5.989 

5.996 

[from  Bocharov  et  al.  (1989)  and  Vergino  (1989)] 

The  maximum-likelihood  calibration  curves  at  Shagan  River  region  using  mb  values  in 
Table  3  ( (n0,  n,,  n2,  n3)  =  (4, 2, 0,1) )  are  listed  as  follow: 

(#10) 

Marshall’s  mb  =  4.476[±0.090]  +  0.741[±0.052]  log(W)  with  95%  confidence  factor  1.605 
and  o  0.076. 

(#11) 

Sykes'  mb  =  4.525[±0.096]  +  0.698[±0.054]  log(W)  with  95%  confidence  factor  1 .577  and 
o  0.069. 

(#12) 

NEIS’  mb  =  4.807[±0.164]  +  0.614[±0.093]  log(W)  with  95%  confidence  factor  2.671  and 
o  0.131. 

Bocharov  et  al.  (1989)  and  Vergino  (1989)  also  listed  the  yields  of  Soviet  nuclear  explo¬ 
sions  in  Konystan  (Murzhik)  and  Degelen  regions.  Using  Marshall’s  mb  measurements,  the 
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maximum-likelinood  calibration  curves  are: 

(#13) 

Marshall's  mb  =  4.535[±0.045]  +  0.768[±0.039]  log(W)  at  Konystan;  a  =  0.069;  95% 
confidence  factor  =  1.516.  (n0,  n,.  n2,  n3)  =  (6, 7, 0,1). 

(#14) 

Marshall's  mb  =  4.370[±0.020]  +  0.869[±0.017]  log(W)  at  Degelen  with  a  0.076  and  95% 
confidence  factor  1.494.  (n0,  n,.  n2,  n3)  =  (9,46,0,15). 

It  remains  to  examine  the  second  claim  we  made  earlier,  namely  that  the  mb  values 
computed  with  Geotech’s  GLM  method  are  better  than  those  computed  with  other  methods. 
We  have  separately  regressed  Marshall’s  (1988)  and  our  mb  on  the  announced  yields  of  the 
high-coupling  shots  detonated  at  NTS,  using  the  same  maximum-likelihood  regression 
scheme.  The  results  in  Table  4  clearly  indicate  that  for  each  event  in  common,  our  predicted 
yield  is  systematically  closer  to  the  announced  yield  than  that  based  on  Marshall's  mb  values. 
About  half  of  Geotech’s  NTS  events  in  common  with  Nuttli's  have  the  announced  yields  closer 
to  the  predictions  based  on  our  formula  derived  with  Nuttli’s  mb(Lg) ,  in  accordance  with 
Nuttli's  claim  that  mb(Lg)  could  provide  yield  estimates  as  good  as  those  based  on  the  “good" 
mb. 
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Table  4.  Comparison  of  NTS  Yield  Estimates 

Date 

Event  Code 

Announced  W 

Estimated  W 

[KT] 

mb(Lg) 

mb  (Marshall) 

(^max  >  TG) 

621005 

MISSISSIPPI 

110 

87.8 

630913 

BILBY 

235 

205.8 

171.5 

631026 

SHOAL 

12.2 

12.0 

14.4 

641105 

HANDCAR 

12 

7.3 

7.5 

10.7 

660224 

REX 

16 

16.0 

8.3 

13.7 

660414 

DURYEA 

65 

53.0 

28.0 

660506 

CHARTREUSE 

70 

72.6 

44.6 

56.5 

660527 

DISCUSTHROWER 

21 

14.5 

11.9 

660602 

PILEDRIVER 

56 

93.5 

113.9 

93.7 

660630 

HALFBEAK 

300 

351.9 

412.1 

661220 

GREELEY 

825 

727.2 

644.9 

670520 

COMMODORE 

250 

175.7 

229.3 

670523 

SCOTCH 

150 

199.4 

131.5 

146 

670526 

KNICKERBOCKER 

71 

70.4 

51.5 

680426 

BOXCAR 

1200 

1096 

825 

1293 

681219 

BENHAM 

1100 

1205 

899 

1122 

691029 

CALABASH 

110 

106.1 

113.0 

700526b 

FLASK 

105 

99.6 

98.9 

701217 

CARPETBAG 

220 

240.9 

183.3 

701218 

BANEBERRY 

10 

12.8 

21.9 

710708 

MINIATA 

80 

124.2 

82.4 

730426 

STARWORT 

85 

96.5 

100.0 

_ 

#  of  Events 

22+14+1+30 

19+13+1+27 

9+2+1 +2 

6mle  (mt>) 

0.086 

0.186 

0.091 

2o  Factor 

1.717 

2.799 

W (mb(Lg) )  estimated  with  the  formula  #8. 
W(m6  .Marshall)  estimated  with  the  formula  #7. 
W(mb,TG)  estimated  with  the  formula  #1. 
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Patton  (1988)  repeated  Nuttli's  (1986a)  procedure  to  estimate  the  yields  of  69  NTS  high- 
coupling  shots  recorded  at  LLNL’s  high-quality  digital  network.  Based  on  his  regression 
result,  the  predicted  mb(Lg)  at  explosive  yields  of  10,  50,  100,  150KT  are  5.159,  5.687,  5.914, 
and  6.047,  respectively.  Nuttli’s  (1986a)  original  regression  with  22  NTS  shots  recorded  at 
WWSSN  stations  gave  5.072,  5.607,  5.837,  and  5.972,  respectively  (cf.  formula  *5  in  Table 
2).  Our  formula  (#8),  which  is  based  on  Nuttli’s  (1986a)  mb(Lg)  measurements  exclusively, 
gives  5.132,  5.642,  5.861,  and  5.990  at  10,  50,  100,  and  150KT,  respectively.  It  is  obvious 
that  our  maximum-likelihood  scheme  gives  mb(Lg)  estimates  closer  to  Patton’s  results  at  all 
levels  of  explosive  yield.  In  other  words,  including  the  censored  information  in  the  regression 
as  proposed  in  this  study  does  improve  the  determination  of  the  calibration  curve,  regardless 
of  what  type  of  magnitude  is  used. 


1.4  DISCUSSION  AND  CONCLUSIONS 


Officially  announced  yields  of  underground  nuclear  explosions  are  often  truncated  or 
incomplete.  In  this  study,  we  have  presented  a  maximum-likelihood  regression  scheme  which 
takes  all  the  censored  yields  into  account  to  refine  the  estimated  mb  ryield  relationship  with  an 
attempt  to  make  the  maximum  use  of  the  available  data.  Preliminary  applications  of  this 
scheme  to  events  from  several  test  sites  of  different  geology  show  that  it  is  indeed  a  superior 
procedure,  as  compared  to  the  conventional  least-squares  approach.  The  same  algorithm  can 
be  applied  to  other  magnitude  measurements  such  as  Ms,  mb(Lg)  or  RMS  Lg  values  etc. 

Nuttli’s  Lg  work  (1986a,  1986b)  proposed  that  careful  analysis  of  Lg  peak  amplitude  data 
from  explosions  could  produce  yield  estimates  nearly  as  accurate  as  the  best  teleseismic  esti¬ 
mates.  Based  on  the  assumption  that  his  mb{Lg)  yield  formulae  are  site  independent,  he 
obtained  a  mb  bias  estimate  (relative  to  NTS)  of  0.35  and  0.54  at  Shagan  River  and  Degelen 
Mountain,  respectively.  The  combination  of  these  two  values  would  seem  to  be  consistent 
with  our  preliminary  mb  bias  estimates  of  0.40  (10KT)  and  0.435  (100KT)  based  on  events 
from  Eastern  Kazakhstan  including  Shagan  River  and  Degelen  Mountain. 

Our  regression  with  Marshall’s  (1987)  mb  values  suggests  that  there  is  a  mb  bias  of 
0.112  and  0.150  at  Konystan  and  Degelen,  respectively,  relative  to  Shagan  River  for  100KT 
shots.  At  150KT,  the  bias  becomes  0.117  and  0.173,  respectively.  Marshall’s  mb  values  gen¬ 
erally  have  better  quality  than  the  ISC  (International  Seismological  Centre,  Newbury,  U.K.) 
bulletin  data  which  Nuttli  (1987)  used.  Thus  combining  this  mb  (Marshall)  bias  estimate  with 
Nuttli's  mb  -  mb(Lg)  offset,  0.23,  would  imply  that  there  is  a  mb(Lg)  bias  of  approximately  0.23 
-  0.15  =  0.08  and  0.23  -  0.173  =  0.057  at  100KT  and  150KT,  respectively,  between  Shagan 
River  and  Degelen  Mountain.  Linear  finite-difference  calculations  by  Jih  and  McLaughlin 
(1988)  and  Jih  et  al.  (1989)  also  suggest  that  there  should  be  observable  coupling  variations 
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affecting  Lg  amplitude.  We  are  currently  expanding  Geotech’s  mb  database  to  investigate 
such  spatial  variation  among  three  subregions  of  Eastern  Kazakhstan  (Jih  et  at. ,  1990a; 
1990b).  At  any  rate,  our  preliminary  analysis  using  Nuttli's  mb(L.g)  values  tends  to  suggest  that 
the  regionalized  calibration  curves  should  provide  a  better  result.  For  instance,  formulae  (#8) 
and  (#9)  would  give  a  better  fit  than  the  formula  (*5)  in  Table  2  for  NTS  events.  In  principle, 
this  should  be  true  not  just  for  mb(Lg)  alone.  Porting  any  empirical  magnitude:yield  calibration 
curve  from  one  site  to  another  could  be  unreliable  in  some  cases.  The  difference  between  for¬ 
mula  (#10)  for  Shagan  River  and  formula  (#14)  for  Degelen  Mountain  is  an  example. 

Recent  theoretical  studies  on  Lg  (Lilwall,  1988;  Jih  eta!.,  1989;  Frankel,  1989)  seem  to 
agree  that  in  a  medium  where  the  velocity  increases  with  depth  a  smaller  and  smaller  focal 
sphere  of  pS  will  be  trapped  as  depth  increases,  thus  decreasing  the  Lg  amplitude.  Since  the 
larger  shots  are  buried  more  deeply,  this  would  imply  that  in  general  the  slope  in  mb{Lg)  :yield 
relationship  would  be  less  than  that  in  the  mb:yield  relationship,  as  indicated  in  formulae  (#1) 
through  (#9). 

Special  purpose  magnitudes,  like  mQ  in  Marshall  et  at.  (1979)  which  include  corrections 
for  source  depth  and  source  region  attenuation  should  be,  in  principle,  superior  to  mb  for 
estimating  the  explosive  yield.  However,  the  present  study  has  shown  that  this  may  not  be  the 
case  (cf.  Table  2).  The  success  of  the  pP  and  t*  corrections  depends  on  the  accuracy  of  the 
corrections.  In  our  examples,  the  network  mb  (or,  m2  in  Marshall  etal.,  1979),  which  were 
only  corrected  for  the  instrument  gain,  geometrical  spreading  (Veith  and  Clawson,  1972)  as 
well  as  the  station  terms,  would  give  fairly  good  yield  estimates.  Finally,  the  results  in  Table  2 
seem  to  indicate  that  the  phase  "b”  ( i.e . ,  the  first  peak  to  the  first  trough)  of  the  teleseismic  P 
wave  could  give  the  yield  estimate  equally  well  as  does  the  phase  “max”.  However,  further 
investigation  is  necessary. 
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SECTION  II 

MAGNITUDE:YIELD  RELATIONSHIP  AT  VARIOUS  TEST  SITES 

Rong-Song  Jih,  Robert  A.  Wagner,  and  T.  W.  McEifresh 

Teledyne  Geotech  Alexandria  Laboratories 
314  Montgomery  Street 
Alexandria,  VA  22314-1581 

11.1  SUMMARY 

We  have  conducted  a  systematic  analysis  of  the  magnituderyield  relationship  at  several 
test  sites  using  miscellaneous  magnitudes.  The  main  tool  of  this  study  is  a  linear-regression 
scheme  “MLE-CY"  (Jih  etal.,  1990a;  1990b)  which  takes  all  censored  yields  [e.g. ,  yield  <  20 
KT  or  100  KT  <  yield  <  150  KT)  into  account  to  refine  the  determination  of  the  calibration 
curve.  The  majority  of  the  recently  published  96  Soviet  explosive  yields  (Bocharov  etal., 
1989;  Vergino,  1989)  and  the  U.S.  announced  yields  (Springer  and  Kinaman,  1971,  1975) 
were  heavily  truncated  or  rounded,  and  hence  the  maximum-likelihood  approach  would  seem 
ideal  to  make  full  use  of  the  yield  information.  The  regression  routine  we  use  is  very  similar 
to  the  maximum-likelihood  estimator  used  in  computing  the  optimal  network  mb  values  based 
on  the  censored  station  amplitude  measurements  due  to  clipping  and  to  noise  (Blandford  and 
Shumway,  1982;  Jih  and  Shumway,  1989).  In  the  non-censored  case,  it  gives  results  identi¬ 
cal  to  those  derived  by  the  standard  least  squares,  corresponding  to  the  two  extreme  cases  of 
Ericsson's  (1971)  curve-fitting  method  which  puts  different  variances  in  both  the  independent 
and  the  dependent  variables. 

In  the  following  sections,  we  will  tabulate  the  maximum-likelihood  mb: yield  calibration 
curves  which  symbolically  correspond  to  c(mb)/c( Y)  =  0  and  °°,  respectively.  Several 
noteworthy  results  are  summarized  here: 

[1]  Including  the  censored  yields  in  the  regression  does  improve  the  accuracy  of  the  esti¬ 
mates  (cf.  Tables  2C  and  2D).  In  reality,  both  the  magnitude  and  the  yield  measure¬ 
ments  are  subject  to  error.  Pending  the  determination  as  to  which  of  the  two  extreme 
hypotheses,  namely  o(m6)/o(Y)  =  0  and  o(mb)/a(Y)  =  °°,  is  closer  to  the  real  situation, 
we  also  included  the  results  based  on  Ericsson’s  method  with  various  o(mb)/o( Y).  As 
expected,  we  can  see  the  smooth  transition  of  estimated  parameters  (i.e. ,  the  slope  and 
the  intercept)  as  o(mb)/a(Y)  varies  [cf.  Tables  2A,  4A,  5C,  6C,  and  7B).  Thus  the  cen¬ 
sored  cases  with  nontrivial  o[mb  )/o(Y)  values  could  be  “interpolated".  Our  maximum- 
likelihood  regression  scheme  and  Ericsson's  method  represent  two  different  directions  in 
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extending  the  standard  least  squares.  In  the  future,  Efron's  bootstrap  (Efron,  1979;  Efron 
and  Tibshirani,  1985)  or  other  resampling  techniques  could  be  incorporated  into 
Ericsson's  curve-fitting  routine  to  estimate  the  confidence  interval. 

[2]  For  Shagan  events,  Ringdal's  RMS  Lg  provides  the  smallest  scatter  around  the  calibra¬ 
tion  curve,  provided  that  low-yield  events  with  mb(RMS  Lg)  <  5.5  or  yield  <  40KT  (e.g. 
the  explosion  on  10  Feb  72)  are  excluded.  Geotech's  GLM  method  (Blandford  and 
Shumway,  1982)  gives  network  mb  values  better  than  almost  all  other  magnitudes  based 
on  teleseismic  P  waves  and  log(¥«,) .  in  terms  of  both  the  yield  estimation  (cf.  Tables  2B 
and  5C)  and  the  mb  scaling  against  Ringdal's  RMS  Lg  (cf.  Tables  5F  and  9A).  For  all 
five  test  sites  we  have  compared,  mb  measurements  reported  by  ISC  and  NEIS  are 
biased  high  systematically  at  low  yields  (cf.  Tables  2C,  4D,  5E,  and  5D). 

[3]  A  direct  estimation  of  the  test  site  bias  (cf.  Tables  9A  and  9B)  suggests  that  Nuttli’s 
(1987,  1988)  Degelen  puzzle  could  be  invalid  simply  because  of  the  relatively  poorer 
quality  mb( ISC)  used.  Our  data  indicate  that  the  Shagan  River  Test  Site  is  more  efficient 
in  exciting  teleseismic  P  waves  than  Degelen  Mountain,  consistent  with  our  previous 
modeling  study  (Jih  and  McLaughlin,  1988).  Also,  the  test  site  bias  is  yield  dependent,  in 
agreement  with  other  observational  study. 

[4]  We  present  an  alternative  approach  to  derive  the  mb  adjustment  converting  cratering 
shots  to  contained  explosions  of  the  same  yield  (cf.  Tables  8A  and  8B).  The  correction 
derived  by  this  approach  seems  to  match  that  by  the  multichannel  deconvolution  method 
(Der  et  al. ,  1985)  rather  well. 

[5]  Degelen  Mountain  is  the  only  test  site  that  has  a  decreasing  log (PmMlPa)  and 
log (PblPa)  with  increasing  yields  (cf.  Tables  6D,  8A,  and  8B).  It  is  also  the  only  test  site 
for  which  the  phase  "a”  (i.e.,  zero-crossing  to  first  peak)  shows  the  smallest  scatter 
around  the  calibration  curve,  as  compared  to  the  phases  “b”  (i.e.,  first  peak  to  first 
trough)  and  "max”  (i.e.,  max  peak-to-trough  or  trough-to-peak  in  the  first  5  seconds). 
Both  the  mountainous  topography  (which  causes  complex  pP  interference)  as  well  as  the 
testing  practice  (e.g. ,  the  relatively  shallow  and  abnormal  shot  depths)  could  be  respon¬ 
sible.  At  Shagan  River,  the  phase  "b"  has  the  smallest  scatter  around  the  calibration 
curve  (cf.  Tables  5C  and  5D).  These  observations  confirm  the  conjecture  (DARPA, 
1981)  that  in  a  proper  environment  the  first  cycle  could  give  better  results  than  does 
"max"  phase. 

[6]  The  scale  depth  for  Konystan  explosions  is  1 46±1  meters,  and  the  depth  of  burial  [DOB] 
is  roughly  proportional  to  the  quartic  root  of  the  yield,  rather  than  the  cubic  root  as  fre¬ 
quently  cited  at  NTS  (cf.  Table  7B).  This  empirical  scaling  rule  is  applicable  to  Shagan 
River  test  site,  but  not  Degelen  Mountain.  For  Konystan  and  Shagan  regions,  the  yields 
estimated  using  depth  scaling  have  accuracy  comparable  to  those  using  mb. 
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11.2  NTS 


Table  2A.  m^Yield  Relation  of  NTS  High-Coupling  Shots 


(Earlier  Studies) 


#  of  Events1 

Magnitude 

<x™b) 

o(Y) 

Slope 

Intercept 

o(mb) 

2o  Factor2 

Method 

22+0+0+0 

Nuttli  (1986a) 

? 

0.765±0.027 

4.307+0.067 
_ _ , _ 

LS 

69+0+0+0 

Patton  (1988) 

oo 

0.755±0.022 

4.404+0.048 

0.098 

1.818 

LS 

(This  Study) 

#  of  Events 

Magnitude 

<XY) 

Slope 

Intercept 

o(mb) 

2 a  Factor 

Method 

22+14+1+30 

Nuttli,  mb{Lg) 

0 

0.761  ±0.033 

4.336+0.193 

0.116 

2.019 

MLE-CY 

22+14+1+30 

Nuttli,  mb(Lg) 

oo 

0.730±0.018 

4.402+0.038 

0.086 

1.717 

MLE-CY 

19+13+1+29 

ISC 

0 

0.787±0.067 

4.006+0.379 

0.190 

3.036 

MLE-CY 

19+13+1+29 

ISC 

oo 

0.693+0.035 

4.199+0.074 

0.136 

2.475 

MLE-CY 

19+13+1+27 

Marshall 

0 

0.982+0.062 

3.581+0.351 

0.210 

2.672 

MLE-CY 

19+13+1+27 

Marshall 

oo 

0.833+0.049 

3.892+0.105 

0.186 

2.799 

MLE-CY 

9+2+1 +2 

TG,  Pa 

0 

0.893+0.088 

3.165+0.450 

0.204 

2.863 

MLE-CY 

9+2+1 +2 

TG,  Pa 

oo 

0.835+0.065 

3.283+0.147 

0.175 

2.632 

MLE-CY 

9+2+1 +2 

TG,  Pb 

0 

0.887+0.052 

3.441+0.279 

0.124 

1.901 

MLE-CY 

9+2+1 +2 

TG.  Pb 

oo 

0.866+0.040 

3.484+0.089 

0.108 

1.775 

MLE-CY 

9+2+1 +2 

TG,  Pm ax 

0 

0.872+0.045 

3.716+0.253 

0.105 

1.744 

MLE-CY 

9+2+1 +2 

TG,  Pmax 

oo 

0.857+0.034 

3.747+0.075 

0.091 

1.630 

MLE-CY 

22+0+0+0 

Nuttli,  mb{Lg) 

0 

0.760+0.040 

4.340+0.232 

0.082 

1.646 

LS 

22+0+0+0 

Nuttli,  mb(Lg) 

0.1 

0.759+ _ 

4.342+ _ 

_ 

Ericsson 

Nuttli,  mb(Lg) 

1 

0.745+ _ 

4.370+ _ 

Ericsson 

22+0+0+0 

Nuttli,  mb(Lg) 

5 

0.737+ _ 

4.385+ _ 

Ericsson 

22+0+0+0 

Nuttli,  mb(Lg) 

100 

0.737+ _ 

4.386+ _ 

Ericsson 

22+0+0+0 

Nuttli,  mb(Lg) 

oo 

0.737+0.029 

4.386+0.060 

0.081 

1.659 

LS 

1)  #  of  "exact”  yields,  #  of  left-censored  yields,  #  of  right-censored  yields,  and  #  of  bounded 
yields. 

2)  the  multiplicative  uncertainty  factor  in  the  yield  [KT]  at  95%  confidence  level. 


-  23- 


m. -Yield  Calibration  Curves 

O 


Table  2A.  mb:Yield  Relation  of  NTS  High-Coupling  Shots  (Continued) 

(This  Studies) 

#  of  Events 

Magnitude 

oT'V 

<XY) 

Slope 

Intercept 

a(mb) 

2o  Factor 

Method 

9+0+0+0 

TG,  Pmax 

0 

0.870±0.052 

3.71 9±0. 284 

0.098 

1.684 

LS 

9+0+0+0 

TG,  Pmax 

0.1 

0.869± _ 

3.720+ _ 

. _ 

_ 

Ericsson 

9+0+0+0 

TG,  Pmax 

1 

0.860± _ 

3.738± _ 

_ 

Ericsson 

9+0+0+0 

TG,  Pmax 

5 

0.854+ _ 

3.750± _ 

_ 

Ericsson 

9+0+0+0 

TG,  Pmax 

100 

0.854± _ 

3.751  ± _ 

_ 

Ericsson 

9+0+0 +0 

TG,  Pmax 

oo 

0.854±0.044 

3.751+0.093 

0.097 

1.692 

LS 

For  purposes  of  estimating  explosion  yields,  the  media  are  divided  into  three  types: 
unsaturated  material,  e.g. ,  alluvium  and  dry  tuff;  water-saturated  rock;  and  granite.  The  U.S. 
granite  shots  PILEDRIVER  and  SHOAL  will  be  discussed  in  the  next  section  again. 

If  we  ignore  the  different  corner  frequencies  between  events  of  large  and  small  yields, 
and  put  all  high-coupling  shots  in  one  single  regression,  then  both  Geotech’s  mb(Pmax)'s  and 
Marshall’s  mb  give  a  slope  matching  Murphy’s  (1977)  theoretical  prediction,  0.85,  rather  well 
(cf.  Table  2A). 

At  NTS,  yields  estimated  from  mb  alone  have  a  random  uncertainty  factor  of  1 .45  at  the 
95%  confidence  (i.e. ,  2o)  level,  provided  the  best  “official”  mb  values  are  used  (U.S. 
Congress/OTA,  1988).  None  of  the  magnitudes  listed  in  Table  2A  reaches  such  a  precision. 
However,  it  is  also  clear  that  the  mb  based  on  our  Pmax  is  relatively  more  precise  than  other 
unclassified  mb  measurements.  The  phase  “a”  has  much  larger  variance  than  the  phase 
“max”  at  NTS,  possibly  because  of  the  small  amplitudes  measured  were  near  the  noise. 
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Table  2B.  Maximum-Likelihood  Yield  Estimates  of  NTS  Shots 

Date 

Event  Code 

Announced  W 

Estimated  W 

[KT] 

mb(Lg) 

mb  (Marshall) 

max  •  TG) 

621005 

MISSISSIPPI 

110 

87.8 

630913 

BILBY 

235 

205.8 

171.5 

631026 

SHOAL 

12.2 

12.0 

14.4 

641105 

HANDCAR 

12 

7.3 

7.5 

10.7 

660224 

REX 

16 

16.0 

8.3 

13.7 

660414 

DURYEA 

65 

53.0 

28.0 

660506 

CHARTREUSE 

70 

72.6 

44.6 

56.5 

660527 

DISCUSTHROWER 

21 

14.5 

11.9 

660602 

PILEDRIVER 

56 

93.5 

113.9 

93.7 

660630 

HALFBEAK 

300 

351.9 

412.1 

661220 

GREELEY 

825 

727.2 

644.9 

670520 

COMMODORE 

250 

175.7 

229.3 

670523 

SCOTCH 

150 

199.4 

131.5 

146 

670526 

KNICKERBOCKER 

71 

70.4 

51.5 

680426 

BOXCAR 

1200 

1096 

825 

1293 

681219 

BENHAM 

1100 

1205 

899 

1122 

691029 

CALABASH 

110 

106.1 

113.0 

700526b 

FLASK 

105 

99.6 

98.9 

701217 

CARPETBAG 

220 

240.9 

183.3 

701218 

BANEBERRY 

10 

12.8 

21.9 

710708 

MINIATA 

80 

124.2 

82.4 

730426 

STARWORT 

85 

96.5 

100.0 

_ 

#  of  Events 

22+14+1+30 

19+13+1+27 

9+2+1 +2 

dMlE  (mb) 

0.086 

0.186 

0.091 

2o  Factor 

1.717 

2.799 

1.630 

p' 

0.990 

0.942 

0.994 

*  p:  the  correlation  coefficient  between  the  magnitudes  and  the  log  yields. 
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For  each  event  in  common  with  Marshall's  in  Table  2B,  the  yield  predicted  with 
Geotech’s  mb(Pmax)  is  always  closer  to  the  announced  value  than  that  based  on  Marshall’s 
mb  values.  As  noted  in  Jih  etai  (1990a),  this  would  strongly  suggest  that  Geotech's  mb 
values  have  smaller  systematic  error,  since  the  same  regression  methodology  was  used. 

Patton  (1988)  utilized  Nuttli’s  procedure  to  estimate  the  yields  for  69  high-coupling  shots 
at  NTS.  The  NTS  explosions  Patton  used  were  clustered  around  mb(Lg)  =  5.8.  Beyond  that 
level,  the  difference  in  yield  estimates  between  Nuttli’s  and  Patton's  predictions  are  by  no 
means  negligible.  For  mb(Lg)  =  6.0,  they  predict  the  yield  to  be  163KT  (N)  and  130KT  (P), 
respectively.  At  mb(Lg)  =  6.5,  the  predictions  are  736KT  (N)  and  597KT  (P),  respectively. 

Since  Patton  (1988)  did  not  release  the  individual  yields  or  mb(Lg)  in  his  paper,  we  need 
an  alternative  approach  to  make  the  comparison.  The  data  recorded  at  LLNL’s  regional  digi¬ 
tal  network  have  quality  better  than  those  WWSSN  film  chips  which  Nuttli  (1986a)  read.  It 
would  seem  reasonable  to  assume  that  the  mb  predicted  by  Patton’s  regression  is  more  accu¬ 
rate  than  Nuttli’s.  Table  2c  below  indicates  that  regressing  Nuttli’s  mb(Lg)  measurements 
against  the  censored  yields  with  our  maximum-likelihood  scheme  gives  mb  estimates  very 
close  to  Patton’s  results  at  all  levels  of  explosive  yield.  In  other  words,  including  the  censored 
information  in  the  regression  as  proposed  in  Jih  etai  (1990a,  1990b)  does  improve  the 
determination  of  the  calibration  curves,  regardless  of  what  magnitude  is  used. 


-  26- 


m -Yield  Calibration  Curves 

D 

Table  2C.  Expected  magnitudes  of  NTS  High-Coupling  Explosions 
(Regressing  the  magnitudes  on  the  yields) 

(Earlier  Studies) 


mb  :Y  Curve 

If  of  Events 

10KT 

50KT 

100KT 

150KT 

Nuttli1 

22+0+0+0 

5.072 

5.607 

5.837 

5.972 

Patton2 

694-0+0+0 

5.159 

5.687 

5.914 

6.047 

(This  Study) 

mb  :Y  Curve 

#  of  Events 

10KT 

50KT 

100KT 

150KT 

Nuttli3 

22+0+0+0 

5.123 

5.638 

5.860 

5.989 

Nuttli4 

22+14+1+30 

5.132 

5.642 

5.861 

5.990 

ISC 

19+13+1+26 

4.892 

5.376 

5.585 

5.707 

Marshall 

19+13+1+27 

4.725 

5.307 

5.558 

5.704 

TG.  Pa 

9+2+1 +2 

4.118 

4.701 

4.953 

5.100 

TG.  Pb 

9+2+1 +2 

4.350 

4.955 

5.216 

5.368 

TG,  Pm ax 

9+2+1 +2 

4.604 

5.202 

5.460 

5.610 

1)  Nuttli  (1986a):  mb(Lg)  =  4.307[i0.067]  +  0.765[±0.027]log(W)  for  5.2  <  mb(Lg)  <  6.7. 

2)  Patton  (1988):  mb(Lg)  =  4.404(±0.048]  +  0.755[±0.022]log(W)  for  4.22  <  m^Lg)  <  6.7. 

3)  Nuttli's  mb(Lg)  values  regressed  with  the  least  square,  o(mb)/a{Y)  =  °°  (cf.  Table  2A). 

4)  Nuttli's  mb(Lg)  values  regressed  with  MLE-CY,  o(mb )/a(Y)  =  oo  (cf.  Table  2A). 


Table  2C  raises  a  question  as  how  to  evaluate  different  calibration  curves.  Apparently 
the  trade  off  between  a  and  (5  should  be  taken  into  account.  Judging  on  the  slope,  p,  alone 
could  be  very  misleading.  For  instance,  in  comparison  with  the  2  slopes  which  we  obtained 
with  Nuttli's  mb(Lg)  measurements,  his  original  slope  is  closer  to  that  of  Patton’s  (cf.  Table 
2A),  and  yet  our  formulae  actually  predict  the  yields  as  well  as  the  magnitudes  closer  to  those 
of  Patton's  (Tables  2C  and  2D). 
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1)  Nuttli  (1986a):  mb(Lg)  =  4.307[±0.067]  +  0.765[±0.027]log(W)  for  5.2  <  mb(Lg)  <  6.7. 

2)  Patton  (1988):  mb{Lg)  =  4.404[±0.048]  +  0.755[±0.022]log(W)  for  4.22  <  mb(Lg)  <  6.7 

3)  Nuttli’s  mb{Lg)  values  regressed  with  the  least  square  (cf.  Table  2A). 

4)  Nuttli’s  mb(Lg)  values  regressed  with  MLE-CY  (cf.  Table  2A). 


Due  to  the  different  yield  relationships  for  teleseismic  P  and  Lg  at  NTS,  the  yield  esti¬ 
mates  at  the  same  “magnitude"  level  are  very  different.  We  will  compare  the  mb(P)  -  mb(Lg) 
offset  of  various  test  sites  in  a  later  section  (cf.  Table  9A). 

In  comparing  with  Nuttli's  regression  results,  we  noticed  that  his  original  formula  (Equa¬ 
tion  1  in  Table  2C)  seems  not  reproducible.  His  data  set  (cf.  Nuttli,  1986a,  page  2144) 
included  the  Pahute  Mesa  event  HANDLEY  which  had  a  bounded  yield  of  >1000KT.  How¬ 
ever,  Nuttli  seemed  to  have  treated  the  yield  as  exactly  1000KT  in  his  calculations  (cf.  Fig¬ 
ures  7  and  9  of  Nuttli,  1986a).  Different  symbols  for  the  2  granite  events  PILEDRIVER  and 
SHOAL  were  used  in  his  figures  (cf.  Nuttli,  1986a,  pages  2145  and  2147).  Also,  Nuttli 
imposed  a  mb(Lg)  range  of  applicability  (from  5.2  to  6.7)  on  his  calibration  curve. 

We  have  tested  eight  possible  combinations  with  Nuttli’s  mb(Lg)  measurements: 

•  including  NTS  granite  events  or  not, 

•  limiting  mb(Lg)  to  (5.2, 6.7]  or  not, 

•  assuming  HANDLEY  was  1000KT  or  deleting  HANDLEY  from  the  regression. 
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None  of  the  eight  extra  experiments  could  give  an  "exactly  identical”  formula  to  that 
given  by  Nuttli  (1986a),  even  if  the  computer’s  "machine  e”  is  accounted  for.  It  seems  very 
likely  that  Nuttli  was  using  the  “Y'-regression"  models,  i.e.,  o(mb)/a( Y)  =  0,  with  some 
unspecified  constraint  on  the  data  set.  However,  for  all  cases  we  have  tested,  the  comparis¬ 
ons  of  MLE  results  (using  Nuttli’s  data)  against  Patton’s  result  confirmed  consistently  that 
including  the  censored  data  would  improve  the  regression. 


11.3  U.S.  AND  FRENCH  SAHARA  SHOTS  IN  GRANITE 


Table  3A.  mb:Y\e\d  Relation  of  French  Sahara  and  NTS  Events  in  Granite* 


It  of  Events 

mb 

o(m0) 

o00 

Slope 

Intercept 

o{mb) 

2 o  Factor 

Method 

4+0+0+0 

TG,  Pa 

0 

0.875±0.056 

3.365±0.270 

0.035 

1.203 

LS 

4+0+0+0 

TG,  Pa 

oo 

0.869±0.049 

3.374±0.083 

0.035 

1.203 

LS 

4+0+0+0 

1 

0 

1.048±0.064 

3.334±0.325 

0.048 

1.234 

LS 

4+0+0+0 

msm 

oo 

1.040±0.066 

3.348±0.1 13 

0.048 

1.235 

LS 

4+0+0+0 

TG,  Fmax 

0 

1.011  ±0.058 

3.657±0.310 

0.042 

1.211 

LS 

4+0+0+0 

TG,  F’max 

oo 

1.004±0.058 

3.668±0.099 

0.042 

1.212 

LS 

4+4+1 +0 

TG,  Pa 

0 

0.928±0.044 

3.258±0.195 

0.061 

1.353 

MLE-CY 

4+4+1 +0 

TG,  Pa 

oo 

0.905±0.036 

3.296i0.048 

0.056 

1.328 

MLE-CY 

4+6+1 +0 

TG,  Pb 

0 

1.049±0.020 

3.334±0.092 

0.035 

1.169 

MLE-CY 

4+6+1 +0 

TG,  Pb 

oo 

1.040±0.022 

3.348±0.028 

0.037 

1.178 

MLE-CY 

4+6+1 +0 

TG,  Pmax 

0 

1.014±0.018 

3.658±0.084 

0.032 

1.154 

MLE-CY 

4+6+1 +0 

TG,  Pmax 

oo 

1.008±0.018 

3.659±0.022 

0.032 

1.157 

MLE-CY 

*)  2  NTS  events  in  granite  and  9  French  Sahara  explosions;  no  Pa  for  EMERAUDE  and 
TURQUOISE. 
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Table  3B.  Yield  Estimates  of  French  &  U.S.  Shots  in  Granite 


Event 


Official  W 


[KT] 


^ b  (^max ) 


BERYL 

>20.0 

4.986 

CORUNDON 

<20.0 

4.214 

EMERAUDE 

<20.0 

4.569 

GRENAT 

<20.0 

4.766 

OPALE 

<20.0 

3.894 

RUBIS 

52.0 

5.432 

SAPHIR 

120.0 

5.720 

TOURMALINE 

<20.0 

4.646 

TURQUOISE 

<20.0 

4.223 

SHOAL 

12.2 

4.739 

PILEDRIVER 

56.0 

5.436 

LS 


[KT] 


20.6 


3.5 


.9 


12.4 


1.7 


57.2 


110.7 


mb(Pb) 


3.9 


4.263 


4.497 


3.853 


5.170 


5.468 


4.429 


3.942 


00 

3. 

m 

3. 

12.7 


3.1 


56.4 


109.2 


10.9 


3. 


.6 


12.7 


3.1 


56.5 


109.2 


11.0 


3. 


155 

11.6 

11.6 

95 

59.7 

59.7 

4+0+0+0 

4+6+1 +0 

0.042 

0.032 

1.212 

1.157 

0.997 

0.999 

4+0+0 +0 

4+6+1 +0 

0.048 

0.037 

1.235 

1.178 

0.996 

0.999 

#  of  Events 


dMLE  (mb) 


2a  Factor 


Table  3C.  Expected  mb  of  U.S.  and  French  Shots  in  Granite 


mb  :Y  Curve 

#  of  Events 

10KT 

50  KT 

100KT 

150KT 

TG,  Pa 

4+4+1 +0 

4.201 

5.106 

5.265 

TG  ,Pb 

4+6+1 +0 

4.388 

5.428 

5.611 
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11.4  EASTERN  KAZAKHSTAN  AREA 


Table  4A.  m6:Yield  Calibration  Curve  at  Eastern  Kazakhstan’ 


#  of  Events 

mb 

o<Y! 

Slope 

Intercept 

o(mb) 

2o  Factor 

Method 

19+55+0+17 

ISC 

0 

0.715+0.029 

4.532+0.156 

0.105 

1.972 

MLE-CY 

19+55+0+17 

ISC 

CO 

0.687±0.014 

4.570+0.018 

0.076 

1.660 

MLE-CY 

4 

19+55+0+17 

NEIS 

0 

0.745±0.039 

4.607+0.213 

0.157 

2.639 

MLE-CY 

19+55+0+17 

NEIS 

oo 

0.655±0.019 

4.725+0.023 

0.113 

2.222 

19+55+0+17 

Sykes 

0 

0.717±0.024 

4.535+0.129 

0.088 

1.755 

ppisi^j 

19+55+0+17 

E9 

oo 

0.696+0.012 

4.563+0.015 

0.063 

1.520 

MLE-CY 

19+0+0+0 

Marshall 

0 

0.823+0.050 

4.419+0.279 

0.098 

1.728 

LS 

19+0+0+0 

Marshall 

0.1 

0.822+ _ 

4.420+ _ 

Ericsson 

19+0+0+0 

Marshall 

1 

0.802+ _ 

4.448+ _ 

_ 

Ericsson 

19+0+0+0 

Marshall 

5 

0.791+ _ 

4.466+ _ 

_ 

Ericsson 

19+0+0+0 

Marshall 

oo 

0.789+0.039 

4.466+0.060 

0.096 

1.748 

LS 

19+55+0+17 

Marshall 

0 

0.798+0.025 

4.462+0.133 

0.109 

1.872 

MLE-CY 

19+55+0+17 

Marshall 

oo 

0.759+0.015 

4.516+0.018 

0.087 

1.696 

MLE-CY 

12+3+0+5 

TG,  Pa 

0 

0.951+0.042 

3.497+0.208 

0.094 

1.577 

MLE-CY 

12+3+0+5 

TG,  Pa 

oo 

0.926+0.037 

3.537+0.059 

0.088 

1.552 

MLE-CY 

13+3+0+5 

0 

0.951+0.042 

3.795+0.220 

0.096 

1.594 

MLE-CY 

13+3+0+5 

TG,  Pb 

oo 

0.924+0.037 

3.837+0.059 

0.091 

1.571 

MLE-CY 

13+3+0+5 

TG,  Pmax 

0 

0.921+0.047 

4.064+0.257 

0.102 

1.666 

MLE-CY 

_ i 

13+3+0+5 

TG,  Pmax 

oo 

0.892+0.039 

4.110+0.062 

0.093 

*)  including  Shagan  River  (Baiapan),  Konystan  (Murzhik),  and  Degelen  Mountain. 


In  Table  4A,  we  regressed  all  Eastern  Kazakh  explosions  with  announced  yields 
(Bocharov  et  al. ,  1989)  against  various  mb  values  of  Marshall  (1987),  ISC,  NEIS,  and  ours. 
Detailed  descriptions  of  the  explosions  are  listed  in  later  sections  according  to  the  subregion 
they  belong  to. 
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Table  4B.  Least-Squares  Yield  Estimates  of  E.  Kazakh  Shots 


Event,  Region 

Official  W 

ISC 

NEIS 

Sykes 

Marshall 

TG,  Prnax 

[KT] 

[KT] 

[KT] 

[KT] 

[KT] 

[KT] 

651121,  D 

29.0 

32.0 

48.0 

31.0 

27.7 

25.6 

660213,  D 

125.0 

159.6 

188.5 

155.8 

159.0 

660320,  D 

100.0 

115.7 

188.5 

112.8 

98.6 

88.2 

660507,  D 

4.0 

2.4 

1.6 

2.3 

2.2 

2.8 

670922,  M 

10.0 

8.8 

8.7 

8.5 

7.6 

680929,  D 

60.0 

60.8 

48.0 

59.1 

58.5 

50.1 

690723,  D 

16.0 

16.8 

17.2 

16.2 

20.6 

15.6 

691130,  S 

125.0 

115.7 

95.2 

97.2 

100.9 

121.2 

691228,  M 

40.0 

44.1 

34.1 

42.8 

47.7 

710425,  D 

90.0 

83.9 

67.6 

92.9 

109.5 

69.5 

710606,  M 

16.0 

23.2 

17.2 

21.0 

22.0 

711009,  M 

12.0 

12.2 

12.2 

12.5 

14.0 

711021,  M 

23.0 

23.2 

24.3 

23.1 

25.8 

__ 

720210,  S 

16.0 

16.8 

17.2 

14.7 

14.0 

22.1 

720328,  D 

6.0 

6.4 

6.2 

7.0 

8.0 

9.2 

72081 G,  D 

8.0 

4.6 

6.2 

6.8 

6.4 

7.6 

720902,  M 

2.0 

3.4 

4.4 

3.0 

2.6 

721102,  S 

165.0 

159.6 

188.5 

202.4 

168.6 

207.6 

721210,  S 

140.0 

115.7 

95.2 

108.8 

86.7 

133.4 

#  of  Events 

19+0+0+0 

19+0+0+0 

19+0+0+0 

19+0+0+0 

13+0+0+0 

0.080 

0.120 

0.070 

0.096 

0.097 

2a  Factor 

1.669 

2.278 

1.570 

1.748 

1.638 

P 

0.983 

0.957 

0.987 

0.980 

0.984 

D  =  Degelen,  S  =  Shagan  (Balapan),  M  =  Murzhik  (Konystan). 
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Table  4C. 

Maximum-Likelihood  Yield  Estimates  of  E.  Kazakh  Shots 

Event,  Region 

Official  W 

ISC 

NEIS 

Sykes 

Marshall 

TG,  Pmax 

[KT] 

[KT] 

[KT] 

[KT] 

[KT] 

[KT] 

651121,  D 

29.0 

31.7 

43.9 

30.8 

27.3 

25.4 

660213,  D 

125.0 

169.5 

179.1 

160.8 

196.8 

162.9 

660320,  D 

100.0 

121.2 

179.1 

115.5 

102.1 

89.5 

660507,  D 

4.0 

2.2 

1.3 

2.2 

1.9 

2.7 

670922,  M 

10.0 

8.3 

7.6 

8.2 

7.1 

680929,  D 

60.0 

62.0 

43.9 

59.6 

59.3 

50.4 

690723,  D 

16.0 

16.2 

15.3 

15.9 

20.1 

15.4 

691130,  S 

125.0 

121.2 

88.6 

99.2 

104.7 

123.6 

691228,  M 

40.0 

44.3 

30.9 

42.9 

48.0 

710425,  D 

90.0 

86.7 

62.4 

94.8 

113.9 

70.2 

710606,  M 

16.0 

22.7 

15.3 

20.7 

21.5 

711009,  M 

12.0 

H  .6 

10.7 

12.2 

13.4 

711021,  M 

23.0 

22.7 

21.7 

22.9 

25.3 

720210. S 

16.0 

16.2 

15.3 

14.4 

13.4 

21.9 

720328,  D 

6.0 

5.9 

5.3 

6.7 

7.4 

9.0 

720816,  D 

8.0 

4.2 

5.3 

6.5 

6.0 

7.4 

720902,  M 

2.0 

3.0 

3.7 

2.8 

2.3 

721102,  S 

165.0 

169.5 

179.1 

210.2 

178.5 

213.6 

721210,  S 

140.0 

121.2 

88.6 

111.4 

89.4 

136.4 

#  of  Events 

19+55+0+17 

19+55+0+17 

19+55+0+17 

19+55+0+17 

13+3+0+5 

^MLE  (mb) 

0.076 

0.113 

0.063 

0.087 

0.093 

2g  Factor 

1.660 

2.222 

1.520 

1.696 

1.617 

P 

0.993 

0.985 

0.996 

0.994 

0.988 

D  =  Degelen,  S  =  Shagan  (Balapan),  M  =  Murzhik  (Konystan). 
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Table  4D.  Expected  mb  of  Eastern  Kazakhstan  Explosions 

mb  y  Curve 

#  of  Events 

10KT 

50  KT 

100KT 

150KT 

ISC 

19+55+0+17 

5.256 

5.736 

5.943 

6.064 

NEIS 

19+55+0+17 

5.380 

5.837 

6.034 

6.150 

Sykes 

19+55+0+17 

5.260 

5.747 

5.956 

6.079 

Marshall 

19+55+0+17 

5.274 

5.805 

6.033 

6.167 

TG,  Pa 

12+3+0+5 

5.109 

5.388 

5.551 

TG,  Pb 

13+3+0+5 

4.761 

5.407 

5.685 

5.848 

TG,  Pmax 

13+3+5+4 

5.626 

5.895 

6.052 

11.5  SHAGAN  RIVER  TEST  SITE 


Table  5A.  Nuclear  Explosions  at  Shagan  River  (Balapan)  Region 


Date 

Lat 

Long 

Depth 

Yield 

Rock 

[N] 

[E] 

[m] 

[KT] 

650115 

49.9350 

79.0094 

178 

100-150 

Sa 

680619 

49.9803 

78.9855 

316 

<20 

Sa 

691130 

49.9243 

78.9558 

472 

125 

Co 

710630 

49.9460 

78.9805 

217 

<20 

Co 

720210 

50.0243 

78.8781 

295 

16 

Al 

721102 

49.9270 

78.8173 

521 

165 

Al 

721210 

50.0270 

78.9956 

478 

140 

TS 

Sa  =  Sandstone,  Al  =  Aleurolite  (Siltstone), 

Co  =  Conglomerate,  TS  =  Tuffaceous  Sandstone 
[from  Bocharov  et  al.  (1989)  and  Vergino  (1989)] 
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Table  5B.  Reported  mb  of  Shagan  River  Explosions 


Date 

ISC 

NEIS 

Sykes 

Marshall 

Stewart* 

Stewart* 

Stewart* 

mb 

mb 

171 

mb 

mb  \ 

log(4*J 

M0 

650115 

5.8 

6.3 

5.905 

5.931 

5.96 

3.87 

15.80 

680619 

ea 

5.5 

5.350 

5.354 

5.60 

3.31 

15.24 

691130 

6.0 

6.0 

5.954 

6.048 

6.14 

4.00 

15.93 

710630 

5.2 

5.4 

5.290 

5.027 

5.29 

2.98 

14.91 

720210 

WSM 

5.5 

5.370 

5.370 

5.58 

3.22 

15.15 

721102 

6.1 

6.2 

6.181 

6.224 

6.39 

4.38 

16.31 

721210 

6.0 

6.0 

5.989 

5.996 

6.06 

4.38 

16.31 

*)  Averaged  over  EKA,  VKA,  GBA,  and  WRA. 


Table  5B.  Reported  mb  of  Shagan  River  Explosions  (Continued) 

Date 

EKA 

Nuttli 

Ringdal 

TG 

TG 

TG 

/7?b 

mb{Lg) 

RMS  Lg 

™b(Pa) 

mb(Pb) 

^b(Pmax ) 

650115 

5.98 

5.87 

5.950* 

5.495 

5.734 

5.882 

680619 

5.70 

_ 

_ 

4.620 

5.002 

5.263 

691130 

6.30 

6.043 

5.380 

5.770 

5.977 

710630 

5.34 

4.472 

4.768 

5.041 

720210 

5.58 

5.55 

5.4” 

4.805 

5.074 

5.306 

721102 

6.41 

6.04 

6.118 

5.592 

5.940 

6.189 

721210 

6.08 

6.09 

6.095 

_ 

5.786 

6.015 

710425** ***)' 

N/A 

N/A 

5.862 

N/A 

N/A 

N/A 

*)  Inferred  indirectly  from  Nuttli’s  mb(Lg)  =  5.87  (Ringdal  and  Marshall,  1989). 

**)  Low  SNR  for  Lg  phase  (see  text). 

***)  A  Degelen  event  used  in  Ringdal  (1989). 
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Table  5C.  Magnitude:Yield  Calibration  Curve  at  Shagan  River  Area 


#  of  Events 

Magnitude 

<KY) 

Slope 

Intercept 

a(mb) 

2a  Factor 

Method 

4+2+0+1 

ISC 

0 

0.655+0.132 

4.584+0.753 

0.116 

2.267 

MLE-CY 

4+2+0+1 

ISC 

oo 

0.628+0.055 

4.645+0.097 

0.077 

1.753 

MLE-CY 

4+2+0+ 1 

NEIS 

0 

0.722±0.709 

4.621+1.225 

0.189 

3.346 

MLE-CY 

4+2+0+ 1 

NEIS 

oo 

0.614±0.093 

4.807+0.131 

0.131 

2.671 

MLE-CY 

4+2+0+1 

Sykes 

0 

0.720±0.105 

4.481+0.600 

0.095 

1.833 

MLE-CY 

4+2+0+1 

Sykes 

oo 

0.698±0.054 

4.525+0.096 

0.069 

1.577 

MLE-CY 

4+0+0+0 

Marshall 

0 

0.795±0.137 

4.385+0.811 

0.088 

1.669 

LS 

4+0+0+0 

Marshall 

0.1 

0.795+ _ 

4.387+ _ 

_ 

Ericsson 

4+0+0+0 

Marshall 

1 

0.777+ _ 

4.420+ _ 

__ 

Ericsson 

4+0+0+0 

Marshall 

CJ1 

4.439+ _ 

__ 

Ericsson 

4+0+0+0 

Marshall 

oo 

0.767±0.105 

4.441+0.206 

0.087 

1.685 

LS 

4+2+0+1 

Marshall 

0 

0.768+0.089 

4.421+0.510 

0.098 

1.803 

MLE-CY 

4+2+0+1 

Marshall 

oo 

0.741+0.052 

4.476+0.090 

0.076 

1.606 

MLE-CY 

4+2+0+1 

EKA,  mb 

0 

0.705+0.246 

4.724+1 .457 

0.236 

4.654 

MLE-CY 

4+2+0+1 

EKA,  mb 

oo 

0.568+0.104 

4.983+0.181 

0.165 

3.809 

MLE-CY 

4+2+0+1 

Stewart,  mb 

0 

0.667+0.220 

4.738+1 .293 

0.207 

4.170 

MLE-CY 

4+2+0+1 

Stewart,  mb 

oo 

0.570+0.087 

4.926+0.152 

0.138 

3.044 

MLE-CY 

4+2+0+1 

Stewart,  log^J 

0 

1.098+0.124 

1.865+0.467 

0.172 

2.061 

MLE-CY 

4+2+0+1 

Stewart,  log{4/00) 

oo 

0.953+0.135 

2.134+0.237 

0.189 

2.498 

MLE-CY 

4+2+0+1 

Stewart,  M0 

0 

1.099+0.124 

13.795+1.942 

0.172 

2.061 

j 

MLE-CY 

4+2+0+ 1 

Stewart,  M0 

oo 

0.953+0.135 

14.064+0.237 

0.189 

2.497 

MLE-CY 

3+0+0+1 

Nuttli,  mb{Lg) 

0 

0.587+0.346 

4.742+2.035 

0.160 

3.516 

MLE-CY 

3+0+0+1 

Nuttli,  mb(Lg) 

oo 

0.546+0.106 

4.835+0.207 

0.087 

2.085 

MLE-CY 

4+0+0+ 1 

0 

1.075+0.123 

3.768+0.741 

0.026 

1.119 

MLE-CY 

4 +0+0+1 

oo 

1.025+0.134 

3.873+0.281 

0.027 

1.130 

MLE-CY 

*)  Degelen  event  710425  was  used  instead  of  Shagan  event  720210. 
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Table  5C.  Magnitude:Yield  Calibration  Curve  at  Shagan  River  Area  (Continued) 


#  of  Events 

Magnitude 

0 <rn0) 
°<Y) 

Slope 

Intercept 

o(mb) 

2o  Factor 

Method 

3+2+0+1 

TG,  Pa 

0 

0.759±0.075 

3.873±0.382 

0.082 

1.640 

MLE-CY 

3+2+0+1 

TG,  Pa 

oo 

0.738±0.041 

3.910±0.069 

0.060 

1.456 

MLE-CY 

4+2+0+1 

TG,  Pb 

0 

0.812±0.044 

4.083±0.238 

0.051 

1.332 

MLE-CY 

4+2+0+1 

TG,  Pb 

oo 

0.803±0.028 

4.101±0.050 

0.041 

1.264 

MLE-CY 

4+2+0+1 

TG,  Pmax 

0 

0.81 1  ±0.074 

4.298±0.422 

0.082 

1.593 

MLE-CY 

4+2+0+1 

TG,  Pmax 

oo 

0.788±0.039 

4.336±0.068 

0.067 

1.475 

MLE-CY 

/n^(NEIS)  are  biased  high  at  low  yields  for  Shagan  explosions  simply  because  NEIS 
averages  the  signals  reported.  Consequently  their  mb  vs.  log(W)  slope  is  underestimated, 
which  in  turn  causes  the  yields  of  the  high-yield  explosions  to  be  overestimated.  The  yields 
estimated  by  Geotech's  Pb  seem  to  have  accuracy  at  least  as  good  as  that  based  on  . 

In  Tables  5B  and  5C,  Stewart’s  mb,  log^) ,  and  M0  are  those  which  are  averaged  over 
four  arrays:  Eskdalemuir  (EKA)  Scotland,  Yellowknife  (YKA)  Canada,  Gauribidanur  (GBA) 
India,  and  Warramunga  (WRA)  Australia.  The  scatter  is  slightly  reduced  as  compared  to  that 
based  on  a  single  array  EKA.  Marshall’s  mb  values  are  based  on  the  ISC  bulletin  recordings 
(Marshall,  personal  communication). 

Apparently  the  RMS  Lg  averaged  over  the  bandpassed  multi-channel  signals  recorded  at 
NORSAR  fit  the  announced  yields  very  well.  However,  more  data  may  be  needed  to  further 
quantify  its  performance  (Ringdal  and  Hansen,  1989)  ( cf .  Table  5D).  If  Shagan  event  720210 
(which  had  poor  Lg  SNR  at  NORSAR)  is  included,  the  results  would  show  a  slightly  greater 
scatter  (a  =  0.056  and  0.040  for  cases  0  and  respectively).  In  Tables  5C  through  5E,  we 
have  excluded  this  event  at  Ringdal’s  suggestion  (Ringdal,  personal  communication). 

Zavadil  and  Eisenhauer  conjectured  that  the  first  or  “b”  phase  could  replace  the  phase 
"max.’’  However,  these  AFTAC  researchers  and  many  others  did  not  find  convincing  evidence 
to  support  their  argument  (DARPA,  1981).  It  seems  this  conjecture  could  well  be  valid  at 
least  for  Shagan  River.  Among  the  three  phases  we  measured,  the  phase  "b”  has  the  smal¬ 
lest  scatter  (cf.  Table  5C),  and  it  gives  the  best  yield  estimates  (cf.  Table  5D).  At  NTS  and 
Sahara,  the  phase  *‘b”  has  precision  much  better  than  the  phase  "a”.  At  Degelen  Mountain, 
phase  "a”  shows  the  smallest  scatter,  possibly  because  phases  "b”  and  "max"  are  severely 
contaminated  by  the  scattering  at  the  free-surface  topography  (cf.  Table  6C). 


-  37  - 


nv-Yield  Calibration  Curves 

O 


Table  5D.  Maximum-Likelihood  Yield  Estimates  of  Shagan  Explosions 


Date  &  Official  W 


[KT] 


650115,  100-150 


680619,  <20 


691130,  125 


710630,  <20 


720210,  16 


721102,  165 


721210,  140 


dMLE  imb) 


2 a  Factor 


Sykes 

Marshall 

[KT] 

[KT] 

94.9 

92.0 

15.2 

15.3 

111.6 

132.4 

12.5 

5.5 

16.5 

16.1 

235.9 

226.0 

125.2 

112.7 

0.069 

0.076 

1.577 

1.606 

0.989 

0.991 

[KT] 


65.0 


15.2 


134.5 


.3 


14.0 


369.2 


97.4 


0.138 


3.044 


0.957 


69.2  269.6 


16.0 


208.1  185.3 


144.2  87.6 


13.8 


227.6 


227.6 


0.189 


2.498 


0.962 


13.8 


227.6 


227.6 


0.189 


2.497 


0.963 


0.077  0.131 


1.753  2.671 


0.986  0.958 


Table  5D.  Maximum-Likelihood  Yield  Estimates  of  Shagan  Explosions  (Continued) 


Date  &  Official  W 

EKA 

Nuttli* 

Nuttli*’ 

Ringdal*** 

[KT] 

[KT] 

[KT] 

[KT] 

[KT] 

max 


[KT] 


89.3 


14.6 


117.8 


.6 


16.6 


218.9 


131.7 


650115,  100-150 

56.9 

680619,  <20 

18.3 

691130,  125 

208.5 

710630.  <20 

4.2 

720210,  16 

11.2 

721102,  165 

325.8 

721210,  140 

85.4 

106.3 


155.0 


147.2 


140.6 

108.2 

9.2 

13.3 

98.2 

119.9 

5.8 

6.8 

16.3 

16.3 

190.3 

195.3 

0.087 

0.027 

0.060 

2.085 

1.130 

1.456 

0.965 

0.979 

0.996 

125.6 


0.041 


dMLE  (mb) 


2o  Factor 


0.165 


3.809 


0.939 


*)  Nuttli  (1986b):  mb(Lg)  =  4.307[±0.067]  +  0.765[±0.027]log(W)  for  5.2  <  mb{Lg)  <  6.7. 

**)  Nuttli’s  mb{Lg)  regressed  by  our  maximum-likelihood  code. 

***)  Degelen  event  710425  was  used  instead  of  Shagan  event  720210. 


0.067 


1.475 


0.994 
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In  1988  the  United  States  and  the  Soviet  Union  signed  a  bilateral  agreement  whereby 
each  country  was  permitted  to  monitor  at  close  distance  an  underground  nuclear  explosion  at 
the  others  main  test  site.  The  Soviet  JVE  (Joint  Verification  Experiment)  shot  was  detonated 
on  September  14,  1988,  near  the  southern  edge  of  the  Shagan  River  Test  Site.  The  New 
York  Times  states  that  the  American  and  Soviet  on-site  measurements  are  said  to  give  yields 
of  115KT  and  122KT,  respectively,  for  the  Soviet  JVE  explosions  (Sykes  and  Ekstrom,  1989). 
NORSAR’s  RMS  Lg  measurement  for  this  event  was  5.969  (Ringdal  and  Marshall,  1989). 
Assuming  that  the  actual  yield  was  between  100  and  150KT,  as  suggested  by  P.  G.  Richards, 
the  regression  using  NORSAR’s  RMS  Lg  data  including  this  event  would  give  an  estimate  of 
111.2KT,  which  is  very  close  to  Sykes'  113KT  based  on  the  average  of  mb  and  Ms  (Sykes 
and  Ekstrom,  1989).  The  o(mb)  and  the  95%  factor  in  yield  associated  with  NORSAR’s  data 
reduce  from  0.027  and  1.130  (cf.  Table  5D)  to  0.026  and  1.122,  respectively. 


Table  5E.  Expected  Magnitudes  of  Shagan  Explosions 

mb  :Y  Curve 

#  of  Events 

10KT 

50  KT 

100KT 

150KT 

ISC 

4+2+0+1 

5.273 

5.711 

5.900 

6.011 

NEIS 

4+2+0+1 

5.421 

5.850 

6.035 

6.144 

Sykes 

4+2+0+ 1 

5.223 

5.711 

5.921 

6.044 

Marshall 

4+2+0+ 1 

5.217 

5.735 

5.958 

6.088 

EKA,  mb 

4+2+0+ 1 

5.551 

5.948 

6.119 

6.219 

Stewart,  mb 

4+2 +0+1 

5.496 

5.895 

6.067 

6.167 

Stewart,  logf'f'J 

4+2+0+ 1 

3.087 

3.753 

4.040 

4.208 

Stewart,  Ma 

4+2+0+1 

15.017 

15.683 

15.970 

16.138 

Nuttli,  mb(Lg) 

3+0+0+1 

5.381 

5.762 

5.926 

6.023 

4+0+0+1 

_ 

5.614 

5.923 

6.103 

TG,  Pa 

3+2+0+1 

4.648 

5.164 

5.386 

5.516 

TG  ,Pb 

4+2+0+1 

4.904 

5.465 

5.707 

5.848 

TG,  Pmax 

4+2+0+1 

5.133 

5.684 

5.922 

6.062 

In  Table  5D,  we  have  listed  two  sets  of  yield  estimates  based  on  Nuttli’s  (1986b)  mb(Lg) 
measurements.  Although  Nuttli’s  mb{Lg)  database  for  Shagan  had  only  four  events  in  com¬ 
mon  with  that  of  Bocharov,  it  is  clear  that  regressing  the  mb{Lg)  (or  mb)  on  each  test  site 
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separately,  whenever  the  data  are  available,  would  give  better  results  than  the  global  calibra¬ 
tion  curve  as  recommended  by  Nuttli. 

Table  5E  indicates  that  our  mb(Pmax )  matcnes  Ringdal's  RMS  Lg  very  well  (except  at 
low  yields).  Note  that  mb(Lg)  is  defined  to  be  equal  to  mb  in  eastern  North  America,  which 
has  geology  similar  to  Eastern  Kazakhstan.  Thus  relative  to  mb(Lg)  scaling,  our  mb(Pmax) 
seem  to  have  the  smallest  bias,  as  compared  to  other  mb.  We  have  also  regressed  various 
magnitudes  on  Ringdal’s  RMS  Lg  with  slope  fixed  at  1  (Table  5F).  As  expected,  our  mb{Pb) 
and  mb(Pmax)  possess  the  strongest  correlation,  the  smallest  scatter  around  the  fitted  straight 
line,  as  well  as  the  smallest  standard  error  in  the  estimated  intercept. 


Table  5F.  Various  Magnitudes  Versus  Ringdal’s  RMS  Lg  for  Shagan  Events* 

Magnitude 

a(mb) 

Intercept 

P 

ISC 

0.068 

-0.034±0.028 

0.968 

NEIS 

0.143 

0.066±0.058 

0.840 

Sykes 

0.094 

0.002±0.039 

0.960 

Marshall 

0.112 

0.030±0.050 

0.926 

EKA,  mb 

0.136 

0.149±0.068 

0.907 

Stewart,  mb 

0.125 

0.105±0.062 

0.912 

Stewart,  log(*FJ 

0.211 

-1.951±0.106 

0.955 

Stewart,  M0 

0.211 

9.979±0.106 

0.955 

Nuttli,  mb(Lg) 

0.085 

-0.016±0.049 

0.986 

TG,  Pa 

0.077 

-0.560±0.039 

0.968 

TG,  Pb 

0.058 

-0.266±0.026 

0.988 

TG,  Pmax 

0.064 

-0.057*0.029 

0.982 

*)  Regressed  on  RMS  Lg  with  slope  1  and  free  intercept. 


It  should  not  be  surprising  that  M0  and  log('F„)  reported  by  the  four  British  arrays  give 
identical  slope,  o,  p,  and  yield  estimates  etc.  (Tables  5C,  5D,  and  5F)  since  Stewart  (1988) 
computed  the  seismic  moment,  M0 ,  as 

M0  =  4  n  d  Vp2 
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where  d  =  2.4  g/cc  and  Vp  =  5.0  km/sec  are  the  presumed  density  and  the  P-wave  velocity  of 
the  source  material. 


11.6  DEGELEN  MOUNTAIN 


Table  6A.  Nuclear  Explosions  at  Degelen  Mountainous  Region 


1 - 

Date 

Lat 

Long 

Depth 

Yield 

Rock 

[N] 

!.   ra..  

[m] 

[KT] 

611011 

49.77272 

77.99500 

116 

<20 

Gr 

620202 

49.77747 

78.00164 

238 

<20 

Gr 

640315 

49.81597 

78.07517 

220 

20-150 

Gr 

640516 

49.80772 

78.10197 

253 

20-150 

Gr 

640719 

49.80908 

78.09292 

168 

<20 

Gr 

641116 

49.80872 

78.13344 

194 

20-150 

QP 

650303 

49.82472 

78.0526 7 

196 

<20 

Gr 

65051 1 

49.77022 

77.99428 

103 

<20 

Gr 

650617 

49.82836 

78.06686 

152 

<20 

Gr 

650729 

49.77972 

77.99808 

126 

<20 

Gr 

650917 

49.81158 

78.14669 

156 

<20 

QP 

651008 

49.82592 

78.11144 

204 

<20 

QP 

651121 

49.81919 

78.06358 

278 

29 

Gr 

651224 

49.80450 

78.10667 

213 

<20 

QP 

660213 

49.80894 

78.12100 

297 

125 

QP 

660320 

49.76164 

78.02389 

294 

100 

QP 

660421 

49.80967 

78.10003 

178 

<20 

Gr 

660507 

49.74286 

78.10497 

274 

4 

QP 

660629 

49.83442 

78.07336 

187 

20-150 

Gr 

660721 

49.73667 

78.09703 

170 

<20 

QP 

Gr  =  Granite.  QP  =  Quartz  Porphyrite,  Po  =  Porphyrite,  QS  =  Quartz  Syenite 
[from  Bocharov  et  al.  (1989)  and  Vergino  (1989)] 
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Table  6A.  Nuclear  Explosions  at  Degelen  Mountainous  Region  (Continued) 


Date 

Lat 

Long 

Depth 

Yield 

Rock 

[N] 

[E] 

[m] 

[KT] 

660805 

49.76431 

78.04242 

171 

<20 

Gr 

660819 

49.82708 

78.10875 

134 

<20 

QP 

660907 

49.82883 

78.06375 

117 

<20 

Gr 

661019 

49.74711 

78.02053 

185 

20-150 

Gr 

661203 

49.74689 

78.03336 

153 

<20 

Gr 

670130 

49.76744 

77.99139 

| 

<20 

QS 

670226 

49.74569 

78.08231 

20-150 

QP 

670325 

49.75361 

78.06300 

152 

<20 

Gr 

670420 

49.74161 

78.10542 

225 

20-150 

QP 

670528 

49.75642 

78.01689 

262 

<20 

QP 

670629 

49.81669 

78.04903 

195 

<20 

Gr 

670715 

49.83592 

78.11817 

161 

<20 

QP 

670804 

49.76028 

78.05550 

160 

<20 

Gr 

671017 

49.78089 

78.00383 

181 

20-150 

Gr 

671030 

49.79436 

78.00786 

173 

<20 

Gr 

671208 

49.81714 

78.16378 

150 

<20 

QP 

680107 

49.75442 

78.03094 

237 

<20 

Gr 

680424 

49.84519 

78.10322 

127 

<20 

QP 

68061 1 

49.79300 

78.14508 

149 

<20 

QP 

680712 

49.75469 

78.08994 

172 

<20 

Gr 

680820 

49.82264 

78.07447 

208 

<20 

Gr 

680905 

49.74161 

78.07558 

162 

<20 

Gr 

680929 

49.81197 

78.12194 

290 

60 

QP 

681109 

49.80053 

78.13911 

125 

<20 

QP 

681218 

49.74594 

78.09203 

194 

<20 

Gr 

690307 

49.82147 

78.06267 

214 

20-150 

Gr 

Gr  =  Granite,  QP  =  Quartz  Porphyrite,  Po  =  Porphyrite,  QS  =  Quartz  Syenite 
(from  Bocharov  et  al.  (1989)  and  Vergino  (1989)] 


-42- 


mb-Yield  Calibration  Curves 


Table  6A.  Nuclear  Explosions  at  Degelen  Mountainous  Region  (Continued) 

Date  Let  Long  Depth  Yield  Rock 

[N]  [i]  M  [KT] 


690516 

49.75942 

78.07578 

184 

<20 

Gr 

690704 

49.74603 

78.11133 

219 

<20 

QP 

690723 

49.81564 

78.12961 

175 

16 

QP 

69091 1 

49.77631 

77.99669 

190 

<20 

■a 

691001 

49.78250 

78.09831 

144 

<20 

ma 

691229 

49.73367 

78.10225 

86 

<20 

QP 

700129 

49.79558 

78.12389 

214 

20-150 

Po 

700327 

49.74781 

77.99897 

138 

<20 

Gr 

700527 

49.73131 

78.09861 

66 

<20 

QP 

700628 

49.80150 

78.10681 

332 

20-150 

Gr 

700724 

49.80972 

78.12839 

154 

<20 

QP 

700906 

49.75975 

78.00539 

212 

<20 

Gr 

701217 

49.74564 

78.09917 

193 

<20 

Gr 

710322 

49.79847 

78.10897 

283 

20-150 

Gr 

710425 

49.76853 

78.03392 

296 

90 

Gr 

710525 

49.80164 

78.13883 

132 

<20 

Gr 

711129 

49.74342 

78.07850 

<20 

Gr 

711215 

49.82639 

77.99731 

115 

<20 

Gr 

711230 

49.76003 

78.03714 

249 

20-150 

Gr 

720310 

49.74531 

78.11969 

171 

<20 

720328 

49.73306 

78.07569 

124 

6 

El 

720607 

49.82675 

78.11547 

208 

20-150 

QP 

720706 

49.73750 

78.11006 

81 

<20 

QP 

720816 

49.76547 

78.05883 

139 

8 

Gr 

721210 

49.81939 

78.05822 

264 

20-150 

Gr 

721228 

49.73919 

78.10625 

132 

<20 

QP 

Gr  =  Granite,  QP  =  Quartz  Porphyrite,  Po  =  Porphyrite,  QS  =  Quartz  Syenite 
(from  Bocharov  et  al.  (1989)  and  Vergino  (1989)] 
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fable  68.  Reported  mb  of  Degeten  Mountain  Explosions 


ISC 

NEIS 

Sykes 

Marshall 

TG,  Pa 

TG,  Pb 

TG,  Pr 

mb 

mb 

mb 

mb 

mb 

fTlfo 

mb 

mb-Yield  Calibration  Curves 


Table  6B.  Reported  mb  of  Degelen  Mountain  Explosions  (Continued) 


ISC 

NEIS 

Sykes 

Marshall 

TG,  Pa 

TG,  Pb 

TG,  P max 

mb 

mb 

mb 

mb 

mb 

mb 

[Tlfo 

670420 


670528 


670629 


670715 


670804 


671017 


671030 


671208 


680107 


680424 


68061 1 


680712 


680820 


680905 


680929 


681109 


681218 


690307 


690516 


690704 


690723 


690911 


5.500 


5.400 


5.300 


5.400 


5.300 


5.600 


5.300 


5.400 


5.100 


5.000 


5.200 


5.300 


4.800 


5.400 


5.800 


4.900 


5.000 


5.600 


5.200 


5.200 


5.400 


5.000 


5.556 


5.464 


5.336 


5.387 


5.316 


5.629 


5.413 


5.314 


4.977 


4.911 


5.240 


5.169 


4.761 


5.439 


5.861 


4.751 


5.044 


5.664 


5.264 


5.241 


5.127 


5.434 


5.629 


4.596 

4.922 

5.169 

3.977 

4.236 

4.578 
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Table  6B.  Reported  mb  of  Degelen  Mountain  Explosions  (Continued) 
ISC  NEIS  Sykes  Marshall  TG,  Pa  TG.  Pb 
mb  mb  mb  mb  mb  mb 


700724 


700906 


701217 


710322 


710425 


710525 


711129 


711215 


711230 


720310 


720328 


720607 


720706 


720816 


721210 


721228 


5.130 


5.600 


4.900 


TG,  Pa 
mb 


TG-  Pma« 

mb 


15 

4.339 

4.622 

4.887 

5 

4.977 

5.355 

5.534 

Table  6B  indicates  that  for  Degelen  events,  all  other  mb's  are  systematically  larger  than 
ours  by  a  A mb  of  approximately  0.2  to  0.3.  The  mb  offset  is  less  significant  for  Shagan 
events  (cf.  Tables  68  and  5E),  however.  This  is  possibly  due  to  the  different  focusing  and 
defocusing  patterns  between  Shagan-Europe  and  Degelen-Europe  paths.  In  our  WWSSN 
database,  there  were  8  and  10  European  stations  which  detected  the  Shagan  event  650115 
(100-150KT)  and  Degelen  event  710425  (90KT),  respectively.  The  averaged  mb  residuals  of 
the  European  WWSSN  stations  for  these  two  events  are  0.07±0.133  and  0.122±0.057,  respec¬ 
tively.  Marshall’s  database  has  a  heavy  clustering  of  ISC  stations  in  Europe,  and  hence  the 
resulting  mb  offset  may  just  be  reflecting  the  even  more  severe  path  focusing  effects 
enhanced  by  the  ISC  clustering  in  the  western  Europe. 
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Table  6C.  mb: Yield  Calibration  Curve  at  Degelen  Mountain 


#  of  Events 

mb 

<KY) 

Slope 

Intercept 

o(mb) 

2c  Factor 

Method 

9+46+0+15 

ISC 

0 

0.833+0.022 

4.362+0.114 

0.069 

1.466 

MLE-CY 

9+46+0+15 

ISC 

oo 

0.809+0.015 

4.392+0.018 

0.058 

1.392 

MLE-CY 

9+46+0+15 

NEIS 

0 

0.861+0.036 

4.445+0.195 

0.129 

1.997 

MLE-CY 

9+46+0+15 

NEIS 

oo 

0.773+0.024 

4.556+0.028 

0.110 

1.920 

MLE-CY 

9+46+C+15 

Sykes 

0 

0.803+0.019 

4.414+0.102 

0.062 

1.426 

MLE-CY 

9+46+0+15 

Sykes 

oo 

0.786+0.012 

4.438+0.015 

0.051 

1.345 

MLE-CY 

9+0+0+0 

Marshall 

0 

0.912+0.067 

4.300+0.377 

0.097 

1.635 

LS 

9+0+0-*-0 

Marshall 

0.1 

0.912+ _ 

4.300+ _ 

Ericsson 

9+0+0+0 

Marshall 

1 

0.897+ _ 

4.322+ _ 

_ 

Ericsson 

9+0+0+0 

Marshall 

5 

0.885+ _ 

4.338+ _ 

Ericsson 

9+0+0+0 

Marshall 

100 

0.884+ _ 

4.340+ _ 

_ 

Ericsson 

9+0+0+0 

Marshall 

oo 

0.884+0.059 

4.340+0.090 

0.096 

1.647 

LS 

9+46+0+15 

Marshall 

0 

0.908+0.022 

4.318+0.115 

0.083 

1.525 

MLE-CY 

9+46+0+15 

Marshall 

oo 

0.869+0.017 

4.370+0.020 

0.076 

1.494 

MLE-CY 

9+1 +0+3 

TG,  Pa 

0 

0.981+0.048 

3.449+0.226 

0.087 

1.505 

MLE-CY 

£+1+0+3 

TG,  Pa 

oo 

0.959+0.044 

3.479+0.066 

0.084 

1.499 

MLE-CY 

9+1 +0+3 

[ 

TG,  Pb 

0 

0.972+0.058 

3.752+0.297 

0.108 

1.665 

MLE-CY 

9+1 +0+3 

TG,  Pb 

oo 

0.939+0.052 

3.798+0.079 

0.103 

1.654 

MLE-CY 

9+1+0+3 

TG,  Pmax 

0 

0.931+0.062 

4.033+0.333 

0.108 

1.709 

MLE-CY 

9+1+0+3 

TG,  P max 

oo 

0.899+0.051 

4.079+0.078 

0.099 

1.660 

MLE-CY 
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Table  6D.  Expected  mb  of  Degelen  Explosions 


mb: Y  Curve 

#  of  Events 

10KT 

50KT 

ISC 

9+46+0+15 

5.201 

5.767 

NEIS 

9+46+0+15 

5.329 

5.870 

Sykes 

9+46+0+15 

5.223 

5.772 

Marshall 

9+46+0+15 

5.239 

5.846 

TG  ,Pa 

9+1 +0+2 

4.438 

5.108 

TG.  Pb 

9+1 +0+2 

4.737 

5.393 

TQ.  Pmax 

9+1 +0+2 

4.978 

5.606 

100KT 


6.010 


6. 


6. 


6.108 


5.397 


5.676 


5.876 


150KT 


6.147 


6.261 


5.566 


5.841 


6.034 
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11.7  KONYSTAN  (MURZHIK)  AREA 


Table  7 A.  Explosions  at  Konystan  (Murzhik)  Region 

Date 

Lat 

Long 

Depth 

Yield 

Rock 

ISC 

NEIS 

Sykes 

Marshall 

[N] 

[E] 

[m] 

[KT] 

mb 

mb 

mb 

651014 

49.9906 

77.6357 

048 

1.1 

Al 

661218 

49.9246 

77.7472 

427 

20-150 

Po 

5.8 

5.9 

5.800 

5.922 

670916 

49.9372 

77.7281 

230 

<20 

Sa 

5.3 

5.3 

5.300 

5.245 

670922 

49.9596 

77.6911 

229 

10 

Al 

5.2 

5.3 

5.200 

5.160 

671122 

49.9419 

77.6868 

227 

<20 

Al 

4.8 

_ 

4.800 

4.410 

681021 

49.7279 

78.4863 

31 

0.2 

Ar 

_ 

_ 

681112 

49.7124 

78.4613 

31 

0.2x3 

Gs 

690531 

49.9503 

77.6942 

258 

<20 

Al 

5.3 

5.4 

5.300 

5.290 

691228 

49.9373 

77.7142 

388 

46 

Al 

mm 

5.700 

5.791 

700721 

49.9524 

77.6729 

225 

<20 

Sa 

5.4 

5.400 

5.376 

701104 

49.9892 

77.7624 

249 

<20 

Po 

5.4 

5.400 

5.439 

710606 

49.9754 

77.6603 

299 

16 

Al 

5.5 

5.5 

5.480 

5.526 

710619 

49.9690 

77.6408 

290 

<20 

Po 

mu 

5.5 

5.410 

5.538 

711009 

49.9779 

77.6414 

237 

12 

Al 

5.3 

5.4 

5.320 

5.371 

711021 

49.9738 

77.5973 

324 

23 

Sa 

5.5 

5.6 

5.510 

5.580 

720826 

49.9820 

77.7166 

285 

<20 

Al 

5.3 

5.5 

5.370 

5.363 

720902 

49.9594 

77.6409 

185 

2 

Sa 

4^_ 

5.1 

4.880 

4.788 

Sa  =  Sandstone,  Al  =  Aleurolite  (Siltstone),  Po  =  Porphyrite,  Gs  =  Gritstone 
[from  Bocharov  et  al.  (1989)  and  Vergino  (1989)) 
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Table  7B.  mb:Yield  Calibration  Curve  at  Konystan  Area 

#  of  Events 

mb 

o^o) 

<*Y) 

Slope 

Intercept 

<*(mb) 

2o  Factor 

Method 

6+7+0+1 

ISC 

0 

0.632±0.097 

4.659+0.518 

0.093 

1.962 

MLE-CY 

6+7+0+1 

ISC 

oo 

0.602+0.036 

4.691+0.042 

0.057 

1.542 

MLE-CY 

6+7+0+ 1 

NEIS 

0 

0.500+0.106 

4.894+0.573 

0.099 

2.490 

MLE-CY 

6+7+0+1 

NEIS 

oo 

0.472+0.023 

4.920+0.027 

0.046 

1.562 

MLE-CY 

6+7+0+1 

Sykes 

0 

0.638+0.080 

4.650+0.426 

0.077 

1.740 

MLE-CY 

6+7+0+ 1 

Sykes 

oo 

0.617+0.031 

4.671+0.036 

0.048 

1.429 

MLE-CY 

6+0+0+0 

Marshall 

0 

0.791+0.102 

4.498+0.547 

0.081 

1.598 

LS 

6+0+0 +0 

Marshall 

0.1 

0.790+ _ 

4.498+ _ 

Ericsson 

6+0+0 +0 

Marshall 

1 

0.772+ _ 

4.519+ _ 

_ 

Ericsson 

6+0+0 +0 

Marshall 

5 

0.761+ _ 

4.531+ _ 

_ 

_ 

Ericsson 

6+0+0 +0 

Marshall 

100 

0.760+ _ 

4.532+ _ 

_ 

_ 

Ericsson 

6+0+0+0 

Marshall 

oo 

0.760+0.077 

4.532+0.091 

0.079 

1.613 

LS 

6+7+0+1 

Marshall 

0 

0.806+0.065 

4.495+0.347 

0.089 

1.666 

MLE-CY 

6+7+0 +1 

Marshall 

oo 

0.768+0.039 

4.535+0.045 

0.069 

1.516 

MLE-CY 

6+7+0 +1 

DOB* 

0 

0.278+0.400 

2.136+0.972 

0.143 

10.713 

MLE-CY 

6+7+0+1 

DOB 

oo 

0.245+0.024 

2.164+0.028 

0.035 

1.915 

MLE-CY 

*)  Depth  of  Burial. 


Our  maximum-likelihood  regression  routine  can  be  applied  to  estimate  the  depth  scaling 
rule  as  well  (Jih,  1990).  The  result  in  Table  7B  indicates  that  the  scale  depth  for  Konystan 
explosions  is  146±1  meters.  Furthermore,  The  depth  of  burial  [DOB]  is  proportional  to  the 
quartic  root  of  the  yield,  rather  than  the  cubic  root  as  frequently  cited  at  NTS  (e.g. ,  Evernden 
and  Marsh,  1987).  For  Konystan  test  site,  the  yields  estimated  using  DOB  seem  to  have 
accuracy  comparable  to  mb  (Table  7C).  This  is  not  the  case  for  Degelen  region,  however  (Jih. 
1990). 
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Table  7C.  Maximum-Likelihood  Yield  Estimates  of  Konystan  Explosions 


Official  W 


[KT] 


ISC 

NEIS 

Sykes 

Marshall 

DOB 

[KT] 

[KT] 

[KT] 

[KT] 

[KT] 

661218 

20-150 

69.5 

670916 

<20 

10.3 

670922 

10 

7.0 

671122 

<20 

1.5 

690531 

<20 

10.3 

691228 

46 

47.4 

700721 

<20 

15.0 

701104 

<20 

15.0 

710606 

16 

22.1 

710619 

<20 

15.0 

711009 

12 

10.3 

711021 

23 

22.1 

720826 

<20 

10.3 

720902 

2.0 

2.2 

dMLE (mb)  or  dMLE(°0B) 


2a  Factor 


0.057 

0.046 

0.048 

0.069 

0.035 

1.542 

1.562 

1.429 

1.516 

1.915 

0.990 

0.993 

0.993 

0.992 

0.971 

Table  7D.  Expected  mb  &  DOB  of  Konystan  Explosions 


mb  :Y  Curve 

#  of  Events 

10KT 

50  KT 

100KT 

150KT 

ISC 

6+7+0+1 

5.293 

5.714 

5.895 

6.001 

NEIS 

6+7+0+1 

5.392 

5.723 

5.865 

5.948 

Sykes 

6+7+0+1 

5.289 

5.720 

5.906 

6.014 

Marshall 

6+7+0+1 

5.302 

5.839 

6.070 

6.205 

DOB  (meter) 

6+7+0+1 

257 

380 

451 

498 
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11.8  CRATERING  VERSUS  NON  CRATERING  EXPLOSIONS 

Tables  8A  and  8B  list  the  A,mb  s  mb(Pmax )  -  mb(Pa )  and  A zmb  &  mb{Pb)  -  mb{Pa)  at 
four  different  test  sites.  As  the  yield  increases  from  10KT  to  150KT,  the  &mb  decreases 
steadily,  except  at  Degelen  Mountain.  This  could  be  yet  another  indication  that  the  D.O.B.  at 
Degelen  does  not  quite  follow  the  depth  scaling.  Note  that  Sahara  Test  Site  has  the  same 
trend  as  Shagan. 


Table  8A.  Expected  mb(Pmax)  -  mb(Pa)  At  4  Test  Sites 

Test  Site 

10KT 

50  KT 

100KT 

150KT 

NTS 

0.486 

0.501 

0.507 

0.510 

Sahara 

0.467 

0.539 

0.569 

0.588 

KTS 

0.541 

0.517 

0.507 

0.501 

Shagan  River 

0.485 

0.520 

0.536 

0.546 

Degelen 

0.540 

0.498 

0.479 

0.468 

Table  8B.  Expected  mb(Pb)  -  mb(Pa)  At  4  Test  Sites 

Test  Site 

10KT 

50KT 

100KT 

150KT 

NTS 

0.232 

0.254 

0.263 

0.268 

Sahara 

0.187 

0.282 

0.322 

0.346 

KTS 

0.299 

0.298 

0.297 

0.297 

Shagan  River 

0.256 

0.301 

0.321 

0.332 

Degelen 

0.299 

0.285 

0.279 

0.275 

McLaughlin  et  al.  (1985)  studied  the  ratio  of  the  Pa  phase  and  Pmax  phase  of  presumed 
Shagan  River  contained  and  cratering  explosions  by  comparing  the  WWSSN  station  mb’s. 
The  motivation  was  that  the  logarithm  of  amplitude  ratio  of  Pmax/Pa  of  event  650115  was 
significantly  smaller  than  other  presumed  contained  explosions  in  the  vicinity.  Assuming  the 
phase  Pa  is  unaffected  by  the  influence  of  the  non-linear  free-surface  interference,  then  an 
adjustment  to  the  mb{Pmax)  should  be  able  to  convert  that  to  a  contained  explosion  of  the 
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same  yield.  McLaughlin  et  at.  (1985)  concluded  that  a  correction  between  0.17  and  0.27  is 
needed  for  this  conversion,  assuming  a  yield  of  125KT. 

Based  on  46  Shagan  River  explosions  recorded  at  EKA,  Ringdal  and  Marshall  (1989) 
derived  a  value  of  0.75  as  their  mean  log(Pmax/Pa)  across  the  EKA  array  using  the  same 
techniques  as  used  in  McLaughlin  et al.  (1985).  The  cratering  event  650115  had  mb(Ptnax)  - 
mb(Pa)  =  0-62  at  EKA,  and  hence  they  apply  a  correction  of  5.87  +  0.75  -  0.62  =  6.00  for  a 
hypothetical  contained  explosion  with  equivalent  yield.  Both  Ringdal  and  Marshall  (1989)  and 
McLaughlin  et  al.  (1985)  have  the  same  methodological  drawback  in  that  they  did  not  take 
the  yields  of  those  reference  contained  explosions  into  account,  due  to  the  lack  of  data  at  the 
time. 


We  utilize  the  statistics  in  Tables  8A  and  8B  to  illustrate  that  the  correction  by  Ringdal 
and  Marshall  (1989)  might  be  slightly  more  accurate  than  that  in  McLaughlin  et  al.  (1985).  In 
Table  8A,  we  have  rnb(Pmax)  -  mb(Pa)  =  0.536  at  100KT,  and  0.546  at  150KT.  Since  for 
event  650115  our  mb{Pmax  ,TG)  -  mb(Pa  ,TG)  =  0.387  (Table  5B),  this  would  imply  an  adjust¬ 
ment  of  0.149  (lOOKT)  and  0.159  (150KT),  and  a  corrected  mb  of  about  6.031  (100KT)  and 
6.041  (150KT),  respectively.  Note  that  the  adjusted  mb  at  150KT,  6.041,  is  almost  identical  to 
the  “expected  mb(Pmax)"  of  6.062  ( cf .  Table  5E).  The  corrected  mb  at  100KT  would  match 
that  of  Ringdal’s  rather  well  if  the  standard  error  in  the  uncorrected  mb(Pmax),  5.882±0.046,  is 
taken  into  account. 

Der  etal.  (1985)  deconvolved  four  contained  and  the  cratering  Shagan  events  [650115] 
recorded  at  EKA,  and  then  they  convolved  the  Green’s  functions  with  an  appropriate  attenua¬ 
tion  operator  as  well  as  the  source-time  function  of  various  yields  of  interest.  By  comparing 
the  phases  Pa  and  Pmax  of  the  synthetics,  they  obtained  a  cratering-to-contained  correction  of 
0.15,  0.15,  and  0.18  at  60,  125,  and  300KT,  respectively.  The  match  with  our  result  is 
remarkably  good.  This  approach  would  seem  very  attractive  if  the  database  can  be  expanded 
to  events  covering  a  wide  range  of  yields  (and  hence  depths)  and  then  the  method  can  be 
applied  to  events  in  the  same  yield  range. 
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11.9  TEST  SITE  BIAS 


Table  9A.  Expected  mb (P)  -  mb{Lg)  at  Various  Test  Sites 


(Earlier  Studies”) 


Test  Site 

Description 

10KT 

50KT 

100KT 

150KT 

NTS 

mb( ISC)  -  Nuttli's  mb(Lg) 

-0.31 

-0.31 

-0.31 

-0.31 

Shagan  River 

mb( ISC)  -  Nuttli’s  mb(Lg) 

0.04 

0.04 

0.04 

0.04 

Degelen 

mb(\SC)  -  Nuttli’s  mb(Lg) 

0.27 

0.27 

0.27 

0.27 

Degelen 

mb  (Sykes)  -  Nuttli's  mb(Lg ) 

0.23 

0.23 

0.23 

0.23 

Novaya  Zemlya 

mb( ISC)  -  Nuttli’s  mb(Lg) 

-0.11 

-0.11 

-0.11 

-0.11 

(This  Study) 


Test  Site 

Description 

10KT 

50  KT 

100KT 

150KT 

NTS 

/^(Marshall)  -  Patton's  mb(Lg) 

-0.434 

-0.380 

-0.356 

-0.343 

NTS 

mb(pmax.  TG)  -  Patton's  mb(Lg) 

-0.555 

-0.485 

-0.454 

-0.437 

Shagan  River 

mb  (Marshall)  -  RingdaPs  RMS  Lg 

0.121 

0.035 

-0.015 

Shagan  River 

mb(P, max.  TG)  -  RingdaPs  RMS  Lg 

_ 

0.070 

-0.001 

-0.041 

Degelen 

mb  (Marshall)  -  RingdaPs  RMS  Lg 

_ 

0.232 

0.185 

0.158 

Degelen 

mb(PmM,  TG)  -  RingdaPs  RMS  Lg 

_ 

-0.008 

-0.047 

-0.069 

*)  Nuttli  (1987,  1988). 


At  Degelen  and  Shagan,  our  results  show  that  the  mb(P)  -  mb(Lg)  has  a  decreasing  ten¬ 
dency  with  increasing  yield,  contrary  to  the  increasing  trend  at  NTS.  Results  based  on 
Marshall's  mb  are  consistent  with  ours. 
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Table  9B.  Expected  Test  Site  Bias 


(Earlier  Studies) 


Test  Sites 

Description 

10KT 

50KT 

100KT 

150KT 

Shagan  -  NTS 

Nuttli  (1987) 

0.35 

0.35 

0.35 

0.35 

Degelen  -  NTS 

Nuttli  (1987) 

0.54 

0.54 

0.54 

0.54 

Shagan  -  Degelen 

Nuttli  (1987) 

-0.19 

-0.19 

-0.19 

-0.19 

Novaya  Zemlya  -  NTS 

Nuttli  (1987) 

0.20 

0.20 

0.20 

0.20 

(This  Study) 

Test  Sites 

Description 

10KT 

50KT 

100KT 

150KT 

Ringdal’s  RMS  Tg(Shagan)  -  mb(Pmax,  TG,  Sahara) 

0.243 

0.253 

0.251 

Ringdal’s  RMS  Lg  (Shagan)  -  mb(Lg){ 

Patton,  NTS) 

— 

-0.073 

0.009 

0.056 

KTS  -  NTS 

mb  (Marshall) 

0.549 

0.498 

0.475 

0.463 

KTS  -  NTS 

mb(Pa,  TG) 

0.344 

0.408 

0.435 

0.451 

KTS  -  NTS 

mb(Pb,  TG) 

0.411 

0.452 

0.469 

0.480 

KTS  -  NTS 

^b  (,Dmax »  TG) 

0.399 

0.424 

0.435 

0.442 

Shagan  -  NTS 

mb  (Marshall) 

0.492 

0.428 

0.400 

0.384 

Shagan  -  NTS 

mb{Pa,  TG) 

0.530 

0.463 

0.433 

0.416 

Shagan  -  NTS 

mb{Pbl  TG) 

0.554 

0.510 

0.491 

0.480 

Shagan  -  NTS 

max «  TG) 

0.529 

0.482 

0.462 

0.452 

Degelen  -  NTS 

mb  (Marshall) 

0.514 

0.539 

0.550 

0.557 

Degelen  -  NTS 

mb(Pa,  TG) 

0.320 

0.407 

0.444 

0.466 

Degelen  -  NTS 

mb(Pb ,  TG) 

0.387 

0.438 

0.460 

0.473 

Degelen  -  NTS 

Wb  {Pmax  i  TG) 

0.374 

0.404 

0.416 

0.424 

Sahara  -  NTS 

mb(Pa,  TG) 

0.083 

0.132 

0.153 

0.165 

Sahara  -  NTS 

mb(Pb,  TG) 

0.038 

0.160 

0.212 

0.243 

Sahara  -  NTS 

^biPmax  •  TG) 

0.063 

0.168 

0.214 

0.240 

KTS  -  Sahara 

mb{Pgi  TG) 

0.261 

0.276 

0.282 

0.286 

KTS  -  Sahara 

mb{Pb.  TG) 

0.373 

0.292 

0.257 

KTS  -  Sahara 

mb(Pmax  ’  TG) 

0.335 

0.254 

0.220 

0.199 
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Table  9B.  Expected  Test  Site  Bias  (Continued) 

(This  Study) 


Test  Sites 

Description 

10KT 

50KT 

100KT 

150KT 

Shagan  -  Sahara 

mb{Pa,  TG) 

0.331 

0.280 

0.251 

Shagan  -  Sahara 

mb(Pb,  TG) 

| 

0.350 

0.279 

0.237 

Shagan  -  Sahara 

0.465 

0.312 

0.247 

0.209 

Degelen  *  Sahara 

mb(Pa,  TG) 

0.237 

0.275 

0.291 

0.301 

Degelen  -  Sahara 

mb{Pb,  TG) 

0.349 

0.278 

0.248 

0.230 

Degelen  -  Sahara 

fflb(P max  > 

0.310 

0.234 

0.201 

0.181 

Shagan  -  Degelen 

hit,  (ISC) 

0.072 

-0.056 

-0.110 

-0.142 

Shagan  -  Degelen 

mb(  NEIS) 

0.092 

-0.020 

-0.068 

-0.095 

Shagan  -  Degelen 

mb  (Sykes) 

0.000 

-0.061 

-0.088 

-0.103 

Shagan  -  Degelen 

mb  (Marshall) 

-0.022 

-0.111 

-0.150 

-0.173 

Shagan  -  Degelen 

mb{Pa,  TG) 

0.210 

0.056 

-0.011 

-0.050 

Shagan  -  Degelen 

m„(Pb,  TG) 

0.167 

0.072 

0.031 

Shagan  -  Degelen 

( Pmax  i  TG) 

0.155 

0.078 

0.046 

mmm 

Konystan  -  Degelen 

mb(ISC) 

0.092 

-0.053 

-0.115 

-0.152 

Konystan  -  Degelen 

mb(NEIS) 

0.063 

-0.147 

-0.238 

-0.291 

Konystan  -  Degelen 

mb  (Sykes) 

0.066 

-0.052 

-0.103 

-0.133 

Konystan  -  Degelen 

mb  (Marshall) 

0.063 

-0.007 

-0.038 

-0.056 

Konystan  -  Shagan 

mb(  ISC) 

0.020 

0.003 

-0.005 

-0.010 

Konystan  -  Shagan 

mb(NEIS) 

-0.029 

-0.127 

-0.170 

-0.196 

Konystan  -  Shagan 

mb  (Sykes) 

0.066 

0.009 

-0.015 

-0.030 

Konystan  -  Shagan 

mb  (Marshall) 

0.085 

0.104 

0.112 

0.117 

The  mb  bias  between  Sahara  and  Degelen  is  interesting  in  that  different  phases  exhibit 
opposite  tendency  of  bias  change  with  yields.  The  bias  determined  with  phase  "a"  increases 
with  yields,  while  that  of  phases  "b"  and  "max''  decrease. 

Based  on  Geotech's  mb:yield  calibration  curves,  Shagan  River  would  have  more  efficient 
P-wave  coupling  than  does  Degelen  River  by  an  offset  of  about  0.155  m.u.  (magnitude  unit) 
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and  0.046  m.u.  at  10KT  and  100KT,  respectively.  This  mb  bias  can  be  explained  by  the  pro¬ 
found  topography  at  Degelen  Mountain  which  could  cause  strong  P-to-S  conversion,  as  illus¬ 
trated  by  the  linear  finite -difference  calculations  (Jih  and  McLaughlin,  1988).  The  bias  value 
currently  used  by  the  U.S.  government  is  intended  to  be  the  most  appropriate  value  for  yields 
near  the  150KT  threshold  of  the  1974  Threshold  Test  Ban  Treaty  (TTBT).  Our  results  in  Table 
9B  provide  a  direct  clues  of  how  different  the  bias  could  be  at  lower  yields. 

A  brief  review  of  earlier  work  on  test  site  bias  may  be  interesting.  Based  on  the  P-wave 
seismograms  of  three  granite  explosions  (PILEDRIVER,  Shagan  680619,  and  Shagan 
710630)  recorded  at  EKA,  Douglas  (1987)  concluded  that  the  Shagan-NTS  bias  is  about  0.5, 
which  is  very  close  to  what  we  got  with  phases  Pb  and  Pmax  (Table  9B).  Stewart  (1988) 
predicted  a  bias  of  0.37  at  mb=  5.0  and  of  0.32  at  mb= 6.5,  based  on  the  mb  averaged  across 
four  arrays:  Eskdalemuir  (EKA)  Scotland,  Yellowknife  (YKA)  Canada,  Gauribidanur  (GBA) 
India,  and  Warramunga  (WRA)  Australia.  His  predicted  bias  is  yield-dependent,  and  it  has  a 
decreasing  trend  with  increasing  yield,  which  is  consistent  with  our  maximum-likelihood  results 
in  Table  9B.  This  tendency  should  not  be  surprising.  Large-yield  explosions  generate 
predominantly  low-frcqucncy  signals  and  low-yield  explosions  are  relatively  richer  in  higher 
frequencies,  so  a  relatively  large  amount  of  energy  is  removed  by  the  attenuation  from  low 
yield  tests,  hence  the  bias  is  greater  for  such  low  yield  explosions.  Furthermore,  the  bias 
between  two  sites  is  made  up  of  more  than  just  the  attenuation  in  the  mantle  beneath  the  test 
site.  A  difference  in  depth  containment  laws  and  up-hole  velocities  between  the  test  sites  can 
have  an  effect  on  the  observed  amplitudes  and  hence  on  the  final  value  of  bias  (Marshall, 
personal  communication).  The  bias  between  any  two  test  sites  should  be  a  sum  of  these 
effects.  Murphy  and  Tzeng  (1982)  estimated  the  bias  by  comparing  signals  recorded  near 
NTS  and  Semipalatinsk  from  Aleutian  Islands  earthquakes.  They  estimated  the  bias  as  0.24 
magnitude  unit.  Priestley  et  al.  (1987)  used  a  similar  approach,  and  they  estimated  the  bias 
as  0.34.  It  should  be  noted,  however,  that  alt  these  earlier  bias  estimates  were  made  before 
the  publication  of  Bocharov  et  at.  (1989)  and  Vergino  (1989). 
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APPENDIX 

GEOTECH’S  MAXIMUM-LIKELIHOOD  NETWORK  mb:  GLM90A 


Short-period  WWSSN  vertical  recordings  (SPZ)  of  body  waves  from  96  nuclear  explo¬ 
sions  detonated  at  the  Semipalatinsk  Test  Site,  Eastern  Kazakhstan,  USSR,  are  being  meas¬ 
ured  and  added  to  our  database  to  determine  the  optimal  network  magnitudes  using  the 
maximum-likelihood  estimator  (MLE),  which  accounts  for  the  effects  of  data  censoring  due  to 
clipping  and  to  noise.1,2  As  of  now,  our  WWSSN  database  has  been  expanded  to  124  events 
(totaling  366  usable  “a”,  ''b",  and  “max"  event  phases)  from  a  variety  of  test  sites.  Only  the 
stations  at  teleseismic  distance  (20  to  95  degrees)  were  used  in  the  network  mb  determina¬ 
tion.  (Therefore,  some  of  the  mb  values  might  be  slightly  different  from  those  in  an  earlier 
report  TGAL-87-05.)  The  8501  good  signals,  5699  noise  measurements,  and  1088  clipped 
recordings  yield  a  dMLE  0.320. 

The  124  events  in  Table  A1  are  grouped  by  test  sites.  The  three  numbers  under  the 
column  "#  of  signals"  represent  the  number  of  signals,  noise,  and  clips  associated  with  the 
Pmax  phase  of  each  event.  Except  for  the  U.S.  and  French  Sahara  explosions  which  have 
specific  code  names,  all  the  remaining  events  are  identified  with  the  dates  and  abbreviated 
test  site  codes  shown  below: 

azg  Azgir,  U.S.S.R. 

pne  “PNE  ',  Urals,  U.S.S.R. 

mek  Murzhik  (Konystan),  E.  Kazakh,  U.S.S.R. 

dek  Degelen  Mountain,  E.  Kazakh,  U.S.S.R. 

sek  Shagan  River  (Balapan),  E.  Kazakh,  U.S.S.R. 

nnz  Northern  Novaya  Zemlya,  U.S.S.R. 

snz  Southern  Novaya  Zemlya,  U.S.S.R. 

tu  Tuamoto  Islands,  France 

raj  Rajasthan,  India 

ch  Lop  Nor,  Sinkiang,  China 


Table  A2  lists  the  station  correction  terms  determined  jointly  along  with  the  network  mb 
values. 


’  Blandford,  R  R  ,  and  R  H  Shumway  (1982)  Magnitude  yield  for  nuclear  explosions  in  granite  at  the  Nevada  Test  Site 
and  Algena  joint  determination  with  station  effects  and  with  data  containing  clipped  and  low-amplitude  signals,  Technical  Report 
VSC  TR  82-12,  Teledyne  Geotech,  Alexandria,  Virginia 

2  Jih,  R  -S  ,  and  R  H  Shumway  (1989)  Iterative  network  magnitude  estimation  and  uncertainty  assessment  with  noisy 
and  clipped  data.  Bull.  Soismo  Soc  Am..  79,  1122-1141 
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Table  A1.  Geotech’s  Maximum-Likelihood  Network  mb 


Event 

#  of  Signals 

™b(Pa) 

mb{Pb) 

^b(Pmax  ) 

a 

ALMENDRO 

26,0,2 

5.730 

6.026 

6.233 

0.060 

BENHAM 

42.1,7 

5.772 

6.103 

6.359 

0.045 

BILBY 

36,3,0 

5.148 

5.404 

5.658 

0.051 

BOURBON 

18,31,0 

4.587 

4.720 

4.904 

0.046 

BOXCAR 

32,0,4 

5.849 

6.189 

6.412 

0.053 

CAMBRIC 

12,34,0 

4.091 

4.340 

4.551 

0.047 

CHANCELLOR 

16,12,1 

4.887 

5.183 

5.338 

0.059 

CHARTREUSE 

31,16,1 

4.884 

5.249 

0.046 

CHATEAUGAY 

17,28,2 

4.478 

4.884 

5.066 

0.047 

CORDUROY 

18,14,0 

4.971 

5.092 

5.287 

0.057 

HANDCAR 

16,33,0 

4.308 

4.495 

4.629 

0.046 

HANDLEY 

41,1,1 

6.062 

6.307 

6.480 

0.049 

HARZER 

31.5,1 

5.011 

5.312 

5.536 

0.053 

KANKAKEE 

24,27,0 

iWl 

4.597 

4.847 

0.045 

MAST 

29,1,0 

5.739 

5.981 

0.058 

NASH 

31,21,0 

4.758 

4.918 

5.149 

0.044 

PILEDRIVER 

38,12,1 

5.194 

5.435 

0.045 

REX 

16,35,1 

msm 

4.376 

4.720 

0.044 

SCOTCH 

38,8,1 

5.079 

5.344 

5.600 

0.047 

CANNIKIN 

49,0,20 

6.663 

6.911 

0.039 

MILROW 

52,0,4 

5.945 

6.195 

6.494 

0.043 

LONGSHOT 

67,4,3 

5.056 

5.428 

5.818 

0.037 

FAULTLESS 

47,1,3 

5.829 

6.157 

6.460 

0.045 

GASBUGGY 

11,37,0 

4.153 

4.412 

4.661 

0.046 

RIO  BLANCO 

15,20,0 

4.068 

4.545 

4.810 

0.054 

RULISON 

9,37,0 

4.108 

4.240 

4.554 

0.047 

SHOAL 

13,27,0 

4.321 

4.455 

4.738 

0.051 

SALMON 

6,33,0 

3.439 

3.974 

4.180 

0.051 
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Table  A1.  Geotech'e  Maximum-Likelihood  Network  mD  (Continued) 


Event 

#  of  Signals 

™b(Pa) 

mb(Pb) 

( Pmax ) 

a 

azg22apr66 

3,10,0 

3.867 

4.101 

4.183 

0.089 

azg22dec7 1 

12.0,2 

5.473 

5.826 

6.164 

0.086 

azg25apr75 

1,16,0 

3.904 

3.944 

0.078 

azg29jul76 

41,5,7 

5.105 

5.579 

5.864 

0.044 

azg30sep77 

21,30,1 

4.049 

4.588 

4.828 

0.044 

azgl  7oct78 

7,0,5 

5.271 

5.724 

6.097 

0.092 

azg18dec78 

9,0,3 

5.374 

5.748 

6.119 

0.092 

azg17jan79 

10,0,4 

5.515 

5.869 

6.153 

0.086 

azg14jul79 

10,0,1 

4.831 

5.371 

5.699 

0.097 

azg24oct79 

3,0,6 

4.848 

5.681 

5.960 

0.107 

pne29aug74 

27,18,0 

3.994 

4.397 

4.722 

0.048 

mek18dec66 

55,9,1 

5.261 

5.493 

5.709 

0.040 

dek21nov65 

48,15,1 

4.875 

5.152 

5.362 

0.040 

dek13feb66 

51,4,10 

5.640 

5.890 

6.082 

0.040 

dek20mar66 

50,9,8 

5.335 

5.624 

5.850 

0.039 

dek07may66 

9,26,1 

3.992 

4.233 

4.486 

0.053 

dek26feb67 

48,9,6 

5.353 

5.597 

5.821 

0.040 

dek29sep68 

50,8,6 

5.125 

5.432 

5.627 

0.040 

dek23jul69 

38,21,1 

4.594 

4.920 

5.167 

0.041 

dekl 1sep69 

19,39,0 

3.975 

4.234 

4.576 

0.042 

dek25apr71 

37,5,0 

5.299 

5.566 

5.756 

0.049 

dek30dec7 1 

16,3,0 

4.982 

5.524 

m 

dek28mar72 

28,17,0 

4.351 

■ESS 

4.959 

0.048 

dek16aug72 

24,23,1 

4.337 

4.620 

4.885 

0.046 

dek10dec72 

30,7,5 

4.975 

5.333 

5.532 

0.049 

dek29mar77 

25,14,0 

4.304 

4.700 

4.981 

0.051 

dek30jul77 

21,16,0 

4.200 

4.604 

4.857 

0.053 

dek26mar78 

25,6,0 

4.948 

5.497 

0.057 

dek22apr78 

21,9,0 

4.466 

4.765 

5.014 

0.058 
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Table  A1.  Geotech's  Maximum-Likelihood  Network  mb  (Continued) 


Event 

#  of  Signals 

mb(Pa ) 

mb(Pb) 

mb  ( Pmax ) 

o 

sek15jan65 

46,1,2 

5.493 

5.732 

5.880 

sek19jun68 

28,3,2 

4.618 

5.000 

5.261 

sek30nov69 

32,0,0 

5.378 

5.767 

5.975 

0.057 

sek30juh71 

31,19,1 

4.470 

4.766 

5.038 

0.045 

sek10feb72 

34,8,2 

4.803 

5.071 

5.304 

0.048 

sek02nov72 

29,1,15 

5.590 

5.938 

6.187 

0.048 

sek10dec72 

29,2,11 

5.784 

6.013 

0.049 

sek23jul73 

38,1,1 

5.753 

5.996 

6.181 

sekl 4dec73 

45,8,6 

5.245 

5.545 

5.770 

0.042 

sek27apr75 

18,1,1 

4.904 

5.242 

5.491 

0.072 

sek04jul76 

14,0,5 

5.598 

5.927 

0.073 

sek07dec76 

17,2.1 

4.961 

5.416 

5.606 

0.072 

sek11jun78 

17,0.1 

5.296 

5.580 

5.889 

gm 

sek15sep78 

30,1,5 

5.431 

5.691 

5.884 

sek23jun79 

38,3,3 

5.615 

5.846 

6.049 

0.048 

sek14sep80 

29,5,6 

5.439 

5.752 

5.987 

0.051 

nnz27oct66 

56,0,14 

6.063 

6.295 

6.436 

0.038 

nnz21oct67 

53,5,3 

5.400 

5.590 

5.765 

0.041 

nnz07nov68 

59,1,5 

5.580 

5.831 

6.025 

0.040 

nnz14oct69 

59,2,7 

5.760 

5.957 

6.129 

nnz14oct70 

35,0,22 

6.424 

6.633 

6.813 

(^29 

nnz27sep71 

23,0,21 

6.259 

6.475 

6.619 

0.048 

nnz28aug72 

32,0,11 

5.989 

6.247 

6.371 

0.049 

nnz12sep73 

23,0,21 

6.347 

6.672 

6.763 

0.048 

nnz29aug74 

25,0,18 

6.126 

6.394 

6.578 

0.049 

nnz21oct75 

23,0,17 

6.095 

6.333 

6.541 

0.051 

nnz23aug75 

28,0,12 

6.112 

6.367 

6.488 

0  051 

nnz20oct76 

25,34,1 

4.031 

4.350 

4.659 

0.041 

nnz01sep77 

26,2,2 

5.099 

5.415 

5.561 

0.058 

nnzIOaug 78 

39,3,18 

5.392 

5.625 

5.856 

0.041 
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Table  A1.  Geotech’s  Maximum-Likelihood  Network  mb  (Continued) 


Event 

#  of  Signals 

mb{Pa) 

mb(Pb) 

mb(Prr\ax ) 

G 

nnzl  1oct80 

42,4,6 

5.181 

5.442 

5.658 

0.044 

nnz01oct81 

43,4,5 

5.226 

5.489 

5.649 

0.044 

nnz18aug83 

30,5,5 

5.321 

5.526 

5.703 

0.051 

nnz25oct84 

22,3,4 

5.154 

5.427 

5.599 

0.059 

snz27sep73 

32,3,1 

5.196 

5.490 

5.729 

0.053 

14,0,24 

6.647 

6.873 

7.092 

0.052 

snz27oc73b 

9,28,0 

_ 

3.999 

4.150 

0.053 

snz27oc73c 

4,34,0 

3.544 

3.886 

3.908 

0.052 

snz02nov74 

12,0,29 

6.497 

6.790 

7.012 

0.050 

snz18oct75 

21,0,21 

6.227 

6.518 

6.834 

0.049 

BERYL 

11,6,0 

4.412 

4.778 

4.985 

0.078 

CORUNDON 

1 1 ,42,0 

3.797 

3.899 

4.212 

0.044 

EMERAUDE 

14,25,0 

4.261 

4.566 

0.051 

GRENAT 

32,32,1 

4.292 

4.494 

4.763 

0.040 

OPALE 

3,51,0 

3.770 

3.855 

3.896 

0.044 

RUBIS 

42,5.0 

4.826 

5.167 

5.429 

0.047 

SAPHIR 

52,5,5 

5.182 

5.464 

5.716 

0.041 

TOURMALINE 

27,39,0 

4.106 

4.427 

4.644 

0.039 

TURQOISE 

1 1 ,55,0 

— 

3.941 

4.221 

0.039 

tu19feb77 

16,28,0 

4.370 

4.622 

0.048 

tu19mar77 

20,6,1 

5.141 

5.438 

5.639 

0.062 

tu24nov77 

33,0,0 

5.051 

5.369 

5.662 

0.056 

tu25jul79 

18,0,0 

5.090 

5.570 

5.864 

0.075 

tu23mar80 

27,14,3 

4.677 

5.105 

5.358 

0.048 

tu19jul80 

38,2,2 

4.891 

5.158 

5.513 

0.049 

tu03dec80 

32,11,0 

4.689 

4.981 

5.331 

0.049 

tu25jul82 

22,13,0 

4.675 

5.034 

5.210 

0.054 

tu19apr83 

22,1.0 

4.993 

5.199 

5.495 

0.067 
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Table  A1.  Geotech’s  Maximum-Likelihood  Network  mb  (Continued) 


#  of  Signals  mb(Pa)  mb(Pb)  mb(Pmax) 


5.150 


4.820 


4.022 


4.325 


mI 


7,23,0 


ch22sep69 

30,12,0 

ch27oct75 

12,24,0 

ch17oct76 

13,33,0 

ch06oct83 

17,13,1 

ch03oct84 

10,12,0 

ch19dec84 

3,10,0 

4.017 


5.455 


5.234 


4.303 


4.742 


4.396 


4.146 


5.029 


4.747 


3.999 


5.785 


5.611 


4.563 


5.133 


a 


0.075 


0.046 


0.058 


0.049 


0.053 


0.047 


0.057 


0.068 


0.089 
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Table  A2.  WWSSN  Station  Corrections 


Code 

#  of  Signals 

Site  Term 

Longitude 

Latitude 

Description 

AAE 

78,93,17 

-0.243±0.023 

38.765556 

9.029166 

Addis  Ababa,  Ethiopia 

AAM 

134,64,6 

0.254±0.022 

-83.656113 

42.299721 

Ann  Arbor,  Michigan 

ADE 

16,25,0 

0.001  ±0.050 

138.708893 

-34.966946 

Adelaide,  Australia 

AFI 

27,65,0 

-0. 143±0.033 

-171.777252 

-13.909333 

Afiamalu,  Samoa  Islands 

AKU 

71,69,0 

-0.093±0.027 

-18.106667 

65.686668 

Akureyri,  Iceland 

ALQ 

99,15,19 

0.039±0.028 

-106.457497 

34.942501 

Albuquerque,  New  Mexico 

ANP 

20,67,0 

-0.327±0.034 

121.516670 

25.183332 

Anpu,  Formosa 

ANT 

41,49,2 

0.056±0.033 

-70.415276 

-23.705000 

Antofagasta,  Chile 

AQU 

66,44,13 

-0.054±0.029 

13.403055 

42.353889 

Aquila,  Italy 

ARE 

83,40,0 

0.1 01  ±0.029 

-71.491280 

-16.462084 

Arequipa,  Peru 

ASP 

1,2,0 

-0.581  ±0.1 85 

133.896667 

-23.683332 

Alice  Springs,  Australia 

ATL 

79,19,2 

0. 164±0.032 

-84.337502 

33.4oJ334 

Atlanta,  Georgia 

ATU 

112,78,16 

0.146±0.022 

23.716667 

37.972221 

Athens  Univ.,  Greece 

BAG 

132,68,8 

-0.028±0.022 

120.579720 

16.410833 

Baguio  City,  Philippine  Islands 

BDF 

10,2.0 

-0.009±0.092 

-47.903332 

-15.663834 

Brasilia  array,  Brazil 

BEC 

45.102,3 

-0.131±0.026 

-64.681114 

32.379444 

Bermuda-Columbia,  Atlantic  Ocean 

BHP 

30,75,0 

-0.176±0.031 

-79.558052 

8.960834 

Balboa  Heights,  Panama 

BKS 

141,65.1 

0.087±0.022 

-122.235001 

37.876667 

Byerly,  California 

BLA 

152,53,12 

0.122±0.022 

-80.420998 

37.211304 

Blacksburg,  West  Virginia 

BOG 

41,76,0 

0.057±0.030 

-74.065002 

4.623055 

Bogota,  Colombia 

BOZ 

44,4,5 

0.238±0.044 

-111.633331 

45.599998 

Bozeman,  Montana 

BUL 

149,34,9 

0.049±0.023 

28.613333 

-20.143333 

Bulawayo,  Rhodesia 

CAR 

92,59,7 

0.154±0.025 

-66.927635 

10.506667 

Caracas,  Venezuela 

CCG 

1.0,0 

-0.186±0.320 

-61.133335 

77.166664 

Camp  Century,  Greenland 

CHG 

97,16,36 

-0.127±0.026 

98.976944 

18.790001 

Chiengmai,  Asia 

CMC 

50,27,0 

-0.140±0.036 

-115.083336 

67.833336 

Copper  Mine,  Canada 

COL 

259,47,28 

0.087±0.018 

-147.793335 

64.900002 

College  Outposta,  Alaska 

COP 

74,101,14 

0. 166±0.023 

12.433333 

55.683334 

Copenhagen,  Denmark 

COR 

0.1 11  ±0.027 

-123.303192 

44.585724 

Corvallis,  Oregon 

CTA 

57,16,4 

0.214+0.036 

146.254440 

-20.088333 

Charters  Towers,  Australia 
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Table  A2.  WWSSN  Station  Corrections  (Continued) 


Code 

#  of  Signals 

Site  Term 

Longitude 

Latitude 

Description 

DAG 

21,13,5 

0.001±0.051 

-18.770000 

76.769997 

Danmarkshavn,  Greenland 

DAL 

17.24,4 

0.191±0.048 

-96.783890 

32.846111 

Dallas,  Texas 

DAV 

24,81,0 

-0.276±0.031 

125.574722 

7.087778 

Davao,  Philippine  Islands 

DUG 

170,17,29 

0.075+0.022 

-112.813332 

40.195000 

Dugway,  Utah 

EIL 

25,3,43 

0.075±0.038 

34.950001 

29.549999 

Eilat,  United  Arab  Republic 

EPT 

29,2,2 

0.005±0.056 

-106.505836 

31.771667 

El  Paso,  Texas-Mexico  border 

ESK 

88,80,2 

0.1 17±0.025 

-3.205000 

55.316666 

Eskdalemuir,  United  Kingdom 

FLO 

80,20,9 

0.064±0.031 

-90.370003 

38.801666 

Florissant,  Missouri 

FVM 

44,8.0 

-0.008±0.044 

-90.426003 

37.984001 

French  Village,  Missouri 

GDH 

154,126,1 

-0.159±0.019 

-53.533333 

69.250000 

Godhavn,  Greenland 

GEO 

88,69,2 

0.021±0.025 

-77.066666 

38.900002 

Georgetown,  Virginia 

GIE 

9,38,0 

-0.188±0.047 

-90.300003 

-0.733333 

Galapagos  Islands 

GOL 

157,24,11 

-0.216±0.023 

-105.371109 

39.700279 

Golden,  Colorado 

GRM 

1,20,0 

-0.093±0.070 

26.573334 

-33.313332 

Grahamstown,  South  Africa 

GSC 

89,22,16 

0.089±0.028 

-116.804611 

35.301666 

Goldstone,  California 

GUA 

78,175.0 

-0.250+0.020 

144.911667 

13.538333 

Guam,  Mariana  Islands 

HKC 

85,84,0 

-0.131±0.025 

114.171890 

22.303556 

Hong  Kong 

HLW 

47,36,32 

-0.047±0.030 

31.341667 

29.858334 

Helwan,  United  Arab  Republic 

HNR 

30,92,0 

0.188±0.029 

159.947113 

-9.432195 

Honiara,  Solomon  Islands 

HON 

6,9,0 

0.051  ±0.083 

-158.008331 

21.321667 

Honolulu,  Hawaii 

HOW 

I  .... 

1,10,0 

0.258+0.097 

88.309166 

22.416666 

Howrah,  India-Bangladesh  border 

1ST 

102,79,25 

0.186+0.022 

28.995832 

41.045555 

Istanbul,  Turkey 

JCT 

59.4,24 

0.159+0.034 

-99.802223 

30.479445 

Junction  City,  Texas 

JER 

89,45,25 

0.039±0.025 

35.197224 

31.771944 

Jerusalem,  Dead  Sea 

KBL 

14,0,46 

0.142±0.041 

69.043167 

34.540833 

Kabul,  Afghanistan 

55,40,0 

-0.1 81  ±0.033 

11.923889 

78.917503 

Kingsbay,  Svalbard 

121,102,4 

-0.123+0.021 

27.006666 

69.755280 

Kevo,  Finland 

KIP 

84,153,0 

0.107±0.021 

-158.014999 

21.423334 

. | 

Kipapa,  Hawaii 

KOD 

107,33,30 

0.100±0.025 

77.466667 

10.233334 

Kodaikanal,  India 

KON 

129,65,70 

0.102±0.020 

9.598222 

59.649082 

Kongsberg,  Norway 
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Table  A2.  WWSSN  Station  Corrections  (Continued) 


Code 

#  of  Signals 

Site  Term 

Longitude 

Latitude 

Description 

KRK 

8,7.0 

-0.171±0.083 

30.062500 

69.724167 

Kirkenes,  Norway-USSR  border 

KTG 

72,75,1 

-0.247±0.026 

-21.983334 

70.416664 

Kap  Tobin,  Greenland 

LAH 

5,17,3 

0.428±0.064 

74.333336 

31.549999 

Lahore,  India-Pakistan  border 

LEM 

54,82,0 

-0.531  ±0.027 

107.616669 

-6.833333 

Lembang,  Java 

LON 

162,44,21 

-0.034±0.021 

-121.809998 

46.750000 

Longmire,  Washington 

LOR 

74,8,16 

0.154±0.032 

3.851389 

47.266666 

Lormes,  France 

LPA 

8,91,0 

0.426±0.032 

-57.931946 

-34.908890 

La  Plata,  Uruguay 

LPB 

58,39,3 

-0.043*0.032 

-68.098358 

-16.532667 

La  Paz,  Peru-Bolivia  border 

LPS 

50,27,3 

-0.071  ±0.036 

-89.161942 

14.292222 

La  Palma,  Quatemala 

LUB 

40,30,3 

0.214±0.037 

-101.866669 

33.583332 

Lubbock,  Texas 

ME 

0.055±0.027 

-4.411111 

36.727501 

Malaga,  straits  of  Gibraltar 

0.316±0.048 

121.076859 

14.662000 

Manila,  Philippine  Islands 

MAT 

145,53,25 

-0.1 12±0  021 

138.206665 

36.541668 

Matsushiro,  Japan 

MDS 

39,19,0 

-0.032±0.042 

-89.760002 

43.372223 

Madison,  Wisconsin 

MHI 

5.2,2 

0.358±0.107 

59.494499 

36.299999 

Meshed,  Iran-USSR  border 

MNN 

8.6,2 

0. 1 79±0.080 

-93.190002 

44.914444 

Minneapolis,  Minnesota 

MSH 

31,19,9 

0.226±0.042 

59.587776 

36.311111 

Meshed,  Iran-USSR  border 

MSO 

46,7,2 

0.061  ±0.043 

-113.940552 

46.829166 

Missoula,  Montana 

MUN 

39,41,0 

0.015±0.036 

116.20833b 

-31.978333 

Mundaring,  Australia 

NAI 

115,46,9 

-0.089±0.025 

36.803665 

-1.273944 

Nairobi,  Kenya 

NAT 

27,27,0 

0.070±0.044 

-35.033333 

-5.116667 

Natal,  Brazil 

NDI 

129,24,25 

0.124±0.024 

77.216667 

28.683332 

New  Delhi,  India 

NHA 

12,3,0 

-0.127±0.083 

109.211670 

12.210000 

Nhatranga,  Asia 

NIL 

14,6,18 

-0.008+0.052 

73.251663 

33.650002 

Nilore,  Pakistan 

NNA 

47,57,0 

-0.162±0.031 

-76.842140 

-11.987556 

Nana,  Peru 

NOR 

78,50,3 

-0.260±0.028 

-16.683332 

81.599998 

Nord,  Greenland 

NUR 

103,82,7 

0.051  ±0.023 

24.651417 

60.508999 

Nurmijarvi,  Finland 

OGD 

135,63,6 

-0.119+0.022 

-74.595833 

41.087502 

Ogdensburg,  New  York 

OXF 

79,10,17 

0.347±0.031 

-89.409164 

34.511806 

Oxford,  Mississippi 

PDA 

31,103,3 

0.050+0.027 

-25.663334 

37.746666 

Ponta  Delgada,  Azores  Islands 
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Table  A2.  WWSSN  Station  Corrections  (Continued) 


Code 

#  of  Signals 

— 

Site  Term 

Longitude 

Latitude 

Description 

PEL 

29,39,3 

0.052±0.038 

-70.685280 

-33.143612 

Peldehue,  Chile-Argentina 

PMG 

82,58,2 

-0.005±0.027 

147.153885 

-9.409166 

. 

Port  Moresby,  New  Guinea 

POO 

133,40,28 

0.037±0.023 

73.849998 

18.533333 

Pcona,  India 

PRE 

65,42,0 

-0.074±0.031 

28.190001 

-25.753334 

Pretoria,  South  Africa 

PTO 

84,52,5 

-0.140±0.027 

-8.602222 

41.138611 

Porto  Serro  Do,  Portugal 

QUE 

82,23,41 

-0.41 2±0. 026 

66.949997 

30.188334 

Quetta,  Pakistan 

QUI 

9,67,0 

0.007+0.037 

-78.500504 

-0.200139 

Quito,  Ecuador 

RAB 

45,135,0 

-0.177±0.024 

152.169830 

-4.191278 

Rabaul,  New  Britain 

RAR 

12,30,0 

-0.070+0.049 

-159.773331 

-21.212500 

Rarotonga,  Cook  Islands 

28,22,3 

0.439+0.044 

-103.208336 

44.075001 

Rapid  City,  South  Dakota 

K£S| 

9,22,0 

0.355±0.057 

151.158340 

-33.829361 

Riverview,  Australia 

SBA 

2,12,0 

-0.61 9±0. 086 

166.756104 

-77.850281 

Scott  Base,  Antarctica 

SCP 

167,68,21 

0.055±0.020 

-77.864998 

40.794998 

State  College,  Pennsylvania 

SDB 

75,17,9 

0.083±0.032 

13.571944 

-14.925834 

Sa  Da  Bandeira,  Angola 

SEO 

97,76,12 

-0.076±0.024 

126.966667 

37.566666 

Seoul  Keizyo,  South  Korea 

SHA 

76,65,0 

0.346±0.027 

-88.142807 

30.694361 

Spring  Hill,  Mississippi 

SHI 

77,14,30 

0.298±0.029 

52.519943 

29.638306 

Shiraz,  Iran 

SHK 

41,76,0 

-0.324±0.030 

132.677505 

34.532223 

Shiraki,  Honshu,  Japan 

SHL 

83,15,41 

0.033±0.027 

91.883331 

25.566668 

Shillong,  India-Bangladesh  border 

SJG 

129,57,0 

-0.248±0.023 

-66.150002 

18.111666 

San  Juan,  Puerto  Rico 

SNA 

6,11,0 

0.108±0.078 

-2.325000 

-70.315002 

Sanae,  Antarctica 

SNG 

44,31 ,3 

-0.072±0.036 

100.620003 

7.173333 

Songkhla,  Malay  Peninsula 

SPA 

13,7,0 

-0.756±0.072 

0.000000 

-90.000000 

South  Pole,  Antarctica 

STU 

172,94,20 

0.094±0.019 

9.195000 

48.771946 

Stuttgart,  Germany 

76,53,5 

0.216i0.028 

46.326668 

38.067501 

Tabriz,  Iran-USSR  border 

EH 

12,14,0 

-0.115+0.063 

147.320419 

-42.909916 

Tasmania  Univ.,  Tasmania 

TOL 

112,52,23 

0.211  ±0.023 

-4.048611 

39.881390 

Toledo,  Spain 

TRI 

128,85,25 

-0.105±0.021 

13.764167 

45.708889 

Trieste,  Italy 

TRN 

112,70,1 

0.1 01  ±0.024 

-61.402779 

10.648916 

Trinidad,  Trinidad 

TUC 

57,3,21 

0.077±0.036 

-110.782219 

32.309723 

Tucson,  Arizona 
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Table  A2.  WWSSN  Station  Corrections  (Continued) 


Code 

#  of  Signals 

Site  Term 

Longitude 

Latitude 

Description 

UME 

123,53,2 

0.170±0.024 

20.236666 

63.814999 

Umea,  Sweden 

UNM 

10,13,1 

-0.236±0.065 

-99.178085 

19.329000 

Nat.  University  of  Central  Mexico 

UPA 

2,1,0 

-0.261±0.185 

-79.533997 

8.981500 

Univ.  de  Panama,  Panama 

VAL 

122,122,12 

0.015±0.020 

-10.244166 

51.939445 

Valentia  Eire 

WEL 

8,12,0 

0.139±0.072 

174.768326 

-41.286110 

Wellington,  New  Zealand 

WES 

135,118,6 

-0.139±0.020 

-71.322083 

42.384693 

Weston,  New  England 

WIN 

32,29,0 

-0.186±0.041 

17.100000 

-22.566668 

Windhoek,  South-West  Africa 

-  73  - 


Geotech  Network  mfa 


(THIS  PAGE  INTENTIONALLY  LEFT  BLANK) 


-74- 


I'ontiji-tois  (United  Stated 


Prof  Thomas  Ahrens 
S-  :.v  illogical  Lab,  252-21 
Do  ision  of  Geological  &  Planetary  Sciences 
California  Institute  of  Technology 
Pasadena,  CA  91125 

Prof.  Charles  B.  Archambeau 
Cl  RES 

University  of  Colorado 
Boulder,  CO  80309 


Prof.  Muawia  Barazangi 

histitute  for  the  Study  of  the  Continent 

t  'orneil  University 

li'-aca.NY  14853 


Dr.  Douglas  R.  Baumgardt 
1- NS  COr  Inc 
5400  Port  Royal  Road 
Springfield,  VA  22151-2388 


Prof.  Jonathan  Berger 
1GPP,  A-025 

Scripps  Institution  of  Oceanography 
University  of  California,  San  Diego 
La  Jolla,  CA  92093 

Dr.  Lawrence  J.  Burdick 
Woodward-Clyde  Consultants 
566  El  Dorado  Street 
Pa  adena.CA  91109-3245 


u;  Karl  Coyner 

Sew  England  Research,  Inc. 

76  Olcott  Dnve 

Write  River  Junction,  VT  05001 
•  ernon  F.  Cormier 

i  Lj ;  nment  of  Geology  &  Geophysics 
U  15,  Room  207 
i  a  University  of  Connecticut 

Storrs.  CT  06268 

Professor  Anton  W.  Dainty 
E'ir'7  Resources  Laboratory 
Massachusetts  Institute  of  Technology 
42  Carieton  Street 
Cambridge,  MA  02142 

Prof.  Steven  Day 

Department  of  Geological  Sciences 
San  Diego  Slate  University 
■an  Diego,  CA  92182 


Dr.  Zoltan  A.  Der 
ENSCO,  Inc. 

5400  Port  Royal  Road 
Springfield,  VA  22151-2388 


Prof.  John  Ferguson 
Center  for  Lithospheric  Studies 
The  University  of  Texas  at  Dallas 
P.O.  Box  830688 
Richardson,  TX  75083-0688 

Prof.  Stanley  Flatte 
Applied  Sciences  Building 
University  of  California 
Santa  Cruz,  CA  95064 


Dr.  Alexander  Florence 
SRI  Lntemational 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025-3493 


Prof.  Henry  L.  Gray 
Vice  Provost  and  Dean 
Department  of  Statistical  Sciences 
Southern  Methodist  University 
Dallas,  TX  75275 

Dr.  Indra  Gupta 
Teledyne  Geotech 
314  Montgomery  Street 
Alexandria,  VA  22314 


Prof.  David  G.  Harkrider 
Seismological  Laboratory 
Division  of  Geological  &  Planetary  Sciences 
California  Institute  of  Technology 
Pasadena,  CA  91125 

Prof.  Donald  V.  Helmberger 
Seismological  Laboratory 
Division  of  Geological  &  Planetary  Sciences 
California  Institute  of  Technology 
Pasadena,  CA  91125 

Prof.  Eugene  Herrin 

Institute  for  the  Study  of  Earth  and  Man 

GeophysicalLaboratory 

Southern  Methodist  University 

Dallas,  TX  75275 

Prof  •  Robert  B.  Herrmann 

Department  of  Earth  &  Atmospheric  Science" 

St  1  .o"is  University 

St.  lxuns,  MO  63156 


- 1  - 


Prof.  Bryan  Isacks 
Cornell  University 
Department  of  Geological  Sciences 
SNEE  Hall 
Ithaca,  NY  14850 

Dr.  Rong-Song  Jih 
Teledyne  Geotech 
314  Montgomery  Street 
Alexandria,  VA  223 14 


Prof.  Lane  R.  Johnson 
Seismographic  Station 
University  of  California 
Berkeley,  CA  94720 


Prof.  Alan  Kafka 

Department  of  Geology  &  Geophysics 

Boston  College 

Chestnut  Hill,  MA  02167 


Dr.  Richard  LaCoss 
MIT-Lincoln  Laboratory 
M-200B 
P.  O.  Box  73 

Lexington,  MA  02173-0073  (3  copies) 
Prof  Fred  K.  Lamb 

University  of  Illinois  at  Urbana-Champaign 
Department  of  Physics 
1110  West  Green  Street 
Urbana,  IL  61801 

Prof.  Charles  A.  Langston 
Geosciences  Department 
403  Deike  Building 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

Prof.  Thome  Lay 

Institute  of  Tectonics 

Earth  Science  Board 

University  of  California,  Santa  Cruz 

Santa  Cruz,  CA  95064 

Prof.  Arthur  Lemer-Lam 
Lamont-Doherty  Geological  Observatory 
of  Columbia  University 
Palisades,  NY  10964 


Dr.  Christopher  Lynnes 
Teledyne  Geotech 
314  Montgomery  Street 
Alexandria,  VA  22314 


Prof.  Peter  Mai  in 

University  of  California  at  Santa  Barbara 
Institute  for  Crustal  Studies 
Santa  Barbara,  CA  93106 


Dr.  Randolph  Martin,  III 
New  England  Research,  Inc. 

76  Olcott  Drive 

White  River  Junction,  VT  05001 


Dr.  Gary  McCanor 
Mission  Research  Corporation 
735  State  Street 
P.O.  Drawer  719 

Santa  Barbara,  CA  93102  (2  copies) 

Prof.  Thomas  V.  McEvilly 
Seismographic  Station 
University  of  California 
Berkeley,  CA  94720 


Dr.  Keith  L.  McLaughlin 
S-CUBED 

A  Division  of  Maxwell  Laboratory 
P.O.  Box  1620 
La  Jolla,  CA  92038-1620 

Prof.  William  Menke 
Lamont-Doherty  Geoiogical  Observatory 
of  Columbia  University 
Palisades,  NY  10964 


Stephen  Miller 
SRI  International 
333  Ravenswood  Avenue 
Box  AF  116 

Menlo  Park,  CA  94025-3493 

Prof.  Bernard  Minster 
IGPP,  A-025 

Scripps  Institute  of  Oceanography 
University  of  California,  San  Diego 
La  Jolla,  CA  92093 

Prof.  Brian  J.  Mitchell 

Department  of  Earth  &  Atmospheric  Sciences 
St.  Louis  University 
St.  Louis,  MO  63156 


Mr.  Jack  Murphy 

S-CUBED,  A  Division  of  Maxwell  Laboratory 
1 1800  Sunrise  Valley  Drive 
Suite  1212 

Reston,  VA  22091  (2  copies) 


-2- 


Dr.  K-u-  Nguyen 
GL/LWH 

Hanscom  AFB,  MA  01731-5000 


Prof.  John  A.  Orcutt 
IGPP,  A-025 

Scripps  Institute  of  Oceanography 
University  of  California,  San  Diego 
La  Joila,  CA  92093 

Prof.  Keith  Priestley 
University  of  Cambridge 
Bullard  Labs,  Dept,  of  Earth  Sciences 
Madingley  Rise,  Madingley  Rd. 
Cambridge  CB3  OEZ,  ENGLAND 

Prof.  Paul  G.  Richards 
L-210 

Lawrence  Livermore  National  Laboratory 
Livermore,  CA  94550 


Dr.  Wilmer  Rivers 
Teledyne  Geotech 
314  Montgomery  Street 
Alexandria,  VA  22314 


Prof.  Charles  G.  Sammis 
Center  for  Earth  Sciences 
University  of  Southern  California 
University  Park 
Las  Angeles,  CA  90089-0741 

Prof.  Christopher  H.  Scholz 
Lamont- Doherty  Geological  Observatory 
of  Columbia  University 
Palisades,  NY  10964 


I'm.  David  G.  Simpson 
Lai.iont-Doherty  Geological  Observatory 
of  Columbia  University 
Palisades,  NY  10964 


Dr  Jeffrey  Stevens 
S- CUBIT) 

A  Division  of  Maxwell  Laboratory 
P.O.  Box  1620 
La  Jolla,  CA  92038-1620 

Pro f.  Brian  Stump 

Institute  for  the  Study  of  Earth  &  Man 
Geophysical  Laboratory 
Southern  Methodist  University 
Dallas,  TX  75275 


Prof.  Jeremiah  Sullivan 

University  of  Illinois  at  Urbana-Champaign 

Department  of  Physics 

1 1 10  West  Green  Street 

Urbana,  1L  61801 

Prof.  Clifford  Thurber 
University  of  Wisconsin-Madison 
Department  of  Geology  &  Geophysics 
1215  West  Dayton  Street 
Madison,  WS  53706 

Prof.  M.  Nafi  Toksoz 
Earth  Resources  Lab 
Massachusetts  Institute  of  Technology 
42  Carleton  Street 
Cambridge,  MA  02142 

Prof.  John  E.  Vidale 

University  of  California  at  Santa  Cruz 

Seismological  Laboratory 

Santa  Cruz,  CA  95064 


Prof.  Terry  C.  Wallace 
Department  of  Geosciences 
Building  #77 
University  of  Arizona 
Tucson,  AZ  85721 

Dr.  Raymond  Willeman 
GL/LWH 

Hanscom  AFB,  MA  01731-5000 


Dr.  Lorraine  Wolf 
GL/LWH 

Hanscom  AFB,  MA  01731-5000 


Prof.  Francis  T.  Wu 
Department  of  Geological  Sciences 
State  University  of  New  York 
at  Binghamton 
Vestal,  NY  13901 


3- 


OTHERS  (United  States) 


Dr.  Monem  Abdel-Gawad 
Rockwell  International  Science  Center 
1049  Camino  Dos  Rios 
Thousand  Oaks,  CA  91360 


Prof.  Keiiti  Aki 

Center  for  Earth  Sciences 

University  of  Southern  California 

University  Park 

Los  Angeles,  CA  90089-0741 

Prof.  Shelton  S.  Alexander 
Geosciences  Department 
403  Deike  Building 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

Dr.  Kenneth  Anderson 
BBNSTC 
Mail  Stop  14/IB 
Cambridge,  MA  02238 


Dr.  Ralph  Archuleta 
Department  of  Geological  Sciences 
University  of  California  at  Santa  Barbara 
Santa  Barbara,  CA  93102 


Dr.  Thomas  C.  Bache,  Jr. 
Science  Applications  Int'l  Corp. 
10210  Campus  Point  Drive 
San  Diego,  CA  92121  (2  copies) 


J.  Barker 

Department  of  Geological  Sciences 
State  University  of  New  York 
at  Binghamton 
Vestal,  NY  13901 

Dr.  T.J.  Bennett 
S-CUBED 

A  Division  of  Maxwell  Laboratory 
1 1800  Sunrise  Valley  Drive,  Suite  1212 
Reston,  VA  22091 

Mr.  William  J.  Best 
907  Westwood  Drive 
Vienna,  VA  22180 


Dr.  N.  Biswas 
Geophysical  Institute 
University  of  Alaska 
Fairbanks,  AK  99701 


Dr.  G.A.  Bollinger 
Department  of  Geological  Sciences 
Virginia  Polytechnical  Institute 
21044  Derring  Hall 
Blacksburg,  VA  24061 

Dr.  Steven  R.  Bratt 
Center  for  Seismic  Studies 
1300  North  17th  St.,  Suite  1450 
Arlington,  VA  22209 


Michael  Browne 
Teledyne  Geotech 
3401  Shiloh  Road 
Garland,  TX  75041 


Mr.  Roy  Burger 
1221  Serry  Road 
Schenectady,  NY  12309 


Dr.  Robert  Burridge 
Schlumberger-Doll  Research  Center 
Old  Quarry  Road 
Ridgefield,  CT  06877 


Dr.  Jerry  Carter 
Rondout  Associates 
P.O.  Box  224 
Stone  Ridge,  NY  12484 


Dr.  W.  Winston  Chan 
Teledyne  Geotech 
314  Montgomery  Street 
Alexandria,  VA  22314-1581 


Dr.  Theodore  Cherry 
Science  Horizons,  Inc. 

710  Encinitas  Blvd.,  Suite  200 
Encinitas,  CA  92024  (2  copies) 


Prof.  Jon  F.  Claerbout 
Department  of  Geophysics 
Stanford  University 
Stanford,  CA  94305 


Prof.  Robert  W.  Clayton 

Seismological  Laboratory 

Division  of  Geological  &  Planetary  Sciences 

California  Institute  of  Technology 

Pasadena,  CA  91125 


-4- 


Pi  or  .  f.  A.  Dahlen 

Geological  and  Geophysical  Sciences 

Princeton  University 

Princeton,  NJ  08544-0636 


Prof.  Adam  Dziewonski 
Hoffman  Laboratory 
Harvard  University 
20  Oxford  St 
Cambridge,  MA  02138 

Prof.  John  Ebel 

Department  of  Geology  &  Geophysics 
Boston  College 
Chestnut  Hill,  MA  02167 


Eric  Fielding 
SNEE  Hall 
INSTOC 

Cornell  University 
Ithaca,  NY  14853 

Prof.  Donald  Forsyth 
Department  of  Geological  Sciences 
Brown  University 
Providence,  RI  02912 


Dr.  Cliff  Frolich 
Institute  of  Geophysics 
8701  North  Mopac 
Austin,  TX  78759 


Prof.  Art  Frankel 
Mail  Stop  922 
Geological  Survey 
790  National  Center 
Ret  ton,  VA  22092 

Dr  /.nthony Gangi 
Texas  A&M  University 
Department  of  Geophysics 
College  Station,  TX  77843 


Dr  Freeman  Gilbert 

Inst,  of  Geophysics  &  Planetary  Physics 
University  of  California,  San  Diego 
P.O.  Box  109 
La  Jolla,  CA  92037 

Mr.  Edward  Giller 
Pacific  Sierra  Research  Corp. 

1401  Wilson  Boulevard 
Arlington,  VA  22209 


Dr.  Jeffrey  W.  Given 
Sierra  Geophysics 
11255  Kirkland  Way 
Kirkland,  WA  98033 


Prof.  Stephen  Grand 
University  of  Texas  at  Austin 
Department  of  Geological  Sciences 
Austin,  TX  78713-7909 


Prof.  Roy  Greenfield 
Geosciences  Department 
403  Deike  Building 
The  Pennsylvania  State  University 
University  Park,  PA  16802 

Dan  N.  Hagedorn 
Battelle 

Pacific  Northwest  Laboratories 
Battelle  Boulevard 
Richland,  WA  99352 

Kevin  Hutchenson 
Department  of  Earth  Sciences 
St.  Louis  University 
3507  Laclede 
St.  Louis,  MO  63103 

Prof.  Thomas  H.  Jordan 
Department  of  Earth,  Atmospheric 
and  Planetary  Sciences 
Massachusetts  Institute  of  Technology 
Cambridge,  MA  02 1 39 

Robert  C.  Kemerait 
ENSCO,  Inc. 

445  Pineda  Court 
Melbourne,  FL  32940 


William  Kikendall 
Teledyne  Geotech 
3401  Shiloh  Road 
Garland,  TX  75041 


Prof.  Leon  Knopoff 

University  of  California 

Institute  of  Geophysics  &  Planetary  Physics 

Los  Angeles,  CA  90024 


Prof.  L.  Timothy  Long 
School  of  Geophysical  Sciences 
Georgia  Institute  of  Technology 
Atlanta,  GA  30332 


5- 


Prof.  Art  McGarr 
Mail  Stop  977 
Geological  Survey 
345  Middlefield  RJ. 
Menlo  Park,  CA  94025 

Dr.  George  Mellman 
Sierra  Geophysics 
1 1255  Kirkland  Way 
Kirkland,  WA  98033 


Prof.  John  Nabelek 
College  of  Oceanography 
Oregon  State  University 
Corvallis,  OR  97331 


Prof.  Geza  Nagy 

University  of  California,  San  Diego 
Department  of  Ames,  M.S.  B-010 
La  Jolla,  CA  92093 


Prof.  Amos  Nur 
Department  of  Geophysics 
Stanford  University 
Stanford,  CA  94305 


Prof.  Jack  Oliver 
Department  of  Geology 
Cornell  University 
Ithaca,  NY  14850 


Prof.  Robert  Phinney 
Geological  &  Geophysical  Sciences 
Princeton  University 
Princeton,  NJ  08544-0636 


Dr.  Paul  Pomeroy 
Rondout  Associates 
P.O.  Box  224 
Stone  Ridge,  NY  12484 


Dr.  Jay  Pulli 
RADIX  System,  Iric. 

2  Taft  Court,  Suite  203 
Rockville,  MD  20850 


Dr.  Norton  Rimer 
S-CUBED 

A  Division  of  Maxwell  Laboratory 
P.O.  Box  1620 
La  Jolla,  CA  92038-1620 


Prof.  Larry  J.  Ruff 
Department  of  Geological  Sciences 
1006  C.C.  Little  Building 
University  of  Michigan 
Ann  Arbor,  MI  48109-1063 

Dr.  Richard  Sailor 
TASC  Inc. 

55  Walkers  Brook  Drive 
Reading,  MA  01867 


Thomas  J.  Sereno,  Jr. 

Science  Application  Int’l  Corp. 
10210  Campus  Point  Drive 
San  Diego,  CA  92121 


John  Sherwin 
Teledyne  Geotech 
3401  Shiloh  Road 
Garland,  TX  75041 


Prof.  Robert  Smith 
Department  of  Geophysics 
University  of  Utah 
1400  East  2nd  South 
Salt  Lake  City,  UT  84112 

Prof.  S.  W.  Smith 
Geophysics  Program 
University  of  Washington 
Seattle,  WA  98195 


Dr.  Stewart  Smith 
IRIS  Inc. 

1616  North  Fort  Myer  Drive 
Suite  1440 

Arlington,  VA  22209 

Dr.  George  Sutton 
Rondout  Associates 
P.O.  Box  224 
Stone  Ridge,  NY  12484 


Prof.  L.  Sykes 

Lamont-Doherty  Geological  Observatory 
of  Columbia  University 
Palisades,  NY  10964 


Prof.  Pradeep  Talwani 
Department  of  Geological  Sciences 
University  of  South  Carolina 
Columbia,  SC  29208 


Pro!'.  Ta-Iiang  Teng 
Center  for  Earth  Sciences 
University  of  Southern  California 
University  Park 
Los  Angeles,  CA  90089-0741 

Dr.  R.B.  Tittmann 

Rockwell  International  Science  Center 

1049  Camino  Dos  Rios 

P.O.  Box  1085 

Thousand  Oaks,  CA  91360 

Dr.  Gregory  van  der  Vink 
IRIS,  Inc. 

1616  North  Fort  Myer  Drive 
Suite  1440 

Arlington,  VA  22209 

Professor  Daniel  Walker 
University  of  Hawaii 
Institute  of  Geophysics 
Honolulu,  HI  96822 


William  R.  Walter 
Seismological  Laboratory 
University  of  Nevada 
Reno,  NV  89557 


Dr.  Gregory  Wojcik 
VVeidlinger  Associates 
4410  El  Camino  Real 
Suite  1 10 

Los  Altos,  CA  94022 

Prof.  John  H.  Woodhouse 
Hoffman  Laboratory 
Harvard  University 
20  Oxford  St. 

Cambridge,  MA  02138 

Dt.  Gregory  B.  Young 
ENSCO,  Inc. 

5400  Port  Royal  Road 
Springfield,  VA  22151-2388 


GOVERNMENT 


Dr.  Ralph  Alewine  III 
DARPA/NMRO 
1400  Wilson  Boulevard 
Arlington,  VA  22209-2308 


Mr.  James  C.  Battis 
GL/LWH 

Hanscom  AFB,  MA  01731-5000 


Dr.  Robert  Blandford 
DARPA/NMRO 
1400  Wilson  Boulevard 
Arlington,  VA  22209-2308 


Eric  Chael 
Division  924 1 
Sandia  Laboratory 
Albuquerque,  NM  87185 


Dr.  John  J.  Cipar 
GL/LWH 

Hanscoin  AFB,  MA  01731-5000 


Mr.  Jeff  Duncan 
Office  of  Congressman  Markey 
2133  Rayburn  House  Bldg. 
Washington,  DC  20515 


Dr.  Jack  Evemden 
USGS  -  Earthquake  Studies 
345  Middlefield  Road 
Menlo  Park,  CA  94025 


Art  Frankel 
USGS 

922  National  Center 
Reston,  VA  22092 


Dr.  T.  Hanks 
USGS 

Nat'I  Earthquake  Research  Center 
345  Middlefield  Road 
Menlo  Park,  CA  94025 

Dr.  James  Hannon 

Lawrence  Livermore  Nat'I  Laboratory 
P.O.  Box  808 
Livermore,  CA  94550 


Paul  Johnson 

ESS-4,  Mail  Stop  J979 

Los  Alamos  National  Laboratory 

Los  Alamos,  NM  87545 


Janet  Johnston 
GL/LWH 

Hanscom  AFB,  MA  01731-5000 


Dr.  Katharine  Kadinsky-Cade 
GL/LWH 

Hanscom  AFB,  MA  01731-5000 


Ms.  Ann  Kerr 
IGPP,  A-025 

Scripps  Institute  of  Oceanography 
University  of  California,  San  Diego 
La  Jolla,  CA  92093 

Dr.  Max  Koontz 
US  Dept  of  Energy/DP  5 
Forrestal  Building 
1000  Independence  Avenue 
Washington,  DC  20585 

Dr.  W.H.K.  Lee 

Office  of  Earthquakes,  Volcanoes, 
&  Engineering 
345  Middlefield  Road 
Menlo  Park,  CA  94025 

Dr.  William  Leith 
U.S.  Geological  Survey 
Mail  Stop  928 
Reston,  VA  22092 


Dr.  Richard  Lewis 

Director,  Earthquake  Engineering  &  Geophysic 
U.S.  Army  Corps  of  Engineers 
Box  631 

Vicksburg,  MS  39180 

James  F.  Lewkowicz 
GL/LWH 

Hanscom  AFB,  MA  01731-5000 


Mr.  Alfred  Lieberman 
ACDA/VI-OA'State  Department  Bldg 
Room  5726 
320  -  21st  Street,  NW 
Washington,  DC  20451 


Stephen  Mangino 
GL/LWH 

Hanscom  AFB,  MA  01731-5000 


Dr.  Robert  Masse 
Box  25046,  Mail  Stop  967 
Denver  Federal  Center 
Denver,  CO  80225 


Art  McGarr 

U.S.  Geological  Survey,  MS-977 
345  Middlefield  Road 
Menlo  Park,  CA  94025 


Richard  Morrow 
ACDA/VI,  Room  5741 
320  21st  Street  N.W 
Washington,  DC  20451 


Dr.  Keith  K.  Nakanishi 
Lawrence  Livermore  National  Laboratory 
P.O.  Box  808,  L-205 
Livermore,  CA  94550 


Dr.  Carl  Newton 

Los  Alamos  National  Laboratory 

P.O.  Box  1663 

Mail  Stop  C335,  Group  ESS-3 
Los  Alamos,  NM  87545 

Dr.  Kenneth  H.  Olsen 

Los  Alamos  Scientific  Laboratory 

P.O.  Box  1663 

Mail  Stop  C335,  Group  ESS-3 
Los  Alamos,  NM  87545 

Howard  J.  Patton 

Lawrence  Livermore  National  Laboratory 
P.O.  Box  808,  L-205 
Livermore,  CA  94550 


Mr.  Chris  Paine 
Office  of  Senator  Kennedy 
*  SR  315 

United  States  Senate 
Washington,  DC  20510 

Colonel  Jerry  J.  Perrizo 
AFOSR/NP,  Building  410 
Bolling  AFB 

Washington,  DC  20332-6448 


Dr.  Frank  F.  Pilotte 

HQ  AFTAQTT 

Patrick  AFB,  FL  32925-6001 


Katie  Poley 
CIA-OSWR7NED 
Washington,  DC  20505 


Mr.  Jack  Rachlin 
U.S.  Geological  Survey 
Geology,  Rm  3  Cl 36 
Mail  Stop  928  National  Center 
Reston,  VA  22092 

Dr.  Robert  Reinke 
WL/NTESG 

Kirtland  AFB,  NM  87117-6008 


Dr.  Byron  Ristvet 

HQ  DNA,  Nevada  Operations  Office 
Attn:  NVCG 
P.O.  Box  98539 
Us  Vegas,  NV  89193 

Dr.  George  Rothe 
HQ  AFTAC/TGR 
Patrick  AFB,  FL  32925-6001 


Dr.  Alan  S.  Ryall,  Jr. 
DARPA/NMRO 
1400  Wilson  Boulevard 
Arlington,  VA  22209-2308 


Dr.  Michael  Shore 
Defense  Nuclear  Agency/SPSS 
6801  Telegraph  Road 
Alexandria,  V A  22310 


Donald  L.  Springer 

Uwrence  Livermore  National  Laboratory 
P.O.  Box  808,  L-205 
Livermore ,  CA  94550 


Mr.  Charles  L.  Taylor 
G1VLWG 

Hanscom  AFB,  MA  01731-5000 


Dr.  Thomas  Weaver 
Los  Alamos  National  Laboratory 
P.O.  Box  1663,  Mail  Stop  C335 
Los  Alamos,  NM  87545 


J.J.  Zucca 

Lawrence  Livermore  National  Laboratory 
Box  808 

Livermore,  CA  94550 


GL/SULL 
Research  Library 

Hanscom  AFB  ,  MA  01731-5000  (2  copies) 


Secretary  of  the  Air  Force 
(SAFRD) 

Washington,  DC  20330 


Office  of  the  Secretary  Defense 
DDR&E 

Washington,  DC  20330 


HQ  DNA 

Attn:  Technical  Library 
Washington,  DC  20305 


DARPA/RMO/RETRffiVAL 
1400  Wilson  Boulevard 
Arlington,  VA  22209 


DARPA/RMO/Security  Office 
1400  Wilson  Boulevard 
Arlington,  VA  ^2209 


Geophysics  Laboratory 
Attn:  XO 

Hanscom  AFB,  MA  01731-5000 


Geophysics  Laboratory 
Attn:  LW 

Hanscom  AFB,  MA  01731-5000 


DARPA/PM 

1400  Wilson  Boulevard 

Arlington,  VA  22209 


Defense  Technical  Information  Center 
Cameron  Station 

Alexandria,  V A  22314  (5  copies) 


Defense  Intelligence  Agency 
Directorate  for  Scientific  & 
Technical  Intelligence  /dtib 
Washington,  IX)  20340-6158 


AFTAC/CA 

(STINFO) 

Patrick  AFB,  FL  32925-6001 


TACTEC 

Battelle  Memorial  Institute 
505  King  Avenue 

Columbus,  OH  43201  (Final  Report  Only) 


-10- 


CONTRACTORS  (Foreign) 


Dr.  Ramon  Cabre,  SJ. 
Observatorio  Sun  Calixto 
Casilla  5939 
La  Puz,  Bolivia 


Prof.  Hans-Peter  Harjes 
Institute  for  Geophysik 
Ruhr  University/Bochum 
P.O.  Box  102148 
4630  Bochum  1,  FRG 

Prof.  Eystein  Husebye 
NTNF/NORSAR 
P.O.  Box  51 

N-2007  Kjeller,  NORWAY 


Prof.  Brian  L.N.  Kennett 
Research  School  of  Earth  Sciences 
Institute  of  Advanced  Studies 
G.P.O.  Box  4 

Canberra  2601,  AUSTRALIA 

Dr.  Bernard  Massinon 

Societe  Radiomana 

27  rue  Claude  Bernard 

75005  Paris,  FRANCE  (2  Copies) 


Dr.  Pierre  Mecheler 
Societe  Radiomana 
27  rue  Claude  Bernard 
75005  Paris,  FRANCE 


Dr  Svein  Mykkeltveit 
NTNF/NORSAR 
P.O.  Box  51 

N-2007  Kjeller,  NORWAY 


FOREIGN  (Others) 


Dr.  Peter  Basham 

Earth  Physics  Branch 

Geological  Survey  of  Canada 

1  Observatory  Crescent 

Ottawa,  Ontario,  CANADA  K1A  0Y3 

Dr.  Eduard  Berg 
Institute  of  Geophysics 
University  of  Hawaii 
Honolulu,  HI  96822 


Dr.  Michel  Bouchon 
I.R.I.G.M.-B.P.  68 
38402  St.  Martin  D'Heres 
Cedex,  FRANCE 


Dr.  Hilmar  Bungum 
NTNF/NORSAR 
P.O.  Box  51 

N-2007  Kjeller,  NORWAY 


Dr.  Michel  Campillo 
Observatoire  de  Grenoble 
I.R.I.G.M.-B.P.  53 
38041  Grenoble,  FRANCE 


Dr.  Kin  Yip  Chun 
Geophysics  Division 
Physics  Department 
University  of  Toronto 
Ontario,  CANADA  M5S  1A7 

Dr.  Alan  Douglas 
Ministry  of  Defense 
Blacknest,  Brimpton 

Reading  RG7-4RS,  UNITED  KINGDOM 


Dr.  Fekadu  Kebede 
Seismological  Section 
Box  12019 

S-750  Uppsala,  SWEDEN 


Dr.  Tormod  Kvaema 
NTNF/NORSAR 
P.O.  Box  51 

N-2007  Kjeller,  NORWAY 


Dr.  Peter  Marshal 
Procurement  Executive 
Ministry  of  Defense 
Blacknest,  Brimpton 

Reading  FG7-4RS,  UNITED  KINGDOM 

Prof.  Ari  Ben-Menahem 
Department  of  Applied  Mathematics 
Weizman  Institute  of  Science 
Rehovot,  ISRAEL  951729 


Dr.  Robert  North 

Geophysics  Division 

Geological  Survey  of  Canada 

1  Observatory  Crescent 

Ottawa,  Ontario,  CANADA  K 1 A  0Y3 

Dr.  Frode  Ringdal 
NTNF/NORSAR 
P.O.  Box  51 

N-2007  Kjeller,  NORWAY 


Dr.  Jorg  Schlittenhardt 

Federal  Institute  for  Geosciences  &  Nat'l  Res. 
Postfach  510153 

D-3000  Hannover  5 1 .  FEDERAL  REPUBLIC  OF 
GERMANY 


Dr.  Roger  Hansen 
NTNF/NORSAR 
P.O.  Box  51 

N-2007  Kjeller,  NORWAY 


Dr.  Manfred  Henger 

Federal  Institute  for  Geosciences  &  Nat'l  Res. 
Postfach  510153 
D-3000  Hanover  51,  FRG 


Ms.  Eva  Johannisson 
Senior  Research  Officer 
National  Defense  Research  Inst. 
P.O.  Box  27322 
S- 102  54  Stockholm,  SWEDEN 


